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PBEFACE. 


IN  March,  1875,  the  usual  biennial  notice  was  issued,  giving  the 
subjects  for  the  Adams  Prize  to  be  adjudged  in  1877.  The 
following  is  the  chief  portion  of  the  notice  : 

The  University  having  accepted  a  Fund  raised  by  several  members 
of  St  John's  College  for  the  purpose  of  founding  a  Prize  to  be  called  the 
Adams  Prize,  for  the  best  essay  on  some  subject  of  Pure  Mathematics, 
Astronomy  or  other  branch  of  Natural  Philosophy,  the  Prize* to  be 
given  once  in  two  years,  and  to  be  open  to  the  competition  of  all 
persons  who  have  at  any  time  been  admitted  to  a  degree  in  this 
University — 

The  Examiners  give  notice  that  the  following  is  the  subject  of  the 
Prize  to  be  adjudged  in  1877  :  The  Criterion  of  Dynamical  Stability. 

To  illustrate  the  meaning  of  the  question  imagine  a  particle  to  slide 
down  inside  a  smooth  inclined  cylinder  along  the  lowest  generating  line, 
or  to  slide  down  outside  along  the  highest  generating  line.  In  the 
former  case  a  slight  derangement  of  the  motion  would  merely  cause 
the  particle  to  oscillate  about  the  generating  line,  while  in  the  latter 
case  the  particle  would  depart  from  the  generating  line  altogether. 
The  motion  in  the  former  case  would  be,  in  the  sense  of  the  question, 
stable,  in  the  latter  unstable. 

The   criterion  of  the  stability  of  the  equilibrium  of  a  system  is, 
that  its  potential  energy  should  be  a  minimum;  what  is  desired  is,  a 
R.  A.  b 
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corresponding  condition  enabling  us  to  decide  when  a  dynamically  pos- 
sible motion  of  a  system  is  such,  that  if  slightly  deranged  the  motion 
shall  continue  to  be  only  slightly  departed  from. 

The  essays  must  be  sent  in  to  the  Vice-Chancellor  on  or  before  the 
16th  December  1876,  &c.,  &c. 

S.  G.  PHEAR,  Vice-Chancellor. 

J.  CHALLIS. 

G.  G.  STOKES. 

J.  CLERK  MAXWELL. 


The  pressure  of  other  engagements  for  some  time  prevented 
me  from  giving  my  attention  to  the  subject.  This  essay  was 
therefore  almost  entirely  composed  during  the  year  1876.  It  is 
now  printed  as  it  was  sent  in  to  the  Examiners,  the  changes  being 
merely  verbal.  Some  few  additions  have  been  made  where  ex- 
planation appeared  to  be  necessary,  but  all  these  have  been 
marked  by  square  brackets,  so  that  they  can  be  at  once  dis- 
tinguished from  the  original  parts  of  the  essay. 

In  order  to  shorten  the  essay  as  much  as  possible  many 
merely  algebraic  processes  have  been  omitted  and  the  results  only 
are  stated.  It  is  hoped  that  this  will  add  clearness  as  well  as 
brevity  to  the  reasoning,  as  the  attention  of  the  reader  will  not 
be  called  from  the  argument  to  follow  a  manipulation  of  symbols 
which  may  not  present  any  novelty. 

The  line  of  argument  taken  may  be  indicated  in  a  general  way 
as  follows.  Chapter  I.  begins  with  some  definitions  of  the  terms 
stable  and  steady  motions.  It  is  then  pointed  out  that  whether 
the  forces  which  act  on  the  system  admit  of  a  force-function  or  not, 
the  stability  of  the  motion,  if  steady,  is  indicated  by  the  nature 
of  the  roots  of  a  certain  determinantal  equation.  The  boundary 
between  stability  and  instability  being  generally  indicated  by  the 
presence  of  equal  roots,  a  criterion  is  investigated  to  determine 
beforehand  whether  equal  roots  do  or  do  not  imply  instability. 
This  case  being  disposed  of,  the  consideration  of  the  determinantal 
equation  is  resumed.  Two  general  methods  are  given  by  which, 
without  solving  the  equation,  it  may  be  ascertained  whether  the 
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character  of  the  roots  imply  stability  or  instability.  These  occupy 
Chapters  II.  and  in.  In  the  first  method  a  derived  equation  is 
made  use  of,  and  it  is  shown  that  a  simple  inspection  of  the  signs 
of  the  coefficients  of  the  several  powers  in  these  two  equations  will 
decide  the  question  of  stability.  In  the  second  method  a  certain 
easy  process  is  found  which  if  performed  on  the  determinantal 
equation  will  lead  to  the  criteria  of  stability.  At  the  end  of  the 
third  Chapter  a  geometrical  interpretation  is  given  to  the  argu- 
ment. 

In  the  fourth  Chapter  the  forces  which  act  on  the  system 
are  supposed  to  have  a  force-function.  The  determinantal  equa- 
tion is  then  much  simplified.  Several  points  are  considered  in 
this  Chapter  which  are  necessary  to  the  argument,  such  as  the 
proper  method  of  choosing  the  steady  co-ordinates  (if  there  be 
any),  the  distinction  between  harmonic  oscillation  about  steady 
motion  and  that  about  equilibrium,  and  the  changes  which  must 
be  made  in  the  determinantal  equation  when  the  equations  of 
Lagrange  become  inapplicable.  A  method  of  modifying  the 
Lagrangian  function  is  also  given  by  which,  in  certain  cases, 
the  fundamental  determinant  may  be  reduced  to  one  of  fewer 
rows  and  columns. 

In  the  fifth  Chapter  a  series  of  subsidiary  determinants  is 
formed,  and  it  is  shown  that  at  least  as  many  of  the  conditions 
of  stability  are  satisfied  as  there  are  variations  of  signs  lost  in  the 
series  in  passing  from  one  given  state. to  another.  It  is  also  shown 
that  this  is  equivalent  to  a  maximum  condition  of  the  Lagrangian 
function. 

In  the  sixth  Chapter  the  energy  test  of  stability  is  considered. 
It  is  also  shown  that,  when  the  motion  is  steady,  this  reduces 
to  the  same  criterion  as  that  indicated  in  Chap.  V. 

In  the  seventh  Chapter  the  question  considered  is  whether 
the  stability  of  a  state  of  motion  can  really  be  determined  by  an 
examination  of  the  terms  of  the  first  order  only.  In  some  cases 
these  are  certainly  sufficient,  and  an  attempt  is  made  to  discrimi- 
nate between  these  cases  and  those  in  which  the  terms  of  the 
higher  order  ultimately  alter  the  character  of  the  motion. 

If  the  Hamiltonian  characteristic  and  principal  functions  be 
given,  the  conditions  of  stability  as  regards  space  only,  or  both 
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space  and  time  may  be  deduced.  But  if  these  be  not  known  as 
expressed  in  the  Hamiltonian  form,  we  may  yet  sometimes  dis- 
tinguish between  stability  and  instability  if  we  can  determine 
whether  a  certain  integral  ceases  to  be  a  minimum  at  some  instant 
of  tne  motion.  This  is  the  subject  of  the  eighth  chapter. 

As  part  of  the  third  edition  of  my  treatise  on  the  Dynamics  of 
Eigid  Bodies  was  written  at  the  same  time  as  this  essay,  there  are 
necessarily  points  of  contact  between  the  two  works.  Thus  the 
subjects  of  the  first  part  of  the  seventh  chapter  and  of  a  portion 
of  the  sixth  will  be  found  discussed  in  the  treatise  on  Dynamics. 
But  as  tne  objects  of  the  two  books  are  not  the  same,  it  will  be 
found  that  in  all  these  cases  there  are  considerable  differences  in 
the  modes  of  demonstration. 


EDWARD  J,  ROTJTH, 


PETEKHOUSE, 

August  14, 1877. 
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CHAPTER  I. 


Definitions  of  the  terms  small  quantity,  stable  motion,  steady  motion. 
Arts.  1,  2. 

A  system  of  bodies  in  steady  motion  is  stable  if  the  roots  of  a  certain 
determinantal  equation  are  such  tJiat  their  real  parts  are  all  negative. 
Art.  3. 

Effect  of  equal  roots,  and  a  test  to  determine  whether  equal  roots  do  or  do 
not  introduce  terms  which  contain  the  time  as  afcwtor.  Arts.  4 — 7. 

Object  of  Chapters  II.  and  III.     Art,  8. 

1.  Let  us  suppose  a  dynamical  system  to  be  set  in  motion 
under  any  forces  and  to  move  in  some  known  manner.  If  any 
small  disturbance  be  given  to  the  system,  it  may  deviate  only 
slightly  from  its  known  motion,  or  it  may  diverge  further  and 
further  from  it.  Let  6,  $,  &c.  be  the  independent  variables  or 
co-ordinates  which  determine  the  position  of  the  system,  and  let 
the  known  motion  be  given  by  0  =  #0,  (f>  =  </>0,  &c.  where  00,  <£0,  &c. 
are  known  functions  of  the  time  t  To  discover  the  disturbance  of 
the  system  we  put  0  =  00  +  xy  $  =  <£0  -f-  ?/,  &c.  These  quantities 
x,  y,  &c.  are  in  the  first  instance  very  small  because  the  disturb- 
ance is  small.  The  quantities  x,  y,  z,  &c.  are  said  to  be  small 
when  it  is  possible  to  choose  some  quantity  numerically  greater 
than  all  of  them,  which  is  such  that  its  square  can  be  neglected. 
This  quantity  may  be  called  the  standard  of  reference  for  small 
quantities. 

If,  after  the  disturbance,  the  co-ordinates  x,  y,  z,  &c.  remain 
always  small,  the  undisturbed  motion  is  said  to  be  stable;  if,  on  the 
other  hand,  anyone  of  the  co-ordinates  become  large,  the  motion  is 
called  unstable. 

It  is  clear  that  the  same  motion  may  be  stable  for  one  kind  of 
disturbance  and  unstable  for  another.  But  it  is  usual  to  suppose 
the  disturbance  general,  so  that  if  the  motion  can  be  made  un- 
stable by  any  kind  of  disturbance  (provided  it  be  small)  it  is  said 
to  be  unstable.  On  the  other  hand,  it  will  be  called  stable  only 
when  it  is  stable  for  all  kinds  of  small  disturbances. 

R.  A.  1 


2  DEFINITIONS.  [CHAP. 

2.  To  determine  whether  x,  y,  z,  &c.  remain  small,  we  must 
substitute  for  0,  </>,  &c.  in  the  equations  of  motion  their  values 
#0  +  #,  fa  +  y,  &c.  Assuming  that  #,  y,  &c.  remain  small,  we  may 
neglect  their  squares,  and  thus  the  resulting  equations  will  be 

linear  in  x,  y,  z,  &c.    The  coefficients  of  x,  -®-  ,  -J£  ,  y,  -~  ,  -~,  &c. 

in  these  equations  may  be  either  constants  or  functions  of  the 
time.  In  the  former  case  the  undisturbed  motion  is  said  to  be 
steady  for  these  co-ordinates,  in  the  latter  unsteady.  In  the  case 
of  a  steady  motion  x,  y,  z,  &c.  are  all  functions  of  the  time 
which  has  elapsed  since  the  disturbance  and  of  certain  constants 
of  integration  which  are  determined  by  the  initial  values  of 

fj  w          //  7/ 

x,  -T-.,yi-jf>  &c.     We  may  therefore  define  a  steady  motion  to  be 
cLv         dt 

such  that  the  same  change  of  motion  follows  from  the  same  initial 
disturbance  at  whatever  instant  the  disturbance  is  communicated 
to  the  system. 

If  all  the  coefficients  in  the  equations  to  find  x,  T/,  z  are  con- 
stant, they  may  be  made  to  contain  t  by  a  change  of  co-ordinates. 
Thus  we  may  write  for  x,  y,  z,  &c. 


where  a,  /3,  &c.  are  any  functions  of  t  we  please.  Conversely,  when 
the  coefficients  are  functions  of  t,  we  may  sometimes  make  the 
coefficients  constant  by  a  proper  change  of  co-ordinates.  But  this 
cannot  always  be  done.  If  there  are  n  co-ordinates,  we  have  n* 
arbitrary  functions  a,  /3,  &c.  at  our  disposal.  In  each  of  the  n 
linear  equations  of  motion  we  may  have  three  terms  for  each 
co-ordinate,  and  thus  we  have  (3  n  —  1)  n  coefficients  to  make  con- 
stants. We  have  therefore  in  general  too  many  equations  to 
satisfy.  The  proper  method  of  choosing  the  co-ordinates  of  refer- 
ence will  be  considered  in  a  future  chapter. 

3.  Let  us  suppose  a  dynamical  system  to  be  making  small 
oscillations  under  the  action  of  any  forces  which  may,  or  may  not, 
possess  a  force  function  and  to  be  subject  to  any  resistances  which 
vary  as  the  velocities  of  the  parts  resisted.  The  general  equations 
of  motion  will  then  be  of  the  form 


&c.  = 


I.]  THE   DETERMINANTAL   EQUATION.  3 

To  solve  these  equations  we  write 

Substituting  we  obtain  a  determinantal  equation  to  find  m. 
If  we  put 

A  =  Azmz  +  Ajn  +  AQ,        B  =  Bzmz  +  BJM  +  BQ, 
A'  =  A^m*  +  A^m  +  A0't  &c. 

this  equation  may  be  written  in  the  simple  form 


A,   B,    C. 
A',  B',  C'. 


We  may  also  write  the  equation  in  the  form  /  (m)  =  0. 

The  coefficients  M,  M',  &c.  are  not  independent,  but  if  we 
represent  the  minors  of  A,  B,  C,  &c.  by  a,  b,  c,  &c.  we  may  easily 
show  that 

MssM  =  ]Mr_ 

a       b         c 
We  also  have 

M    M     M"     , 

—  =-TT  =  -7-=&C. 

a       b        c 

It  may  be  shown  by  properties  of  determinants  that  these 
equations  all  give  the  same  ratios.  If  A2  be  the  second  minor 
obtained  from  the  determinant  f(m)  by  omitting  the  first  and 
second  rows  and  columns,  we  know  that 

A,/  (m)  =  ab'  -  a'b. 

7 

Hence  iff(m)  =  0  we  have  -=r,. 

a      b 

In  the  same  way  we  may  show  that  —,  =  -,,  and  so  on.     This 

a      c 

property  of  Determinants  is  given  in  Dr  Salmon's  Higher  Algebra, 
Lesson  iv.  Ex.  1. 

The  general  solution  of  the  equation  may  therefore  be  written 
in  the  form 


x  = 


where  Lv  L2...  are  arbitrary  constants,  av  at,  &c.  the  values  of 
the  minor  a  when  mv  m2...  are  substituted  for  m  •  bv  lz  ...  the 

1—2 
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values  of  the  minor  b  when  similar  substitutions  are  made,  and 
so  on. 

4.  We  see  that  the  whole  character  of  the  motion  will  de- 
pend on  the  signs  of  the  quantities  mv  mz...  If  any  one  be  real 
and  positive,  x,  y,  &c.  or  some  of  them  will  ultimately  become 
large,  and  the  steady  motion  about  which  the  system  is  oscillating 
will  be  unstable.  If  all  the  roots  are  real,  negative  or  zero  and 
unequal,  the  motion  will  be  stable. 

If  two  of  the  roots  be  imaginary  we  have  a  pair  of  imaginary 
exponentials.  If  these  imaginary  roots  be  a  +  J3*J  —  1,  the  terms 
can  be  rationalized  into 

eat(Nl  cos  fit  +  N2  sin  j3t). 

The  motion  will  be  stable  if  a  be  negative  or  zero,  and  unstable 
if  a  be  positive. 

If  two  roots  be  equal,  the  form  of  the  solution  is  changed.  Let 
m^  =  ml  +  h  where  h  will  be  ultimately  zero,  we  then  have 

x  =  L^el"1**  +  Lz  (ajF*  +  ~l 

If  we  now  make  L^  and  L2  infinite  in  the  usual  manner,  we 
find 


&c.  =  &c., 

where  Mv  Mz  are  two  arbitrary  constants  which  replace  Llt 
In  the  same  way  if  three  roots  are  equal  we  have 


This  rule  will  be  found  convenient  in  practice  to  supply  the 
defect  in  the  number  of  arbitrary  constants  produced  by  equal 
roots.  At  present  we  are  only  concerned  with  their  effect  on  the 
stability  of  the  system.  The  terms  which  contain  t  as  a  factor 
will  at  first  increase  with  t,  but  if  m  be  negative,  the  term  tnemt 

can  never  be  numerically  greater  than  -  -  .     If  m  be  very  small 

C///C'' 
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the  initial  increase  of  the  terms  may  make  the  values  of  x  and  y 
become  large,  and  the  motion  cannot  be  regarded  as  a  small 
oscillation.  But  if  the  system  be  not  so  much  disturbed  that 

n 

M.-       is  large,  the  terms  will  ultimately   disappear   and   the 

6rfv 

motion  may  be  regarded  as  stable.  If,  however,  the  real  parts  of 
the  equal  roots  are  positive  or  zero,  the  terms  will  become  large 
and  the  motion  will  be  unstable. 

5.  In  some  cases,  however,  the  relations  which  exist  between 
the  coefficients  are  such  that  the  terms  which  contain  t  as  a  factor 
are  all  zero.  It  is  of  some  importance  to  discriminate  these  cases, 
for  the  stability  of  the  system  is  then  unaffected  by  the  presence 
of  equal  roots. 

Let  us  suppose  first  that  the  determinantal  equation  has  two 
roots  only  equal  to  mv  and  let  the  terms  depending  on  these  be 


&c.  =  &c. 

Substituting  in  the  equations  of  Art.  (3)  we  have,  following 
the  same  notation  as  before, 


.IT. 


&c.  =  &c. 

To  avoid  entering  more  minutely  than  is  necessary  into  the 
properties  of  linear  equations,  we  shall  assume  that  these  equations 
for  the  given  value  of  m  lead  to  but  one  solution  with  two  of  the 
N's  arbitrary,  unless  the  determinantal  equation  has  more  than 
two  roots  equal  to  mr  If  in  this  unique  solution  the  jV2's  are  all 
zero  we  must  have 


=  o    in. 

=  0 
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Since  two  of  the  constants  N.,  N'  &c.  are  to  be  arbitrary,  let 

N     N' 

them  be  N^  NJ,  then  since  --1  =  -±  we  must  have  the  minors  a 

a        o 

N     N" 

and  b  each  equal  to  zero.     Also  since  —  -  =  ^-J-,  we  shall  have  N'' 

ct>        c 

infinite  unless  c  =  0.     In  the  same  way  we  may  prove  that  all  the 
other  first  minors  are  zero. 

And  if  the  first  minors  are  zero,  we  may  show  that  two  of  the 
equations   may  be   deduced   from    the   others.     Let  the   symbol 

Atjy\  represent  the  second  minor,  with  the  usual  sign,  formed 
by  omitting  the  rows  and  columns  in  which  AB  A'B'  occur. 
Then  since  the  minors  a,  b,  c,  &c.  are  zero,  we  have 


IV. 


VAB^ 
[_A'B'\ 


B''+B        ""+'"  = 

&c.  =  OJ 


Omitting  the  first  line  of  III.  let  us  multiply  the  others  by 
A'T¥\     A"ft"\  ^C>>  resPectively.     Adding  the  results,  we   have 

an  identity.  Hence  the  second  equation  may  be  deduced  from 
the  others  which  follow  it.  In  the  same  way,  the  first  equation 
may  be  deduced  from  the  others. 

"Rejecting  the  first  two  equations,  let  us  transpose  the  arbitrary 
constants  N^  and  N{  to  the  right-hand  sides  of  the  remaining 
equations.  If  there  are  to  be  only  two  arbitrary  constants,  these 
remaining  equations  must  be  independent;  solving,  we  have 


*  _B>G>     *    'tfj 

with   similar   equations   for    the   others.      Hence   the   constants 
Nv  NI,  &c.,  are  connected  by  equations  of  the  form 

c-\    .7/,[AB~\ 

d*  •  U*J"* 

so  that  when  any  two  are  chosen  as  the  arbitrary  ones,  the  others 
may  be  deduced  from  them. 

If  the  determinantal  equation  has  three  roots  equal  to  mlt  and 
if  the  terms  which  contain  t  as  a  factor  are  all  zero,  the  equa- 
tions III.  must  admit  of  a  solution  with  three  of  the  constants 
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Nlt  NJ,  &c.  arbitrary.     If  these  be  tflt  N^  N",  we  see  that  the 

r  BG~\  r  AC~\VAB^\ 
second  minors        '  ''    m  Zer°' 


BG~\  r  AC~\VAB^\ 
Lf<7      ^'i?' 


Next  since 


*  _B'D' 


'-S  J 


we  must  have  and      ,,,,    both  zero  or  .N"  infinite. 


Hence   all  the  second  minors  are  zero.     And  if  the  second 
minors  are  all  zero,  we  have 


A" 


"  ABC  ' 
AB'C' 

+  A" 

"    ABC 
AB'C' 

...  =  o, 


and  by  similar  reasoning  three  equations  may  be  deduced  from 
the  remaining  ones.     We  have  then 


"  BCD  " 

~W 

"  ACD 

+N? 

'  ABD 

-N- 

"  ABC  ' 

B'C'D' 

AC'D' 

ABD' 

AB'C' 

JB'C"D\ 

_A"  C"D"_ 

_A"B"D\ 

_A"B"C", 

with  similar  equations. 

Since  f(m)  is  the  determinant  formed  by  eliminating  Nlt  N^, 
&c.  from  III.  we  have 


dm 


~*~ 


dA   dm~*~    dB    dm 


dB  df(m)  dA' 

*         '        '          ~* 


dA' 


dA 

-r- 
dm 


-r- 
dm 


,dA' 
-j— 

dm 


This  vanishes  when  a,  b,  &c.,  a,  &c.  are  all  zero.     If  therefore  the 
first  minors  of/(m)  all  vanish  when  m  =  m15  the  equation  /(m)  =  0 

has  two  roots  equal  to  mv     In  the  same  way  -=-  vanishes  if  all  its 
first  minors  are  zero.     But 


dm* 


da  dA      0 
a  -T-J  +  -j—  -j-  +  &c. 
dm      dm  am 


vanishes  if  «,  -  —    &c.  are  all  zero.     If  therefore  the  first  and 
dm 

second   minors   of  f(m)  all  vanish   when  m  =  ml>  the  equation 


8  EFFECT  OF  EQUAL  ROOTS  [CHAP. 

f(m)  =  0  has  three  roots  equal  to  ???v     It  is  evident  the  proposi- 
tion may  be  extended  to  any  number  of  roots. 

The  test  that,  when  equal  roots  occur  in  the  determinantal  equa- 
tion, the  terms  in  the  values  of  x,  y,  &c.  which  contain  t  as  a  factor 
should  be  absent  may  be  stated  thus.  If  there  are  two  equal  roots 
all  the  first  minors  must  vanish.  If  three  equal  roots,  all  the  first 
and  second  minors  must  vanish,  and  so  on.  In  these  cases  the 
equal  roots  introduce  merely  a  corresponding  indeterminateness 
into  the  coefficients. 

When  there  are  more  equal  roots  than  there  are  rows  in  the 
determinantal  equation,  it  is  easy  to  see  that  there  must  be  some 
terms  in  the  integrals  which  contain  t  as  a  factor. 

[The  following  simple  example  will  illustrate  the  application 
of  this  test. 

A  particle  is  in  equilibrium  at  the  origin  of  co-ordinates  under 
the  action  of  forces  whose  force  function  U  is  given  by 


U  =  J  Ax*  +  J  By*  +  i  C£  +  Dyz  +  Ezx  +  Fxy. 

If  the  level  surfaces  are  ellipsoids  and  the  force  acts  inwards, 
it  is  clear  that  the  equilibrium  of  the  particle  must  always  be 
stable.  If  then  any  equal  roots  occur  in  the  determinantal 
equation,  the  test  should  show  that  the  terms  which  contain  t  as 
a  factor  are  absent. 

If  T  be  the  semi  vis  viva  of  the  particle  and  if  its  mass  be 
taken  as  unity,  we  have 


Omitting  accents  and  forming  the  discriminant  of  —m2T+  U 
we  have  the  following  determinantal  equation  : 


A-m2  F  E 
F  B-m  D 
E  D  G-mz 


=  0. 


This  is  the  "discriminating  cubic"  which  determines  the  axes 
of  the  quadric  U=  c,  where  c  is  a  constant.  The  conditions  that 
two  of  its  roots  should  be  equal,  i.e.  that  the  quadric  should  be 
a  spheroid,  are  well  known  to  be 

EF  FD  DE 


where  m*  is  equal  to  either  root.     These  are  just  the  conditions 
obtained  by  equating  any  first  minor  of  the  determinant  to  zero. 
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The  conditions  that  three  of  the  roots  of  the  cubic  should  be 
equal,  i.e.  that  the  quadric  should  be  a  sphere,  are 

A  =  B=C,   Z>  =  0,    E=0,    F=0. 
These  are  the  conditions  that  every  second  minor  should  vanish. 

In  this  example  we  have  taken  the  case  of  a  single  particle. 
Similar  remarks  however  apply,  when  any  system  of  bodies  is 
disturbed  from  a  state  of  stable  equilibrium.  The  oscillations 
may  be  found  by  the  method  of  Lagrange.  The  final  determi- 
nantal  equation  may  be  conveniently  formed  by  equating  to  zero 
the  discriminant  of  —  m?T+  U,  where  T  is  the  semi  vis  viva  with 
the  accents  denoting  differentiations  with  regard  to  the  time 
omitted,  and  U  is  the  force  function.  It  is  a  known  theorem  that 
the  existence  of  finite  equal  roots  does  not  affect  the  stability  of 
the  equilibrium.  Hence  the  conditions  for  equal  roots  must  be 
such  as  to  make  all  the  minors  equal  to  zero.  Conversely,  this 
theorem  will  often  conveniently  give  the  conditions  that  Lagrange's 
determinant  has  equal  roots.] 

6.  That  there  should  be  a  difference  in  the  modes  in  which 
equal  roots  affect  the  motion  is  no  more  than  we  should  expect 
a  priori.  Suppose  the  coefficients  of  the  equation  /  (m)  =  0  to  be 
functions  of  some  quantity  n,  and  that  as  n  passes  through  the 
value  nQi  two  roots  become  equal  to  each  other.  Let  the  quadratic 
factor  containing  these  roots  be  m2  +  2a/7i  +  ft  and  let  us  consider 
only  the  case  in  which  a  and  (3  are  real.  We  have  a2  —  ft  =  0 
when  n  =  nQ .  If  a2  —  /3  change  sign  as  n  passes  through  the  value 
n0,  the  roots  will  change  from  a  trigonometrical  to  a  purely  ex- 
ponential form,  which  would  indicate  a  change  from  oscillatory  to 
non- oscillatory  motion.  The  passage  from  one  kind  of  motion  to 
the  other  may  be  effected  through  a  motion  represented  by  expres- 
sions having  the  time  as  a  factor.  But  if  a2  —  /3  does  not  change 
sign,  for  example,  if  it  be  a  perfect  square  for  all  values  of  n,  there 
will  be  no  change  from  one  kind  of  motion  to  the  other,  and  in  this 
case  we  should  expect  that  the  motion  when  the  roots  are  equal 
will  be  represented  by  terms  of  the  same  character  as  before. 
Briefly,  we  may  expect  equal  roots  to  introduce  terms  with  t  as  a 
factor  at  the  boundary  between  stability  and  instability ;  and  to 
introduce  merely  an  indeterminateness  into  the  coefficients  when 
the  motion  is  stable  on  both  sides. 

It  is  easy  to  show  that  in  the  first  of  these  two  cases  the 
minors  could  not  contain  either  of  the  factors  of  m2  +  2aw  +  ft 
For  since  a2  —  0  changes  sign,  these  factors  are  in  one  case  ima- 
ginary ;  and  therefore  if  one  factor  occur  in  any  minor  the  other 
must  also  be  present.  The  minors  would  not  only  vanish,  but 
must  have  equal  roots  also.  But  as  in  Art.  (3),  A2/(???)  =ab'  —db. 
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Hence  if  all  the  first  minors  have  equal  roots  it  is  clear  that  either 
f(m)  has  more  than  two  equal  roots,  or  all  the  second  minors  must 
vanish.  The  latter  is  impossible  unless  f(m)  has  more  than  two 
equal  roots. 

These  general  considerations  are  not  meant  to  replace  the 
proofs  given  in  the  last  article,  but  merely  to  explain  how  a 
difference  in  the  effects  of  the  equal  roots  might  arise. 

7.  Summing  up  what  precedes,  we  see  that  if  a  dynamical 
system  have  n  co-ordinates  its  stability  depends  on  the  nature  of 
the  roots  of  a  certain  equation  of  the  2nih  degree. 

If  the  roots  of  this  equation  are  all  unequal,  the  motion  will 
be  stable  if  the  real  roots  and  the  real  parts  of  the  imaginary  roots 
are  all  negative  or  zero,  and  unstable  if  any  one  is  positive.  If 
several  roots  are  equal  the  motion  will  be  stable  if  the  real  parts 
of  those  roots  are  negative  and  not  very  small,  and  unstable  if 
they  are  negative  and  small,  zero,  or  any  positive  quantity.  But 
if,  as  often  happens  in  dynamical  problems,  the  terms  which  con- 
tain t  as  a  factor  are  absent  from  the  solution,  the  condition  of 
stability  is  that  the  real  roots  and  the  real  parts  of  the  imaginary 
roots  of  the  subsidiary  equation  should  be  negative  or  zero. 

8.  When  the  equation  f  (D)  —  0  is  of  low  dimensions  we  may 
solve   it   or   otherwise   determine   the   nature   of  its   roots ;   the 
stability  or   instability  of  the  system  will  then  become  known. 
But  if  the  degree  of  the  equation  be  considerable  this  is  not  a  very 
easy  problem.     We  shall  devote  the  two  next  Chapters  to  the 
consideration  of  two  methods  by  either  of  which,  without  solving 
the  equation,  we  can  determine  the  conditions  that  the  real  roots 
and  the  real  parts  of  the  imaginary  roots  should  be  all  negative. 
The  determination  of  these  conditions  has,  it  appears,  never  before 
been  accomplished.*    The  consideration  of  the  equations  of  motion 
will  then  be  resumed,  and  the  form  of  the  determinantal  equation 

/  (Z>)  =  0  when  the  forces  admit  of  a  force  function  will  be  more 
particularly  investigated. 

*  [These  conditions  for  the  cases  of  a  biquadratic  and  a  qnintic  had  been  found 
by  the  Author  in  1873,  and  read  before  the  London  Mathematical  Society  in  June, 
1871  See  also  the  third  edition  of  the  Author's  Rigid  Dynamics,  Art.  436.] 


CHAPTER  II. 


Uatement  of  the  theory  by  which  the  necessary  and  sufficient  tests  of 
stability  are  found.     Objections  to  this  theory.     Arts.  1 — 3. 

These  tests  shown  to  be  integral  functions  of  the  coefficients.     Art.  4. 

Method   of  finding  these  tests  when  the  coefficients  of  the  equation  are 
numerical,  or  when  several  terms  are  absent.     Arts.  5,  6. 

All  these  tests  shown  to  be  derivable  from  one  called  the  fundamental 
term.     Arts.  7,  8. 

The  fundamental  term  found  as  an  eliminant.     Art.  9. 

A  method  of  finding  the  fundamental  term  by  derivation  from  the  funda- 
mental term  of  one  degree  lower.     Art.  10. 

Another  and  better  method  of  doing  the  same  by  means  of  a  differential 
equation.     Arts.  11 — 13. 


1.  The  object  of  this  Chapter  has  been  explained  at  the  end 
of  Chapter  I.  Briefly,  the  criterion  that  the  motion  of  a  system 
of  bodies  should  be  stable  is  that  the  roots  of  a  certain  equation 
should  have  all  their  real  parts  negative.  We  propose  to  investi- 
gate these  conditions. 

Let  the  equation  to  be  considered  be 

f(x)=pjf+plarl+...  +*V.1«  +  p»=  0. 
Let  the  real  roots  be  alt  a2...  and  the  imaginary  roots  be 
aiV-l,    a±/3V^l,  &c. 


*),  &c. 


Then 


If  then  alf  aa,  &c.  at,  ora,  &c.  are  all  negative,  every  term  in 
each  factor,  and  therefore  in  the  product,  must  be  positive. 

It  is  therefore  necessary  that  every  term  in  the  equation 
f(x)  =  Q  should  have  the  same  sign.  It  will  be  convenient  to 
suppose  this  sign  to  be  positive. 
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It  is  also  clear  on  the  same  suppositions  that  none  of  the 
coefficients^,^,  ...pn  can  be  zero,  except  when  the  roots  of  the 
equation  are  all  of  the  form  +  /3\/—  1,  or  when  some  of  the  roots 
are  zero. 

2.  Let  us  now  form  the  equation  whose  roots  are  the  sums  of 
the  roots  of  #  taken  two  and  two.  Let  this  be 


F  (x)  =  PQxm  +  PX1'1  +  •  •  •  +  P^x  +  Pm  =  0, 


n  —  1 

where  m  —  n  —  0—  .      The   real    roots   of   this   equation   will  be 

2i 

ttj  +  ttg,  a:  +  a3,  &c.  2aj,  2ar2,  &c.  and  the  imaginary  roots  will  be 
at  4-  ffj  +  /8j  V  —  1,  &c.  It  is  clear  from  the  same  reasoning  as 
before  that  if  at,a8,&c.  or,,  or2,  &c.  are  all  negative,  the  coefficients 
P0,  Pt,  &c.  must  all  have  the  same  sign. 

Conversely,  if  p0,  pl.-.  have  all  the  same  sign,  the  equation 
f(x)  can  have  no  real  positive  root,  and  if  P0,  Pl  ...  Pm  have  all 
the  same  sign  the  equation  F(x)  can  have  no  positive  root,  and 
therefore  f(x)  can  have  no  imaginary  root  with  its  real  part 
positive. 

3.  Our  first  test  of  the  stability  of  a  dynamical  system  is  that 
all  the  coefficients  of  the  dynamical  equation  f  (D)  =  0  and  all  the 
coefficients  of  its  derived  equation  F  (D)  =  0  should  have  the  same 
sign. 

It  should  be  noticed  that  though  these  conditions  are  all 
necessary  and  sufficient,  they  are  not  all  independent.  We  obtain 
too  many  conditions.  In  many  cases,  however,  we  can  at  once 
reduce  them  to  the  proper  number  of  independent  conditions,  and 
when  this  is  difficult  we  can  have  recourse  to  the  second  method, 
to  be  given  in  the  next  Chapter,  which  is  free  from  this  objection. 

In  order  to  apply  this  method  with  success,  it  is  necessary  to 
have  some  convenient  methods  of  calculating  the  coefficients 

P    P       P 
j.  0,  ±  j...  j.  m. 

4.  The  first  method  which  suggests  itself  is  one  similar  to  that 
usually  given  to  determine  the  coefficients  of  the  equation  whose 
roots  are  the  squares  of  the  differences  of  the  roots  of  any  given 
equation. 

If  8lt  $2...  be  the  sums  of  the  first,  second  powers,  &c.  of  the 
roots  of  the  equation  fx  =  0,  we  have  by  Newton's  theorem 


where  p0  has  been  put  equal  to  unity.     If  St,  28...  be  the  sums  of 
the  powers  of  the  equation  F(x)  =  0,  we  have  in  the  same  way 


II.] 

we  may  also  prove 
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and  the  general  relation  can  be  found  without  difficulty. 
In  this  way  we  find 


1.2.3.4 


1.2 


But  the  process  becomes  longer  and  longer  at  every  stage. 
We  shall  therefore  proceed  to  point  out  some  other  methods  of 
obtaining  the  coefficients. 

This  method  of  proceeding  has  indeed  been  stated  only  because 
it  proves  in  a  convenient  way  that  when  p0=  1,  all  the  coefficients 
-Pj,  P2,  P3...  of  the  derived  equation  are  integral  rational  func- 
tions of  the  coefficients^,  pr..  pn. 

5.  The  equation  f  (x)  =  0  being  given,  to  calculate  the  coefficients 
o/F(x)=0. 

Put  x  =  y  ±  z  and  equate  separately  to  zero  the  sums  of  the 
even  and  odd  powers  of  z,  we  have 


Rejecting  the  root  z  =  0,  let  us  eliminate  z.  Then  the  roots 
of  the  resulting  equation  in  ?/  are  the  arithmetic  means  of  the 
roots  ""'  N  " 
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If,  on  the  other  hand,  we  eliminate  y  we  have  an  equation  of 
an  even  degree  to  find  z.  This,  putting  4s2  =£  is  the  equation 
whose  roots  are  the  squares  of  the  differences  of  the  roots  of  the 
given  equation. 

It  may  be  thought  that  this  elimination  may  prove  tedious, 
but  it  will  be  presently  shown  that  only  the  first  and  last  terms  of 
the  result  are  really  wanted.  All  the  others  may  be  omitted  in 
the  process  of  elimination,  and  thus  the  labour  will  be  greatly 
lessened.  The  method  is  however  most  useful  when  the  given 
equation  has  several  of  its  terms  absent. 

6.  Example.  To  determine  the  condition  that  the  roots  of  the 
biquadratic 

x4  +  px3  +  qx2  +  rx  +  s  =  0 

should  indicate  a  stable  motion. 

Applying  the  rule  we  have 
F  (x)  =  xs  +  3px5  +  (3/  +  2^)  a*  +  (4>pq  +  p*)  x3 

<f  —  4s)  x*  +  (pq*  +p*r  -  4<ps)  x  +  pqr  —  r2  -  p*s. 


The  first  four  coefficients  contain  only  positive  terms,  and  need 
not  be  considered.  If  the  last  three  coefficients  be  called  P4,  Pfi, 
Pfi,  we  have 

P2P4  -  4P6  =  (pq  -  2r)2 


If  then  P6  is  positive,  all  the  other  coefficients  are  positive. 

The  necessary  and  sufficient  conditions  of  stability  are  there- 
fore that  p}  q,  r,  s  should  be  finite  and  positive,  and 

P6  —pqr  —  r2  —p*s 
positive  or  zero. 

7.  In  forming  the  derived  equation  F(x)  the  only  difficulty 
is  to  form  the  last  term  Pm.  For  when  this  is  known  the  other 
terms  can  be  at  once  derived  from  it  by  an  easy  process. 

Let  a,  Z>,  c.  .  .  be  the  roots  of  /(#)  =  0  with  their  signs  changed, 
and  let 


Let  A  stand  for  the  operation 

A-£-tl+#f... 

da     db     dc 
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Po 


Then  since  —  =  abc  ...  we  have  obviously 
P*  J 


In  the  same  way  we  have 


And  generally 


and  so  on  up  to 


Po 


Let  us  now  operate  with  A  on  any  expression 


which  has  the  same  number  of  factors  in  every  term.     Let  r  be 
the  number  of  factors,  then  <£  may  be  written 


d 


d 


Let  P0,  Pl  .  .  .  Pm  be  the  coefficients  of  the  derived  equation 
F(x),  and  let  P0  =  1.  Then  since  the  roots  of  F(x)  are  the  sums 
of  the  roots  of  f(x)  taken  two  and  two,  it  is  easy  to  see  that 


&c.  =  &c. 

Thus  when  Pm  is  known,  the  other  terms  may  be  calculated 
without  difficulty.  The  term  Pn  will  be  called  the  fundamental 
term  of  the  equation. 

Example.     Given  in  the  case  of  a  biquadratic 


to  calculate  P.. 
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Performing  the  operation 

d    ,  0       d    t  0       d  d 

Ps  ~J~  +  2^2  T    +  3ft  T-  +  4»0  T—  > 
8^4         2(^3        -Idp9         °dPl' 

on  Pa  we  find  after  division  by  2, 


which  is  the  result  already  given. 

8.     It  should  be  noticed  that  in  the  equation 
f(x)  =  Poxn  H-^oT1  +  .  .  .  +pn  =  0, 


if  we  regard  x  as  a  number,  pQpl  ...  pn  are  all  of  equal  dimensions. 
It  follows  from  the  theory  of  dimensions,  that  if  any  subject  of 
operation  be  the  sum  of  a  number  of  terms  of  the  form 


there  must  be  the  same  number  of  factors  in  every  term.  For 
example,  in  every  term  of  the  expression  for  Pm  we  have 
a  +  /3  +  7  &c.  the  same. 

On  the  other  hand,  we  may  regard  x  as  a  quantity  of  one 
dimension,  and  in  this  case  p0,pl  >»pn  have  their  dimensions  in- 
dicated by  their  suffixes.  We  must  therefore  have  /3+  27+88+... 
as  well  as  a  +  ft  +  7  +  ...  the  same  in  every  term. 

These  two  tests  of  the  correctness  of  our  processes  will  be 
found  convenient. 

9.  The  whole  derived  equation  being  known  when  the  funda- 
mental term  is  known,  it  is  required  to  find  the  fundamental  term. 

[First  Method.] 

If  we  write  —  x  for  x  in  any  equation,  we  have  a  second  equation 
whose  roots  are  equal  and  opposite  to  those  of  the  first  equation. 
If  we  eliminate  x  between  these  two,  we  shall  get  a  result  which 
must  be  zero  when  the  two  equations  have  a  common  root.  The 
eliminant  must  therefore  contain  as  a  factor  the  product  of  the 
sums  of  the  roots  of  the  given  equation  taken  two  and  two. 

It  will  afterwards  be  shown  that  the  last  term  Pm  of  the 
derived  equation  (whenp0  is  put  equal  to  unity)  always  contains 
the  term 


with  a  coefficient  which  is  positive  and  equal  to  unity. 
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Hence  we  have  this  rule,  to  find  Pm,  eliminate  x  between 


and  divide  the  result  by  the  coefficient 

It  is  obvious  that  the  result  may  be  written  down  as  a 
determinant.  On  trial,  however,  it  will  be  found  more  convenient 
to  make  the  elimination  by  the  method  of  eliminating  the  highest 
aad  lowest  terms  than  to  expand  the  determinant. 

10.      Given  the  fundamental  term  of  the  equation  derived  from 
f(x)=p^-l+p^+...  +  p^  =  Q  ............  (1), 

to  find  the  fundamental  term  of  the  equation  derived  from 

Pjr+P^+.-.+Pn-iV+Pn^O    ............  (2). 

[Second  Method.] 

Let  Qn  be  the  product  of  the  sums,  two  and  two,  of  the  roots 
of  the  equation  (2)  taken  with  their  signs  changed,  so  that  Qn  is 
the  same  as  the  fundamental  term  of  the  derived  equation  and 
differs  from  Pm  only  in  having  a  suffix  more  convenient  for  our 
present  purpose. 

Let  Qn  be  expanded  in  a  series  of  powers  of  pn  :  thus 

Qn  =  fa+faPn  +  fapn+- 

where  fa,  fa,  fa,  &c.,  are  all  functions  of  p0,  pv  &c.,  which  functions 
have  to  be  found.  Let  the  roots  of  f(x)  —  0  with  their  signs 
changed  be  a,  b,  c...,  then 


Let  us  introduce  a  new  root,  which,  when  its  sign  is  changed, 
we  shall  call  r.     Then 


• 
Jro 

This  value  of  Qn  must  be  the  same  as  that  given  by  the  series 
when  we  write 


respectively,  for  p0,  p^  p3,  &c.,  pn_v  pn.     Equating  the  coefficients 
of  the  terms  independent  of  r  we  have 


4>  =  Q    p~^ 

1    Po 
R.  A. 


18  THE   FUNDAMENTAL   TERM.  [CHAP. 

Equating  the  terms  containing  the  first  power  of  r  we  have 


Substituting  for  <£0  we  have,  by  a  known  theorem  in  the  Differential 
Calculus, 

d  d  d   \  <)   , 

r  +  Pi  T-  +  •••  +P*  2  ~  -  • 

dp1  ldp,  -fltfW      Po 

Equating  the  terms  containing  the  second  power  of  r  we  have 


d   .       d 
+  (P>-z 


and  so  on.     Thus  we  have 


If  we  examine  this  process,  we  see  that  when  Qfl_1  is  known, 
we  may  at  once  write  down  the  terms  independent  of  pn  and  the 
coefficient  of  pn.  The  process  to  find  the  coefficient  of  p*  is 
longer,  but  it  may  be  much  shortened  by  the  consideration  that 
when  p0  =  1  the  result  must  be  an  integral  function  of  the  coeffi- 
cients. We  may  therefore  omit  all  terms  as  soon  as  they  make 
their  appearance,  which  do  not  contain  the  factor  pzn_l,  for  we 
know  that  such  terms  must  disappear  from  the  result. 

This  method  is  not  so  convenient  as  the  one  which  will  be 
presently  given  to  find  the  coefficients  of  the  higher  powers  of  pn. 
But  it  is  useful  as  showing  that  Qn  contains  the  term 


with  a  positive  integral  coefficient  equal  to  unity.     This  will  be 
clear  from  the  consideration  that  the  term  independent  of  pn  in 

Qn  is  obtained  from  Qn^  by  multiplying  by  ^-^  .    If  therefore  Qn_l 

Po 

contains  PlPz"J?-^  ,  Qn  must  contain  the  term  &&•;:#•-»  .     No 

other  term  can  be  formed  which  is  equal  to  this  with  an  opposite 
sign,  for  the  terms  which  enter  by  the  other  processes  to  be  per- 
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7) 

formed  on  Qn-l  all  contain  pn  as  a  factor.     Now  Q2=  —  '  ;  therefore 

Po 

Qs  contains  the  term  J-22,  Q  contains  CiBjpl  and  so  on. 

Po  Po 

It  has  been  shown  that  every  term  in  p0n~lQn  has  the  same 
number  of  factors  (Art.  8).  It  follows  from  this  reasoning  that 
this  number  of  factors  is  n  —  1. 

11.  To  find  the  fundamental  term  of  the  derived  equation  by 
means  of  a  differential  equation. 

[Third  Method.] 
The  fundamental  term  required  is  a  factor  of  the  eliminant  of 

1"2  +•••  = 


Let 


,  then  we  have 


n> 


=  o 


*  + 


-o 


If  we  write  y  +  dy  for  y  the  result  of  the  elimination  must  be 
same.     Hence  if  we  make 


even, 


-  1 


-  3 


dp3  _n  —  l  dp5     n  —  3 

id  if  E  be  the  eliminant,  we  have 


n  odd, 


It  follows  that  whether  n  be  even  or  odd,  E  must  satisfy  the 
equation 
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dE  dE 


=  0. 


We  may  make  the  elimination  by  multiplying  the  two  equa- 
tions by  y,  yz-..,  until  we  have  as  many  equations  as  we  have 
powers  to  eliminate.  If  in  the  determinant  thus  formed,  we 
multiply  out  the  terms  in  the  diagonal  joining  the  right-hand  top 

—  -i     —  +1 
corner  to  the  left  bottom  corner,  we  get  when  n  is  even  n  2     79  * 

n-l       n-l 

and  when  n  is  odd  pH  2  pQ  2  .     Now  Qn  must  contain  n  —  1  factors 

%  _  ]^ 
and  be  of  the  n  —^  —  th  degree.     Hence  when  n  is  even  E—  c^Q,, 

arid  when  n  is  odd  E=  cQn,  where  c  is  some  constant. 

Now  -j—  does  not  occur  in  the  above   differential   equation. 

Hence  treating  pl  as  a  constant,  we  see  that  Qn  must  satisfy  the 
differential  equation 


12.  We  may  show  that  p0n~lQn  is  a  symmetrical  function  of  the 
coefficients  pot  p^..,  pn  and  the  same  coefficients  read  backwards. 
Let  a,  6,  c...  be  the  roots  of  f(x)  =0  with  their  signs  changed, 
then  Qn  =  (a  +  b)  (a+c).... 

If  now  we  read  the  coefficients  in  the  opposite  order,  the  roots 
of  the  equation  thus  formed  will,  when  their  signs  are  changed,  be 

-  ,  r....     If  Q'n  be  the  fundamental  term  of  the  equation  derived 
a    b 

from  this,  we  have 


Since  -n  =  abc...  we  see  that 
Po 
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We  may  therefore  infer  that^?0n~1Qn  also  satisfies  the  differential 
equation 


13.     We  may  use  either  of  these  differential  equations  to  find 
Qn  when  §n_!  is  given. 

Let  the  first  differential  equation  be  represented  by 
and  let 


where  AQ,  A^  ...  are  functions  of  p0,  plt  &c.     The  value  of  AQ  has 
been  proved  in  Art.  10  to  be 


To  find  the  other  coefficients  of  the  powers  of  pn  substitute 
this  value  of  Qn  in  the  differential  equation  ;  we  have 


+  &c. 
Equating  the  several  powers  of  pn  to  zero,  we  find 


&C.  =  &C. 

Thus  by  one  regular  and   easy  process   each   term   may   be 
derived  from  the  other. 

In  performing  this  process  we  may  omit  every  term  in  the 
subject  of  operation  which  does  not  contain  £>n_2.     For  pn_z  can  be 

introduced  only  by  performing  -=—  ,  and  since  pn  is  absent  from  the 

™P» 
coefficients,  this  operation  yields  nothing. 
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Ill  this  way  we  find 


To  illustrate  this  process,  consider  how  Q5  is  obtained  from  Q4. 
The  first  line  is  formed  by  multiplying  the  line  above  by  p4,  this 
is  AQ.  To  find  the  coefficient  of  —p5  we  operate  with 

d  d.      /       d  d 


on  such  of  the  terms  in  the  line  above  as  contain  p3  and  then 
divide  by  p3.     Performing  the  same  operation  on  the  coefficient  of 

(—  jt?6)  we  obviously  obtain  the  coefficient  of  —     . 


In  M.  Serret's  Cours  d'AlgebreSupfrieure^otelll.,  there  will 
be  found  a  method  of  forming  the  last  term  of  the  equation  to  the 
squares  of  the  differences,  which  suggested  the  method  used  in 
Art.  13,  of  substituting  in  a  differential  equation,  if  only  a  differ- 
ential equation  could  be  found.  [See  also  Dr  Salmon's  Higher 
Algebra,  Arts.  60,  64,  and  72.] 


CHAPTER   III. 


Statement  of  the  theory  by  which  the  conditions  of  stability  of  a  dynamical 
system  with  n  co-ordinates  are  made  to  depend  on  2n  conditions. 
Arts.  1—4. 

A  rule  by  which  these  conditions  may  be  derived  one  from  another,  to- 
gether with  certain  other  true  but  not  independent  conditions.  Arts. 

5—8. 

A  rule  by  which,  when  the  coefficients  of  the  dynamical  equation  are  letters, 
the  2n  conditions  of  stability  may  be  inferred,  one  from  another, 
without  writing  down  any  other  conditions.  Arts.  9,  10. 

A  method  by  which  certain  extraneous  factors  may  be  discovered  and 
omitted.  Arts.  11,  12. 

Consideration  of  the  reserved  case  in  which  the  dynamical  equation  has 
equal  and  opposite  roots.  Arts.  13  —  18. 

A  few  geometrical  illustrations  not  necessary  to  tJie  argument.  Arts. 
19—25. 

Application  to  a  Dynamical  Problem.     Art.  26. 

1.  It  has  been  shown  in  the  first  Chapter  that  the  stability  of 
a  dynamical  system  with  n  co-ordinates  oscillating  about  a  state  of 
steady  motion  depends  on  the  nature  of  the  roots  of  a  certain 
equation  of  the  2/ith  degree  which  we  may  call 


The  system  is  stable  if  the  real  roots  and  the  real  parts  of  the 
imaginary  roots  are  all  negative.  Now  Cauchy  has  given  the  fol- 
lowing theorem  of  which  we  shall  make  some  use. 

Let  2=#+yV—  1  be  any  root,  and  let  us  regard  x  and  y  as 
co-ordinates  of  a  point  referred  to  rectangular  axes.  Substitute 
for  z  and  let 


Let  any  point  whose  co-ordinates  are  such  that  P  and  Q  both 
vanish  be  called  a  radical  point.     Describe  any  contour,  and  let 
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a  point  move  round  this  contour  in  the  positive  direction  and 

p 
notice  how  often  -~  passes  through  the  value  zero  and  changes  its 

sign.  Suppose  it  changes  a  times  from  4  to  —  and  ft  times  from 
—  to  4.  Then  Cauchy  asserts  that  the  number  of  radical  points 
within  the  contour  is  ^  (a  —  ft).  It  is  however  necessary  that  no 
radical  point  should  lie  on  the  contour. 

2.  Let  us  choose  as  our  contour  the  infinite  semicircle  which 
bounds  space  on  the  positive  side  of  the  axis  of  y.  Let  us  first 
travel  from  ?/  =  —  oo  to  y  =  +  <*>  along  the  circumference. 


we  have  changing  to  polar  co-ordinates 

/  (z)  =  P<f  (cos  n&  +  sm  n&  V^l)  4  . 
Hence 

P=p0r"  cos  n6  4  jp/1'1  cos  (n  -  1)  6  4 
Q  =  porn  sin  nO  4  p/1'1  sin  (n  -  1)  6  4 

In  the  limit,  since  r  is  infinite, 

p 

=  cot  nO  ; 


^  vanishes  when  nd  —  (%K  4  1)  —  ,  i.e. 


STT  5  TT 


?r  is  infinite  when  nd  =  %K  -  ,  i.  e. 


6<7r 


The  values  of  ^  in  series  (B)  it  will  be  noticed  separate  those 
in  series  (A). 

p 
When  6  is  small  and  very  little  greater  than  zero,  n-  is  posi- 

^ 

live,  and  therefore  changes  sign  from  4  to  —  at  every  one  of  the 
values  of  6  in  series  (A).  If  n  be  even  there  will  be  n  changes 
of  sign.  If  n  be  odd  there  will  be  n  —  1  changes  excluding 

6  =  ±  g  >  in  this  case  -~-  is  positive  when  6  is  a  little  less  than  —  , 
and  negative  when  9  is  a  little  greater  than  -y  . 
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Let  us  now  travel  along  the  axis  of  y  still  in  the  positive  di- 
rection, viz.  from  y  =  -f  oo  to  ?/  =  —  oo  .  Since  x  =  0  it  will  be 
more  convenient  to  use  Cartesian  co-ordinates,  we  have,  since 


= 
and 


and  -   P-P«-P>.-*?f+P»-j'-~> 

ami  .  •  •  ~>~i  ---  7  2~r      \~  • 

Q     ytpn-i-pn-^y  +  •••) 

The  condition  that  there  should  be  no  radical  point  within  the 
contour  is  that  this  expression  should  change  sign  through  zero 
from  —  to  +  as  often  as  it  before  changed  sign  from  +  to  —  on 
travelling  round  the  semicircle.  If  n  be  even  the  numerator  has 
one  more  term  than  the  denominator,  and  when  p0  and  p^  have 

the  same  sign,   -  ^  begins  when  y  is  very  great  by  being  negative. 

V 

In  order  that  it  should  change  sign  through  zero  n  times,  it  is 
necessary  and  sufficient  that  both  the  equations 


Pn  -Pn-^      + 

y5  -  •  .  .  =  0, 


should  have  their  roots  real,  and  that  the  roots  of  the  latter  should 
separate  the  roots  of  the  former. 

If  n  be  odd,  the  numerator  and  denominator  have  the  same 

p 
number  of  terms,  and  when  p0  and  p^  have  the  same  sign,  -^  begins 

V 

when  y  is  very  great  by  being  positive.  In  order  that  it  should 
change  sign  through  zero  from  —  to  +  n  —  1  times,  it  is  necessary 
and  sufficient  that  the  same  two  equations  as  before  should 
have  their  roots  real,  and  that  the  roots  of  the  former  should 
separate  the  roots  of  the  latter. 

In  order  then  to  express  the  necessary  and  sufficient  conditions, 
that  f(z)  =0  may  have  no  radical  point  on  the  positive  side  of 
the  axis  of  y,  put  z  =  y  V—  1  and  equate  to  zero  separately  the 
real  and  imaginary  parts.  Of  the  tiuo  equations  thus  formed,  the 
roots  of  the  one  of  lower  dimensions  must  separate  the  roots,  of  the 
other.  It  is  also  necessary  that  the  coefficients  of  the  two  highest 
powers  of  z  in  f  (z)  should  have  the  same  sign. 

3.  It  has  been  stated  that  p0  and  pl  the  coefficients  of  the  two 
highest  powers  in  f(z)  must  have  the  same  sign.  It  is  easy  to  see 
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p 

that,  if  they  had  opposite  signs,  -^  would  change  sign  through  zero 

^5 

2n  times  as  we  travel  round  the  contour.  All  the  radical  points 
of  the  equation  would  then  lie  on  the  positive  instead  of  the 
negative  side  of  the  axis  of  y. 

It  has  also  been  assumed  that  no  radical  point  lies  on  the 
contour.  It  has  therefore  been  assumed  that  /  (z)  =  0  has  no  root 
of  the  form  z  =  y  N  —  \.  It  will  be  more  convenient  to  consider 
this  exception  a  little  further  on. 

4.  It  is  required  to  express  in  an  analytical  form  the  conditions 
that  the  roots  of  an  equation  f2  (x)  =  0  may  be  all  real,  and  may 
separate  the  roots  of  another  equation  fj  (x)  =•  0  of  one  degree  higher 
dimensions. 

To  effect  this,  let  us  use  Sturm's  theorem  reversed.  Perform 
the  process  of  finding  the  greatest  common  measure  of  f\  (x)  and 
f^  (x),  changing  the  sign  of  each  remainder  as  it  is  obtained.  Let 
the  series  of  modified  remainders  thus  obtained  be  /3(#),  /4  (#), 
&c.  Then  it  may  be  shown  that  when  any  one  of  these  functions 
vanishes,  the  two  on  each  side  have  opposite  signs.  It  is  also  clear 
that  no  two  successive  functions  can  vanish  unless  f^  (x)  and  f2  (x) 
have  a  common  factor.  This  exception  will  be  considered  pre- 
sently. 

Hence  in  passing  from  #  =  —  oo  to  +  oo  no  variation  of  sign 
can  be  lost  except  when  /x  (x)  vanishes.  If  a  variation  is  lost  it  is 
regained  when  x  has  the  next  greatest  value  which  makes  yj  (x} 
vanish  unless  f2  (x)  =  0  has  a  root  between  these  two  successive 
roots  of  yj  (x}  =  0.  Hence  this  rule : — 

Tlie  roots  of  the  equations  ^  (x)  =  0,  f2  (x)  =  0,  will  be  all  real 
and  the  roots  of  the  latter  will  separate  those  of  the  former,  if  in 
the  series 

fjx),  f2(x),  fs(X)... 

as  many  variations  of  sign  are  lost  in  passing  from  x  =  —  oo  to 
x  =  +  oo  as  there  are  units  in  the  degree  of  the  equation  fj  (x)  =  0. 

We  have  supposed  the  variations  of  sign  to  be  lost  instead  of 
gained  in  passing  from  x  =  —  oc  to  +  oo .  That  this  may  be  the 
case  the  signs  of  the  highest  powers  of  ft  (x)  and  fa  (x)  must  be 
the  same. 

These  functions  are  alternately  of  an  even  and  odd  degree,  the 
condition  that  the  whole  number  of  variations  of  sign  may  be  lost 
in  passing  from  x  =  —  oo  to  -f  oo  may  be  more  conveniently  ex- 
pressed thus : — The  coefficients  of  the  highest  powers  of  x  in  the 
series 

f;(x),  f,(x),  f,(x)  ... 

must  all  have  the  same  sign. 
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5.  The  process  of  finding  the  greatest  common  measure  of 
two  algebraic  expressions  is  usually  rather  long.  We  may  in  our 
case  shorten  it  materially  by  omitting  the  quotients  and  perform- 
ing the  division  in  the  following  manner.  Let 


then,  since  p^  is  positive,  it  easily  follows  by  division  that 

/3  (x}  =  Axn~*  -  A'x"-*  +  A"xn-« 
where  ^  =p1p2-poPz, 

A'=Pd\-PQP»> 
&c.  =  &c., 

so  that  by  remembering  this  simple  cross-multiplication  we  may 
write  down  the  value  of  f3  (x)  without  any  other  process  than  what 
maybe  performed  by  simple  inspection.  In  the  same  way^(#), 
&c.  may  all  be  written  down. 

6.     Ex.  1.     Express  the  conditions  that  the  real   roots   and 
real  parts  of  the  imaginary  roots  of  the  cubic 

i 
may  be  all  negative. 


=      - 


/a 


The  necessary  conditions  are  that 

p,  pq—r  and  r 
must  all  be  positive. 

Ex.  2.      Express   the    corresponding    conditions    for   the   bi 
quadratic 

x*  +px*  +  qx*  +  rx  +  s  =  0, 
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The  conditions  are  that 

p,  pq  —  r,  (pq  —  r)  r—p*s  and  s 
must  be  all  positive. 

These  are  evidently  equivalent  to  the  five  conditions  that 

P>  q,  ^  s,   (pq-r)r-p*s, 
should  be  all  positive. 

In  both  these  examples  all  the  numerical  work  has  been 
exhibited. 

7.  Since  the  coefficients  of  the  highest  powers  of  x  in/x  (x)  arid 
fz  (x)  are  pQ  and  p  we  see  that  the  condition  that  p0  and  p±  should 
have  the   same  sign  is  included  in  the  general  statement  that  all 
the  coefficients  of  the  highest  powers  should  have  the  same  sign. 
If  the  function /(a?)  be  of  n  dimensions  we  thus  obtain  n  necessary 
and  sufficient  conditions. 

On  examining  these  conditions  in  the  cases  of  the  cubic  and 
biquadratic  it  will  be  seen  that  they  cannot  be  satisfied  if  any  one 
of  the  coefficients  of  the  given  equation  should  be  negative. 

8.  Although  the  theorem  in  its  present  form  gives  n  conditions 
as   the  proper  number  for  an  equation  of  the  nih  degree,  yet  it  is 
important  to  notice  that  it  gives  other  conditions  also  which  are 
true  and  may  be  useful.    It  has  been  shown  in  the  second  Chapter 
that  all  the  coefficients  of  the  equation  f(x)  =  0  must  be  positive, 
hence  the  roots  off^x)  =  0  must  all  be  positive.     It  may  be  shown 
also,  that  the  roots  of  each  of  the  functions  f^(x)}  fz(x),  &c.  are 
separated   by  the  roots  of  the  function  next  below  it  in  order. 
Hence  the  roots  of  all  these  functions  must  be  positive,  and  there- 
fore in  every  one  of  the  functions  the  coefficients  of  all  the  powers 
must  be   alternately    positive    and    negative   and    not    one    can 
vanish.     If  however  f^x)  and  fz(x)  have  one  or  more  common 
factors  some  of  the  functions /3(#),/4(V),  &c.  will  wholly  vanish. 

9.  When   the   degree  of  the    equation  is  very   considerable 
there  is  some  labour  in  the  application  of  the  rule  given  in  Art.  5. 
The  objection  is  that  we  only  want  the  terms  in  the  first  column, 
and  to  obtain  these  we  have  to  write  down  all  the  other  columns. 
We  shall  now  investigate  a  method  of  obtaining  each  term  in  the 

first  column  from  the  one  above  it  without  the  necessity  of  writing 
down  any  expression  except  the  one  required. 

We  notice  that  each  function  is  obtained  from  the  one  above  it 
by  the  same  process.  Now 
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The  first  and  second  lines  will  be  changed  into  the  second  and 
third  lines  by  writing  for 

Po»  ?'i'  P*>  Pv  &C* 

the  values 


If  then  in  any  term  of  any  function  we  make  these  changes,  we 
obtain  the  corresponding  term  of  the  function  next  in  order. 

10.     Example.     Express  the   conditions   of  stability   for   the 
quintic 


We  have 


~  Pi   PlP*  ~  P»P'          PlPi  ~ 


~  (PiP*~ PopJ 

11.  On  examining  the  conditions  as  given  in  the  cases  of  a 
biquadratic  and  quintic  it  will  be  apparent  that  several  contain 
the  previous  conditions  as  factors.  Thus  the  analytical  expressions 
are  rendered  much  longer  than  is  necessary.  It  is  now  proposed 
to  investigate  a  method  of  discovering  and  omitting  these  extraneous 
factors  as  they  occur,  and  thus  obtaining  the  required  conditions  in 
their  simplest  forms. 

Let  the  coefficients  of  the  several  powers  of  x  in  the  functions 
be  when  taken  positively 

fi(%)=Po,         pv          p^          Ps-" 

/2(*0  =  Pl>  Ptf  P*  Pi  — 

fz(x}  =  A,       A',       A",       A'"... 

£  I    \        ~D  T)'  7?"  7?//7 

/;(«)=: If,         B ,        £  ,        ±>    ... 

&c.  =  &c. 
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Let  us  first  find  which  of  these  terms  contain  p1  as  a  factor. 
Putting  ^  =  0  and  using  the  rule  in  Art.  5,  the  series  become 

~PoP5>          ~PoPT>  ~~PoP9>         


o,          o,  o,  o,       

o,          o,  o,  o,       

&c. 

Hence  the  C's  and  D's  all  vanish  and  therefore  contain  j^  as 
a  factor.  By  the  rule  in  Art.  5,  the  E's  contain  p*t  the  F's  con- 
tain p*,  the  G'spf,  the  H'sp*,  and  so  on. 

But  since  each  line  is  formed  from  the  preceding  by  a  uniform 
rule,  it  follows  that  the  D's  and  E's  contain  A  as  a  factor,  the 
F's  contain  A*,  the  G's  contain  A3,  the  H's  contain  A5,  and  so  on. 

The  factor  A  in  the  D's  and  E's  has  its  origin  in  the  factor  pl 
which  occurs  in  the  (7s  and  D's  and  would  not  appear  if  that 
factor  had  been  omitted  when  the  C's  and  D's  were  formed.  The 
factor  A  in  the  D's  and  E's  in  the  same  way  gives  rise  to  the 
factor  B  in  the  E's  and  .Fs.  So  that  if  we  take  care  each  time 
we  perform  the  process  described  in  Art.  9  to  omit  the  common 
factor  PI  whenever  it  occurs,  all  these  subsequent  factors  will  never 
make  their  appearance. 

We  shall  now  show  that  if  these  factors  are  omitted,  the 
dimensions  of  the  nih  function  fn(x)  will  be  n  — - —  .  First  consider 

the  actual  dimensions  of  each  function  before  the  factors  are 
omitted.  If  we  examine  the  rule  by  which  each  function  is 
derived  from  the  preceding,  it  will  become  evident  that,  the 
dimensions  of  each  letter  being  indicated  by  its  suffix,  the 
dimensions  of  any  function  are  equal  to  the  sum  of  the  two  pre- 
ceding +  2. 

In  the  following  table  the  first  column  indicates  the  function. 
In  the  second  column  will  be  found  the  dimensions  of  the  leading 
coefficient  of  that  function  when  calculated  by  the  rule  in  Art.  5. 
In  the  third  column  will  be  found  the  dimensions  as  given  by  the 

formula  n .     In  the  remaining  columns  are  the  dimensions 

of  the  extraneous  factors  p^  A,  B,  &c.  introduced  into  each  term. 
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/!« 

0 

1 

3 

0 
1 
3 

/4(#) 

6 

6 

/.(*) 

11 

10 

1 

/,(*) 

19 

15 

1 

3 

/7(x> 

32 

21 

2 

3 

6 

/8(#) 

53 

28 

3 

6 

6    10 

/9(<r) 

87 

36 

5 

9 

12    10    15 

&c. 

&c. 

&C. 

&c. 

&c. 

&c.   &c.   &c.   &c. 

Each  term  in  the  second  column  is  the  sum  of  the  two  terms 
just  above  it  +  2.  The  nih  term  in  the  third  column  is  equal  to 
the  term  just  above  it  +  (n  —  1).  In  all  the  other  columns  each 
term  is  the  sum  of  the  two  terms  just  above  it.  The  last  term  in 
the  wth  row  is  equal  to  the  (n—  3)th  term  in  the  third  column. 
We  wish  to  show  that  any  term  in  the  second  column  is  equal  to 
the  sum  of  all  the  terms  in  the  same  row  to  the  right  of  that  term. 
It  is  not  difficult  to  show  from  the  data  just  given  that  if  this  be 
true  for  any  two  adjacent  rows,  it  is  true  for  all  the  others,  and 
hence  we  may  assume  it  to  be  always  true. 

It  is  clear  that  these  extraneous  factors  may  be  omitted  since 
by  the  conditions  already  expressed  they  are  all  positive.  When 
omitted  as  they  occur,  the  dimensions  of  the  nih  function  has  just 

71—1 

been  shown  to  be  n  — —  .     It  is  easy  to  see  that  the  conditions 

2i 

thus  reduced  must  contain  the  terms 

P»  PiP*  PiP*P*>  PiPiPsPv&C' 

Now  if  we  take  any  one  of  these  as 

PiPtPsPtPiP* (1), 

and  operate  by  the  rule  in  Art.  9,  we  have 


which  contains  the  term 

(2)- 


Thus  we  have  p1  introduced  as  often  as  there  is  a  factor  in  (1) 
with  an  odd  suffix.  But  it  should  be  introduced  only  once.  These 
extra  p^s  are  the  extraneous  factors  to  be  omitted.  Each  of  these, 
if  left,  would  appear  as  the  factor  p^p2  —p0p3  in  the  next  condition, 
and  be  still  more  complicated  in  the  next  after  that. 
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In  order  then  to  obtain  the,  several  conditions  in  their  simplest 
form  it  is  only  necessary  after  performing  the  operation  described  in 
Art.  5  or  Art.  9  to  divide  by  p*,  where  K  is  one  less  than  the  number 
of  factors  with  odd  suffixes  in  the  condition  operated  on. 

12.  Example.  Express  the  conditions  of  stability  for  the 
sextic 


f(x)  =p0x«  +  Plx5  +p2x4  +  pBx3  +  p,x2  +  p5x 
We  have 


/4  H   =   (PlPlP*  -  P»P*  -  PlP* 

/a  0*0  =  {PtPiPsP* 


/8  («)  =  [PlPlPlPiPs  -P*P8*P*P*  -  Pl*P?P*  + 

-  P'P"  + 


X 

y*7  (^)  =  coefficient  of  x  in  f6  x  x  by  p6. 

13.     In  the  preceding  theory  two  reservations  have  been  made. 

1.  In  applying  Cauchy's  theorem  it  has  been  assumed  that 
there  were  no  radical  points  on  the  axis  of?/. 

2.  It   has  been  assumed   that  P  and    Q  have   no   common 
factor,  so  that  none  of  the  functions  /,,  /a,  &c.  vanish  absolutely. 

If  any  radical  point  lie  on  the  axis  of  y,  it  is  clear  thaty*(z)  =•  0 
must  have  a  factor  of  the  form  (z2  -I-  «2)r.  Let  f(z)  =  (z2  +  a2)r  </>(z). 
In  this  case  when  we  put  z  —  y  V  —  1,  we  have 


P=(<?-f)'F 


Thus  P  and  Q  have  a  common  factor,  and  we  are  warned  of  the 
possible  existence  of  radical  points  on  the  contour  by  the  total 
vanishing  of  some  one  of  the  functions  flt  fz,  /3,  &c. 

The  two  reserved  cases  may  therefore  be  included  in  one.  If 
y(s)  =  0  be  the  equation  furnished  by  dynamical  considerations, 
we  form  the  functions  /1?  /2,  /8,  &c.  If  all  these  be  finite,  the 
question  of  the  stability  of  the  system  has  been  answered.  If  any 
one  vanish  absolutely,  fv  (x)  and  fz  (x)  have  a  common  measure, 
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and  we  must  add  some  further  considerations.  It  will  be  con- 
venient to  examine  separately  the  dynamical  effect  of  the  roots 
which  do  not  and  which  do  enter  through  the  greatest  common 
measure.  Let  us  begin  with  the  former. 

14.     Following  the  same  notation  as  before,  we  have 


Q  =  ±/2  (y)  =p^y 

If  then/  (z)  have  two  roots,  viz.  ±(h  +  k  V  -  1),  which  are  equal 
and  opposite,  then  fl  (y}  and  /2  (y)  must  have  two  common  roots, 

viz.  +  -  -  —  .     The  common  measure  therefore  of  f  (y)  and 

V  —  1 

/a  (y)  contains  all  the  roots  of  f(y  V  —  1)  which  are  equal  and 
opposite.  Conversely  the  greatest  common  measure  of  P  and  Q  is 
necessarily  an  even  function  of  y,  and  if  it  be  -equated  to  zero,  its 
roots  are  necessarily  equal  and  opposite.  These  roots  must  also 
satisfy  f(y\l~^T)  =  Q. 

Let  this  greatest  common  measure  be  ^  (?/2)  =  0,  and  let  y*r  be 
the  highest  power  which  enters  into  it.     Also  let 

/(*)  =*(-*•)*«. 

then  (/>  (z)  is  a  function  which,  as  has  just  been  shown,  has  not  got 
two  roots  equal  and  opposite,  and  to  this  function  we  may  apply 
Cauchy's  theorem  without  fear  of  failure.  Putting  z  =  y  V  —  1,  let 


P' 

Then  we  wish  to  express  the  condition  that  -^-,  should  change 

sign  from  —  to  +  through  zero  n  -  2r  times  if  n  be  even  and 
n  —  2r—l  times  if  n  be  odd.     But 


and  f(y  V=T)  = 


pi 

Thus  the  number  of  changes  of  sign  in  —  ,  is  exactly  the  same 

P 

as  that  of  -^  .     The  factor  ty  (y*)  will  run  through  all  the  functions 

/2  (y)>  fs  (y)y  &c-  obtained  from  /j  (y)  by  a  process  which  is  equiva- 
lent to  that  of  finding  the  greatest  common  measure  of  fv  (y)  and 

R.  A.  3 
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fz  (y).     The  changes  of  sign  of  this  factor  will  therefore  not  affect 
the  number  of  variations  of  sign  in  the  series  flt  /2,  /3,  &c. 

The  last  factor  which  is  not  zero  is  /n+1_2r  (y}  if  n  be  the  dimen- 
sions offt  (y). 

Hence  if  we  omit  the  considerations  of  the  vanishing  factors  and 
apply  the  same  rule  as  before  to  the  n  +  1  —  2r  remaining  factors,  we 
can  express  the  condition  that  the  proper  number  of  changes  of  sign 
from  —  to  +  have  been  lost  through  zero  in  the  function  <£  (z),i.e.  that 
the  roots  not  given  by  the  vanishing  of  fn+1_2r  are  all  of  the  character 
to  ensure  stability. 

15.  Let  us  next  consider  the  effect  on  stability  of  the  roots 
indicated  by  the  absolute  vanishing  of  one  of  the  subsidiary  func- 
tions. This  function  must  be  of  the  form 


where  n  is  even.     The  corresponding  factor  of  f(z)  is 

F  0)  =  <ln  +  2«-/  +  <ln-/  +  "" 

The  roots  of  this  equation  are  two  and  two  equal  with  opposite 
signs,  it  is  therefore  necessary  for  stability  that  no  root  should 
have  any  real  part.  To  express  this  condition,  draw  a  straight 
line  parallel  to  the  axis  of  y  at  an  indefinitely  short  distance  from 
it,  viz.  x  =  h.  Let  us  apply,  in  the  same  manner  as  before,  Cauchy's 
theorem  to  the  contour  formed  by  this  straight  line  and  the  infinite 
semicircle  on  its  positive  side.  Putting  z  —  h  +  z,  we  have 


=  qn  + 
The  two  functions  are  therefore,  omitting  the  positive  factor  h, 


j  (y)  and/2  (y)  are  what  P  and  Q  become  when  arranged 
in  descending  powers  of  y  and  the  coefficients  of  their  highest 
powers  made  to  have  the  same  sign.  Hence 


,      -     dy 

The  rule  described  in  Art.  4  will  now  become  the  same  as 
that  usually  called  Sturm's  theorem.  We  are  to  seek  the  greatest 
common  measure  ofyj  (y)  and  its  differential  coefficient,  and  make 
the  coefficients  of  the  highest  powers  of  y  in  these  two  and  in  the 
series  of  modified  remainders  all  positive. 
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That  we  should  have  been  led  to  Sturm's  theorem  in  this  case 
is  just  what  we  might  have  expected.  For  to  express  the  condi- 
tions that  the  roots  of 


are  all  of  the  form  ±  r  V  —  1  is  the  same  thing  as  to  express  the 
conditions  that  the  roots  of 

2  I  „    ^  _      n 

are  all  real. 

16.  There  is  however  another  mode  of  proceeding.  Suppose 
we  have  calculated  the  functions  /j,  /2,  &c.  for  the  general  equa- 
tion 


and  find  when  the  values  of  p0,  plt  &c.  are  substituted  that  some 
one  function  say  fr  of  the  series  absolutely  vanishes,  and  therefore 
also  all  the  functions  which  follow  it.  Then  operate  on  each  of 
^these  vanishing  functions  with 


repeating  the  operation  until  we  obtain  a  result  which  is  not  zero. 
If  we  now  replace  these  vanishing  functions  by  these  results  we 
may  apply  the  rule  of  Art.  4,  just  as  if  these  were  the  functions 
supplied  by  the  process  of  the  greatest  common  measure.  As  this 
process  is  not  so  convenient  as  that  already  given  it  is  unnecessary 
to  consider  it  in  detail*. 

17.     As  a  numerical  example,  let  us  examine  whether  the  roots 
of 


f(x)  =  x8+  2^7  +  4*6  +  4*5  +  6^4+  6x3  +  7*2  +  4a?  +  2  =0 

satisfy  the  conditions  of  stability.  In  order  to  show  the  working 
of  the  method  it  will  be  necessary  to  exhibit  all  the  numerical 
calculations.  We  have  by  Art.  5, 

*  The  function  fr  (x)  vanishes  because  the  equation  /  (2)  =  0  has  two  roots  equal 
and  opposite.  If  we  put  z=d+h,  where  h  is  as  small  as  we  please,  this  peculiarity 
will  disappear.  Thus  if  the  values  of  z  are  of  the  form  ±(a±/3  <J~^V)  the  corre- 
sponding values  of  z'  are  -  7t±  a±£  *J  -  1.  These  values  of  z  will  indicate  stability 
if  a  be  zero  and  instability  if  a  have  any  value  positive  or  negative.  If  h  be  as  small 
as  we  please  and  positive,  the  values  of  z'  will  indicate  stability  or  instability  under 
the  same  circumstances.  We  may  therefore  apply  the  rule  of  Art.  4  to  the  func- 
tion /(z'+fe)  instead  of  f(z),  provided  we  retain  only  the  lowest  powers  of  h  which 
occur.  Hence  all  the  functions  /,  (x),  /2  (a;)...  /r-i(«)  which  do  not  vanish  are 
unaltered.  To  find  what  function  will  replace  fr  (x)  we  must  increase  by  h  all  the 
roots  of  /(z)=0  when  their  signs  have  been  changed.  This  may  be  effected  by 
performing  on  /  (z)  the  operation  represented  by  A  in  Art.  7.  The  rule  in  the  text 
therefore  follows  from  the  one  given  in  that  article. 

3—2 
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f1(x)  =  x*~  4^  +  6^-7^+2, 
/2  (x)  =  2a7  -  4#5  +  Go;3  -  4#  , 
/30*)  =  4a6- 
/4  (a?)  =  4#5  - 


Here  we  fmd/6  (x)  to  be  absolutely  zero,  accordingly  by  Art.  15 
replace  it  by  the  differential 
Sturm's  rule.     We  have  therefore 


we  replace  it  by  the  differential  coefficient  of  /6  (a?),  this  being 
W 


/»~-8M20. 

We  see  that  the  two  last  of  the  coefficients  of  the  highest 
powers  are  negative.  The  roots  therefore  do  not  satisfy  the  con- 
ditions of  stability. 

As  another  example,  take  the  equation 

fx  =  XQ  +  x5  +  6#4  +  5  a?  +  1  1  x*  +  6x  i-  G  . 

Here 


/2  (x)  =  a*  - 

/3(*)  =  tf4-5tf2+6, 

/4(*)  =  0. 

Eeplacing  /4  (#)  by  the  differential  coefficient  of/3(#),  we  have 

/«00- 

/8  (SB)  = 


Here  all  the  coefficients  of  the  highest  powers  are  positive, 
hence  the  roots  satisfy  the  conditions  of  stability. 

It  is  clear  that  when  the  coefficients  are  numerical  the  rule 
given  in  Art.  5  is  the  most  convenient,  but  when  the  coefficients 
are  letters,  the  rule  in  Art.  9  will  be  found  preferable. 

The  process  would  be  simplified  by  omitting  the  alternate  posi- 
tive and  negative  signs  of  the  terms  in  each  line. 
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18.     It  may  be  interesting  to  express  the  two  subsidiary^  func- 
tions ft  (x)  and  fa  (x)  in  terms  of  the  roots  of  the  given  equation. 

Let  ttj,  a2,  a3. .  .an  be  the  roots  of  the  given  equation  f(x)  —  0,  so 
that 

/(a;)  =  OB-  a,)  (x-  «2)  (aj-oj... 


Then  it  is  evident  that 
±/t  (x  v^l)  =  x 

=  J  (x  +  aj  (x  +  aa)  . . .  +  i  (a;  -  aj  (a?  -  a2)  ... 
It  may  be  shown*  that 

±  *  ^~l        = 


19.  The  following  propositions  are  not  necessary  to  the  main 
argument,  but  as  they  illustrate  geometrically  the  propositions  in 
this  chapter  it  has  been  considered  proper  to  state  them  very 
briefly.     The  demonstrations  will  therefore  be  much  curtailed. 

The  equation  being 

we  put  as  in  Art.  1,  z  =  x  +  yJ—~L,  and  thus  obtain  two  curves, 
whose  equations  expressed  in  polar  co-ordinates  are 

P  =pQrn  cos  nO  +plrn'1  cos  (n  -  1)  6  +  . . .  =  0} 
Q  =  p0rnsmn0  +  plrn~lsm  (n  -  1)  0  +  ...  =  Oj  ' 
These  intersect  in  the  radical  points  of  the  equation/  (z)  =  0. 

20.  If  we  trace  these  curves  we  find  that  the  curve  P  =  0 
has  n  asymptotes  whose  directions  are  given  by  cos  nO  =  0,  i.  e. 

3?r        5  TT 


<9=- 

Let  us  assume 


A1(x-a.2)...(x-  an)  +Az(x-al)...(x-an)  +  &G., 

where  Alt  A^,  &c.  are  constants  whose  values  have  to  be  found.     Putting  x=aly 
we  have 


But  since  xn  +plxn~1  +  .  .  .  =  (x  -  a^  (x  -  a2)  .  .  .  (x  -  an)  , 

we  have  by  putting  x=at  and  x=  -  ^ 


Subtracting  the  second  of  these  results  from  the  first  we  find  Al  to  have  the  value 
given  in  the  text. 
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These  asymptotes  all  pass  through  the  same  point  on  the  axis 

of  x,  viz.  x  =  —  2t  .     It  is  also  clear  that  only  one  branch  of  the 
np0 

curve  can  go  to  each  end  of  an  asymptote.  Similar  remarks  apply 
to  the  curve  Q  =  0,  the  directions  of  its  asymptotes  being  given  by 
sin  n0  =  0,  i.  e. 


' 


21.  From  these   simple  propositions   we   might,  if  it    were 
worth  while,  deduce  that  every  equation  must  have  a  root.     The 
asymptotes  of  the  two  curves  P  =  0,  Q  =  0  are  alternate,  and  no 
two  branches  of  the  same  curve  can  approach  the  same  end  of  an 
asymptote.     By  sketching  a  figure,  it  may  be  easily  shown  that 
some  branch  of  the  P  curve  must  cut  some  branch  of  the  Q  curve. 

22.  Let  us  next  consider  the  intersections  of  the  curves  P  =  0, 


If  we  transform  the  origin  to  Ji,  k,  we  put  x  =  h  +  g,  y 
This  is  the  same  as  expanding  f(h  +  &V—  1+f  +  ^V—  1),  and 
collecting  into  two  parcels  the  real  and  imaginary  terms.  Let  the 
expansion  be 


where  AQA^...  are  of  the  form 

c  (cos  a  +  sina  V—  1). 

If  we  put  f  -f  77  V—  l  =  r  (cos  0  +  sin  6  V—  1),  we  have 
P  =  c0  cos  «0  +  c/  cos  (6  -4-  oj  +  c/2  cos  (20  +  a2)  +  .  .  .) 
Q  =  c0  sin  a0  +  cjr  sin  (0  +  aj  +  c2r2  sin  (2(9  +  aj  +  ...]' 

If  the  point  (h,  k)  be  a  point  of  intersection  c0  =  0.  If  the  inter- 
section be  a  double  point  on  either  curve,  the  terms  of  the  first 
degree  must  be  zero,  therefore  cx  =  0,  and  the  origin  is  therefore 
a  double  point  on  the  other  curve  also. 

It  is  not  difficult  to  show  that  if  the  intersection  be  a  multiple 
point  of  any  degree  of  multiplicity  on  one  curve,  it  is  a  point  of 
the  same  degree  of  multiplicity  on  the  other  curve.  The  tangents 
to  these  branches  all  make  equal  angles  with  each  other,  the 
tangents  to  the  P  and  Q  curves  being  alternate  as  we  travel  round 
the  point  of  intersection.  If  the  intersection  be  not  a  multiple 
point  on  either  curve,  ihe  branches  cut  at  right  angles. 

Let  us  travel  round  a  point  of  intersection  along  the  circum- 
ference of  a  small  circle  whose  centre  is  the  point  of  intersection 
in  the  direction  in  which  0  is  measured.  Then  it  may  be  shown 
that  as  we  pass  from  a  P  curve  to  a  Q  curve,  P  and  Q  have 
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opposite  signs,  and  as  we  pass  from  a  Q  curve  to  a  P  curve  the  same 
sign.     This  is  in  fact  merely  Cauchy's  rule  for  the  changes  of  sign 

p 
of  -ft  when  we  travel  round  a  radical  point. 

23.  Let  us  express  the  condition  that  there  is  no  radical  point 
on  the  positive  side  of  the  axis  of  y.      This  is  the  geometrical 
illustration  of  Art  2. 

Draw  a  circle  of  infinite  radius,  and  let  it  cut  the  asymptotes 
of  the  P  curve  in  Plf  P2,  P3,...P?n  and  the  asymptotes  of  the 
Q  curve  in  Qt,  Q2,...(?2n.  These  points  alternate  with  each  other. 
Taking  only  those  points  which  lie  on  the  positive  side  of  the 
axis  of  y,  the  P  and  Q  curves  may  be  said  to  begin  at  these  points 
and  are  to  intersect  each  other  only  on  the  negative  side  of  the 
axis  of  y.  The  branches  of  the  two  curves  must  therefore  remain 
alternate  with  each  other  throughout  the  space  on  the  positive 
side  of  the  axis  of  y.  Their  points  of  intersection  with  the  axis 
of  y  must  be  also  alternate,  and  hence  if  we  put  x  =  0,  in  the 
equations  P  =  0,  Q  =  0,  and  regard  them  as  equations  to  find  y 
the  roots  of  each  must  separate  the  roots  of  the  other. 

Conversely,  we  may  show  that  if  the  intersections  of  the  two 
branches  are  alternate  on  the  axis  of  y,  they  cannot  have  intersected 
on  that  side  of  the  axis  of  y  on  which  the  common  intersection  of 
all  the  asymptotes  is  not.  This  is  the  result  arrived  at  in  Art.  2. 

24.  The  following  diagrams  exhibit  the  forms  of  the  curve 
P=0,    Q  =  0  for  a  biquadratic.     The  dotted  lines  represent  the 
asymptotes. 


A  biquadratic  with  four  imaginary  roots. 
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A  biquadratic  with  equal  imaginary  roots. 

25.     [If  z  —  x  +  y  V  —  1,  we  have 


Differentiating  with  respect,  firstly  to  x  and  secondly  to  y,  we 

/-»         i  /%*   *     \  Cv^r  Ct\\r        I 

find  /W  =  ^;  +  £v- 


It  easily  follows  that  -7-  =  -~=  and  -^  —  —  ^—  ;  so  that  both 
ax      ay  ax          dy 

the  functions  P  and  Q  satisfy  the  equation 


_ 

df     dx*  ~ 

The  equation/^)  =  0  gave  us  two  curves  which  we  have  called 
P  =  0  and  Q  =  0.  In  the  same  way  the  derived  equation  /'(z)  =  0 
will  give  us  two  other  curves,  which  we  may  represent  by  P  =  0 
and  Q'  =  0.  These  we  may  call  the  derived  P  and  Q  curves. 

If  as  in  Art.  22  we  transform  the  origin  to  the  point  (h,  A?)  we 
have 


F  =  Cj  cos  ^  +  2c2r  cos(0+  «2)+  3c/2  cos(2<9  +  a3)  +  ... 
Q'  =  G!  sin  «!  +  2c2r  sin  (6  +  a2)  +  3c3r2  sin  (20  +  cra)  +  .  .  . 

If  the  origin  be  at  a  point  of  intersection  of  the  curves  P  =  0, 
P'  =  0  which  is  not  a  double  point  on  the  first  of  these  curves,  we 
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have  C0  cos  a0  =  0  and  cos  ^  =  0.  Hence  the  tangent  to  the  curve 
P  =  0  at  the  point  of  intersection  is  parallel  to  the  axis  of  x. 
Conversely  if  we  move  the  origin  to  any  point  on  the  curve  P  =  0 
at  which  the  tangent  is  parallel  to  the  axis  of  x,  we  find  that  the 
curve  P'  =  0  passes  through  the  origin.  Hence  the  derived  P 
curve  passes  through  all  those  points  on  the  curve  P  =  0  at  which 
the  tangent  is  parallel  to  the  axis  of  x,  and  all  those  points  on  the 
curve  Q  =  0  at  which  the  tangent  is  parallel  to  the  axis  of  y.  In 
the  same  way  the  derived  Q  curve  passes  through  those  points  on 
the  curve  Q  =  0  at  which  the  tangent  is  parallel  to  the  axis  of  x, 
and  those  points  on  the  curve  P=0  at  which  the  tangent  is 
parallel  to  the  axis  of  y. 

If  c0  cos  Of0  =  0,  and  ct  =  0,  the  origin  is  a  double  point  on  the 
curve  P  —  0  and  the  origin  also  lies  on  both  the  curves  P'  =  0, 
Q'  =  0.  So  that  both  the  derived  curves  cut  the  curves  P  =  0  and 
Q  =  0  in  their  double  points.  In  other  words,  these  double  points 
are  radical  points  for  the  derived  equation  /'(z)  =0.  In  the  same 
way,  any  multiple  point  on  either  curve  is  a  multiple  point  of  one 
degree  less  multiplicity  on  bath  the  derived  curves. 

If  a  finite  straight  line  AB  be  drawn  parallel  to  the  axis  of  x 
joining  two  points  A,  B  on  the  same  or  on  different  branches  of 
the  curve  P  =  0,  this  finite  straight  line  must  cut  one  or  more 
branches  of  the  derived  P  curve.  For  it  is  clear  that  if  P  vanishes 

at  A  and  B,  P  which  is  equal  to  -y-  cannot  keep  one  sign  between 

U/X 

A  and  B,  and  must  therefore  vanish  somewhere  between  A  and  B. 
If  A  and  B  be  adjacent  points,  i.e.  if  there  be  no  other  points 
between  A  and  B  belonging  to  the  curve  P  =  0,  then  the  straight 
line  AB  must  cut  an  odd  number  of  branches  of  the  derived  P 
curve.  In  the  same  way  if  a  straight  line  CD  be  drawn  parallel 
to  the  axis  of  y  joining  two  adjacent  points  on  a  derived  Q  curve, 
this  straight  line  must  also  cut  an  odd  number  of  branches  of  a 
derived  P  curve  between  C  and  D. 

By  considering  a  tangent  as  the  limit  of  secant,  it  again  follows 
that  if  a  tangent  be  drawn  to  the  curve  P=  0  parallel  to  the  axis 
of  x,  the  derived  P  curve  must  pass  through  the  point  of  contact. 

Let  R  be  a  radical  point  on  the  derived  curve  f  (z]  =  0,  and  let 
it  not  be  a  double  point  on  either  of  the  curves  P=  0,  Q  =  0.  Let 
a  straight  line  be  drawn  from  R  in  any  direction  cutting  the 
branches  of  either  of  the  curves  P=  0,  Q  =  0  in  the  points  Alt  Az, 
&c.  Then  we  may  show  that 


so  that  the  polar  line  of  R  with  regard  to  either  of  the  curves 
P  =  0,  Q  =  0  is  at  infinity. 


42  DYNAMICAL   ILLUSTRATION.  [CHAP. 

The  positions  of  the  radical  points  of  the  derived  equation 
f  (z)  —  0  relatively  to  any  branch  or  branches  of  the  curves  P  =  0, 
Q  =  0  may  be  found  by  the  use  of  Cauchy's  theorem.  If  the 
point  spoken  of  in  Art.  1  travel  along  a  branch  of  the  curve  P  =  0, 

Z)'          *J 

it  is  easy  to  see  that  TV  =  -r^  •     If  it  travel  along  a  branch  of  the 

P'         dx 
curve  Q  =  0  we  have  7y=~j~*     If  then  any  contour  be  partly 

bounded  by  branches  of  these  curves,  the  simplest  inspection  of 
the  points  at  which  the  tangents  are  parallel  to  the  axis  will 

P 

determine  the  changes  of  sign  of  -j  as  it  passes  through  zero. 

If  another  part  of  the  contour  be  an  arc  of  a  circle  of  infinite  radius 
whose  centre  is  the  origin,  the  changes  of  sign  through  zero  will 
be  from  +  to  —  and  their  number  will  be  indicated  by  the  number 
of  asymptotes  of  the  derived  P  curve  which  cut  the  arc.] 

26.  The  use  of  Watt's  Governor  in  the  steam  engine  is  too 
well  known  to  need  description.  It  has  however,  as  commonly 
used,  a  great  defect.  It  is  sometimes  of  importance  that  the 
engine  should  continue  to  work  at  the  same  rate  notwithstanding 
great  changes  in  the  resistances.  Suppose  the  load  suddenly 
diminished,  the  engine  works  quicker,  the  balls  diverging  cut  off 
the  steam,  and  the  engine,  after  a  time,  again  works  uniformly, 
but  at  a  different  rate  from  before.  The  balls  as  they  open  out 
or  close  in  are  usually  made  to  describe  circles.  Let  them  now  be 
constrained  to  describe  some  other  curve  which  we  may  afterwards 
choose  so  as  to  correct  the  above  defect.  If  this  curve  be  a  parabola 
and  the  balls  be  treated  as  particles,  it  is  clear  from  very  elemen- 
tary considerations  that  these  will  be  in  relative  equilibrium  only 
when  the  engine  works  at  a  given  rate.  This  principle  is  due 
to  Huyghens,  see  Astronomical  Notices,  December,  1875.  It  is 
now  proposed  to  determine  the  condition  of  stable  oscillation  about 
a  state  of  steady  motion. 

Two  equal  rods  AB,  AB'  are  attached  at  A  by  hinges  to  a 
small  ring  which  can  slide  smoothly  along  a  vertical  axis.  The 
ring  is  attached  by  a  rod  to  the  valve  and  can  thus  govern  the 
amount  of  steam  admitted.  Two  equal  balls  are  attached  at  B 
and  B',  and  the  centre  of  gravity  G-  of  the  rod  AB  and  the  ball  B 
is  constrained  to  describe  some  curve.  To  represent  the  inertia  of 
the  engine  we  shall  suppose  a  horizontal  fly-wheel  attached  to  the 
vertical  axis  whose  moment  of  inertia  about  the  axis  is  /.  Let  the 
excess  of  the  action  of  the  steam  over  the  resistance  of  the  load 
be  represented  by  some  couple  whose  moment  about  the  vertical 
axis  is  /  (0),  where  6  is  the  inclination  of  the  rod  A  G  to  the 
vertical,  and  /  is  a  function  which  depends  on  the  construction  of 
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the  engine.      Since  the  steam  is  cut  off  when  the  balls  open  out, 
it  is  clear  that  —.. ^  -  is  negative. 

There  may  be  also  some  resistances  which  vary  with  the 
velocity.  Let  these  be  represented  by  a  couple  B  -^  tending  to 
retard  the  motion  round  the  vertical  axis  and  a  couple  round  A  in 
the  plane  BAB  equal  to  m  G  -r- ,  where  <p  is  the  angle  the 
vertical  plane  BAB  makes  with  a  fixed  vertical  plane. 

Let  m  be  the  mass  of  either  sphere  and  rod;  k  the  radius  of 
gyration  about  an  axis  through  G  perpendicular  to  the  rod,  and  k' 
that  about  the  rod,  let  I  =  A  G.  Then  the  equation  of  angular 
momentum  gives 


, ^  "£  +  2m  (F  +  ?  sin26>  +  k'2  cos2  0)  ^   =/(0)  -  B^ . 
at  [at  Cut  J  at 

Let  the  steady  motion  be  given  by  0  —  a,  -?  =  n,  and  let  the 
oscillations  be  represented  by  0  =  a  +  oc,  -j^  =  n  +  y.    We  have  then 


The  equation  then  reduces  to 


{/  +  2m  (F+P  sin2  OL  +  k'2  cos2  a))       +  2m  (P  +  tf  -  k'2)  sin  2*n 

=/(a)x- 
This  equation  may  be  briefly  written 


Let  z  be  the  altitude  of  G  above  some  fixed  horizontal  plane. 
Then  if  T  be  the  semi  vis  viva 

2  T=  1$'*  +  2m  {(k2  +  I2)  sin'0  +  k'2  cos20J  </>'2  + 


If  U  be  the  force  function,  omitting  the  couple  of  resistance, 
we  have  17=  —  2mgz.  The  virtual  moment  of  the  couples  of 
resistance  being  %mCff$0,  the  Lagrangian  equation  of  motion 
becomes 


ddT     dT__e>    rff 
dbdff     dd~  dO 
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Substituting  for  T  and  U  we  have  when  the  motion  is  steady 


Since  I  sin  6  is  the  distance  of  G  from  the  vertical  axis,  we  see 
that  the  path  of  G  must  be  a  parabola.     The  semi  latus  rectum  is 


which,  we  notice,  is  independent  of  the  radius  of  the  balls.  The 
length  of  this  latus  rectum  must  of  course  be  adjusted  to  suit  the 
particular  rate  at  which  the  engine  is  intended  to  work. 

When  the  system  is  oscillating  about  the  state  of  steady  motion 
we  have,  putting 


and  rejecting  the  squares  of  as  and  y, 

(k*  +  a2  cos2  a)  ^  +  C^  -  (P  +  k*  -  &'2)  sin  2am/  =  0. 

The  term   #,   it   will   be   noticed,  has  disappeared  from  the 
equation.     This  equation  may  be  briefly  written  in  the  form 


Eliminating  y  from  the  equations  of  motion  we  have 
~*  +  (AC+£H)~ 

The  coefficients  are  all  positive,  the  necessary  and  sufficient 
condition  of  stability  is  therefore 

(AC+BH)(BG+ML)>AFffL. 

In  some  clocks  to  which  Watt's  Governor  is  applied,  there  is  a 
special  arrangement  which  causes  C  to  be  much  greater  than  B. 
See  the  Astronomical  Notes,  xi,,  1851,  and  the  Memoirs  of  the 
Astronomical  Society,  Vol.  xx.  Neglecting  therefore  B,  we  have 


.'.  2  Cmn  (P  +  #  -  F)  sin  2*  >  F(k*  +  a2  cos2  a). 


CHAPTER   IV. 


Formation  of  the  equations  of  steady  motion  and  of  small  oscillation 
where  Lagrange's  method  may  be  used.  Arts.  1  —  5. 

The  equations  being  all  linear  the  conditions  of  stability  are  expressed 
by  the  character  of  the  roots  of  a  determinantal  equation  of  an  even  order. 
Art.  6. 

Mode  of  expanding  the  determinant.     Art.  7. 

A  method  of  finding  the  proper  co-ordinates  to  make  the  coefficients  of 
tJie  Lagrangian  function  constant.  Arts.  8  —  10. 

How  the  Harmonic  oscillations  about  steady  motion  differ  from  those 
about  a  position  of  equilibrium.  The  forces  which  cause  the  difference  are 
of  the  nature  of  centrifugal  forces  produced  by  an  imaginary  rotation 
about  a  fixed  straight  line.  Arts.  11  —  19. 

Reduction  of  the  fundamental  determinant  to  one  of  fewer  rows  by  t/ie 
elimination  of  all  co-ordinates  which  do  not  appear  except  as  differential 
coefficients  in  tJie  Lagrangian  function  ;  with  an  example.  Arts.  20  —  23. 

Formation  of  the  equations  of  Motion  and  of  the  determinant  when 
the  geometrical  equations  contain  differential  coefficients,  so  that  Lagrange's 
met/iod  cannot  be  used  ;  with  an  example.  Arts.  24  —  27. 

1.  Let  the  system  be  referred  to  any  co-ordinates  f  ,  rj,  f,  &c. 
The  general  expression  for  the  kinetic  energy  is 


where  P,  Q,  &c.  are  known  functions  of  f,  r\,  f,  &c.  and  accents 
have  their  usual  meaning.  Let  us  suppose  the  system  to  have 
some  motion  represented  by 

&C. 


and  when  disturbed,  we  wish  to  find  the  oscillations  about  this 
motion.     To  effect  this,  we  put 
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where  0,  <j>,  &c.  are  all  small  quantities.     Substituting  and  expand- 
ing T  in  powers  of  6,  </>,  &c.,  we  find 


+  &c. 

In  the  same  way  we  may  make  an  expansion  for  the  Potential 
Energy  of  the  forces,  viz. 


when  these  two  functions  are  given  the  whole  dynamical  system 
and  the  forces  are  known;  and  we  may  form  the  equations  of 
motion  by  Lagrange's  method. 

2.  We  shall  here  however  limit  the  question  by  supposing 
that  the  motion  about  which  the  system  is  oscillating  is  what  has 
been  called  in  Chap.  i.  steady.     The  analytical  peculiarity  of  such 
a  motion  is  that  when  referred  to  proper  co-ordinates,  every  coeffi- 
cient in  each  of  these  two  series  is  constant,  i.e.  independent  of  t. 
As  already  explained  the  physical  peculiarities  are  that  the  vis 
viva  is  constant  throughout  the  steady  motion  and  the  same  oscil- 
lations follow  from  the  same  disturbance  at  whatever  instant  it 
may  be  applied  to  the  motion.      A  method  of  discovering  the 
proper  co-ordinates,  if  unknown,  will  be  given  a  little  further  on. 

3.  In  order  to  form  the  equations  of  motion  we  must  now 
substitute  in  Lagrange's  equations 

dLdT_dT     dV_ 
dt  d&      dd  +  d6  ~ 

&c.  =  0, 

rejecting  all  the  squares  of  small  quantities.  The  steady  motion 
being  given  by  6,  0,  &c.  all  zero,  each  of  these  must  be  satisfied 
when  we  omit  the  terms  containing  0,  <£,  &c.  We  thus  obtain  the 
equations  of  steady  motion,  viz. 


&c.  =  &c. 
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These  equations  may  be  simply  formed  in  any  case  by  the 
following  rule.  Putting  L=T—V,  so  that  L  is  the  difference 
between  the  kinetic  and  potential  energies,  expand  L  in  powers  of 
the  co-ordinates  0,  (j>,  &c.  regarding  0',  <£',  &c.  as  zero.  The  required 
relations  are  obtained  by  equating  the  coefficients  of  the  first 
powers  to  zero.  This  rule  may  be  also  expressed  thus.  Let  L  be 
the  general  expression  for  the  excess  of  the  kinetic  energy  over  the 

ftential  energy  of  a  dynamical  system  in  terms  of  its  n  co-ordinates 
7j,  &c.     Let  this  system  be  moving  in  steady  motion  with  constant 

A  fc      A 

values  of  -=|  ,  -^  ,  &c.*     Then  substituting  these  constant  values  in 

the  general  expression  for  L,  the  relations  between  the  constants  of 
steady  motion  are  given  by 

dL  dL 


In  this  way  we  obtain  in  general  as  many  equations  as  there 
are  co-ordinates.  Usually  the  coefficients  in  the  expression  for  T 
are  constant  because  some  of  the  co-ordinates  are  constant  in  the 
state  of  steady  motion,  and  the  other  co-ordinates  appear  in  the 
expressions  for  T  and  V  only  as  differential  coefficients.  In  such 
cases  we  have  clearly  fewer  equations  than  co-ordinates  to  con- 
nect the  constants  of  steady  motion.  We  have  then  a  system  of 
possible  steady  motions  which  we  may  conveniently  term  parallel 
steady  motions. 

4.  To  obtain  the  equations  to  the  oscillatory  motion,  we  retain 
the  first  powers  of  6,  </>,  &c.  We  thus  obtain  a  series  of  equations 
of  which  the  following  is  a  specimen  : 


4-  ............  =0. 

To  solve  these  we  write 

&c. 


Substituting  we  obtain  on  eliminating  the  ratios  Ml  :  Mz  :  M3  &c. 
a  determinantal  equation,  viz. 

*  Since  we  may  in  Art.  1  change  the  co-ordinates  from  £,  17,  £,  &c.  to  £„  77^  flt 
<fec.,  where  £&&/(£,)  77  =  ^(77^  &c.,  it  is  clear  that  the  steady  motion  can  be  always 
expressed  by  constant  values  of  the  differential  coefficients  of  the  co-ordinates. 
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&c. 


&c. 


&c. 


&c. 


&c. 


&c. 


&c. 


=  0. 


This  equation  will  be  referred  to  &s  the  Determinantal  equa- 
tion. 

5.  If  we  refer  to  the  equation  formed  by  this  determinant 
and  read  it  in  horizontal  lines,  we  have  of  course  the  several 
equations  of  motion,  each  term  being  the  coefficient  of  0,  <£,  -xjr,  &c. 
in  order.  In  this  form  the  equations  may  be  reproduced  by  the 
following  easy  rule. 

Taking  the  expression  for  T—  V  as  given  in  Art.  1,  let  us 
consider  only  the  terms  of  the  second  order,  those  of  the  first 
having  been  already  used  to  determine  the  steady  motion  as  ex- 
plained in  Art.  3.  Separate  from  the  rest,  the  even  powers  of 
0,  0',  (/>,  <£',  &c.  and  write  for  &\  0'$',  &c.  -  Z>2<92,  -  D*0(f>,  &c.,  so 

that  D  will  stand  either  for  -^  or  for  the  m  in  the  determinant 

at 

of  Art.  4,  when  we  write  0  =  MjS"*t  <j>  =  M2emt,  &c.     Let  the  sum 
of  these  terms  be  called  P,  so  that 

Let  the  remaining  portion  of  the  terms  of  T  —  F,  viz.  those  con- 
taining both  0}  $,  &c.  and  0',  <£',  &c.,  be  called  Q,  so  that 


Then  the  several  equations  may  be  formed  from  the  rule 

dP     dQ_      d 

d0  "*"  d0          d 


0 


dP  +  dQ_J)d_Q^  =  0 

&c.  =  0  J 

In  applying  this  rule  no  accented  letters  will  occur  except  in 
the  second  term  of  each  equation.     If  we  wish  D  to  stand  for  m 
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in  the  determinant,  we  must  regard  ff,  <£',  &c.  as  abbreviations  for 
DO,  Dcf),  &c.  If  we  wish  -to  use  the  equations  themselves,  we 
replace  DO,  D<j>,  &c.  by  fft  </>',  &c. 

[The  determinant  may  also  be  found  by  another  rule.  Taking 
as  before  only  the  terms  of  the  second  order  in  the  Lagrangian 
function  L  =  T—  V,  let  us  separate  the  terms  of  the  form 

Q  =  CJff  +  CJ$  +  C^  +  &c. 

In  the  remaining  terms  put  6'  =  6m  V—  1,  <£>'  =  <t>m  V—  1,  and 
so  on,  and  write  down  the  discriminant  If  the  system  oscillates 
about  a  position  of  equilibrium  the  terms  represented  by  Q  are 
absent  and  the  discriminant  thus  formed  will  be  the  determinantal 
equation  giving  the  required  values  of  m.  But  if  the  system 
oscillate  about  a  state  of  steady  motion  wre  must  modify  the  dis- 
criminant by  adding  some  quantity  derived  from  Q  to  each  term. 
To  find  this,  write  above  the  columns  ff  ,  <£>',  &c.  and  before  the 
rows  0,  </>,  &c.  Consider  any  term,  say  the  term  in  the  column 
<fi  and  row  0.  We  must  add  to  that  term  (  C12  —  (721)  m,  where 
(712  —  C2l  is  the  excess  of  the  coefficient  of  (/>'&  above  the  coefficient 
of  <£#'  in  the  expression  for  Q.  Since  the  determinant  is  un- 
changed by  writing  —  m  for  m,  we  may,  if  preferred,  add  the 
excess  of  the  coefficient  of  6'<f>  above  the  coefficient  of  Qfi  provided 
we  adhere  to  one  order  throughout.] 

6.  If  in  the  determinantal  equation  we  write  —  m  for  m,  the 
rows  of  the  new  determinant  are  the  same  as  the  columns  of  the 
old,  so  that  the  determinant  is  unaltered.  When  expanded,  we 
shall  have  an  equation  which  contains  only  even  powers  of  m. 

The  condition  of  dynamical  stability  is  that  the  roots  of  this 
equation  should  be  all  of  the  form 


If  we  write  mz=—p2,  the  roots  of  the  transformed  equation 
must  be  all  real. 

In  the  case  of  equal  roots  of  the  form  m=  +  /3V—  1,  it  has 
been  shown  in  Art.  5  of  Chap.  I.  that  it  is  necessary  for  stability 
that  the  proper  number  of  minors  in  this  determinant  should 
vanish.  If  there  be  two  equal  roots,  these  roots  must  make  all 
the  first  minors  zero  ;  if  three  equal  roots,  all  the  first  and  second 
minors  must  vanish,  and  so  on. 

7.  When  the  system  depends  on  many  co-ordinates  the  labour 
of  expanding  this  determinant  is  often  considerable.  Methods  of 
evading  this  in  certain  cases  will  be  given  in  the  following  chap- 

K.  A.  4 
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ters.     But  when  the  development  is  necessary  we  may  proceed  in 
the  following  manner.     Let  the  determinant  be  written 


&c. 


&c. 


&c. 


and  let 


F  =  — 

•C        — 


qp' 


We  know  that  the  determinant  when  expanded  is  of  an  even 
order,  hence  all  odd  powers  of  m  must  finally  vanish.  Let  us 
expand  the  determinant  in  powers  of  the  F's.  The  first  term  is 
the  discriminant  of  Bm*  +  a,  this  term  is  independent  of  the  F's. 
The  terms  which  contain  the  first  powers  of  the  F's  are  obtained 
by  erasing  any  one  line  of  this  discriminant  and  replacing  it  by 
the  corresponding  F  terms.  But  these  terms  all  vanish  and  we 
need  not  describe  them  minutely.  The  terms  containing  the 
products  and  squares  of  the  F's  may  be  obtained  by  erasing  every 
two  lines  of  the  discriminant  and  replacing  them  by  the  corre- 
sponding F  terms.  Thus  if  we  erase  the  two  first  lines  we  have 
the  determinant 


0 


&c. 


0 

* 
&c. 


&c. 
&c. 
&c. 


&c. 


and  so  on  for  all  the  other  rows  taken  two  and  two.  The  terms 
which"  contain  the  cubes  and  all  odd  powers  of  the  F's  vanish, 
while  the  terms  which  contain  the  fourth  powers  may  be  obtained 
by  erasing  four  lines  of  the  discriminant  and  replacing  them  by 
the  corresponding  F's. 

When  the  determinant  has  been  expanded,  we  have  an  equa- 
tion of  an  even  order  to  find  the  values  of  m.  We  may  therefore 
employ  the  short  method  of  Art.  5,  Chap,  m.,  to  obtain  the 
Sturmian  functions. 

8.  Necessary  and  sufficient  tests  of  the  stability  of  the  motion 
of  a  system  of  bodies  are  given  in  the  preceding  pages.  But  it 
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is  assumed,  as  explained  in  Art.  2,  that  the  co-ordinates  have  been 
properly  chosen.  They  are  supposed  to  have  been  so  chosen  that 
the  coefficients  in  the  expanded  Lagrangian  function  are  all 
constants.  When  this  is  not  the  case  we  must  discover  the 
proper  co-ordinates  to  which  the  system  must  be  referred  before 
we  can  apply  the  test  of  stability.  But  when  the  motion  is  steady 
this  is  not  difficult. 

There  are  obviously  many  such  systems  of  co-ordinates,  and 
one  set  may  generally  be  found  by  a  simple  examination  of  the 
steady  motion.  If  there  are  any  quantities  which  are  constant 
during  the  steady  motion,  they  will  often  serve  for  some  of  the 
co-ordinates.  Others  may  be  found  by  considering  what  quan- 
tities appear  only  as  differential  coefficients  or  velocities.  Practi- 
cally these  will  be  the  most  convenient  methods  of  discovering 
proper  co-ordinates,  since  no  further  change  will  then  be  necessary 
and  we  may  at  once  form  the  determinant  of  stability.  But  if 
these  methods  fail  we  may  adopt  the  following  analytical  method 
of  transforming  (where  possible)  the  general  Lagrangian  function 
with  variable  coefficients  into  one  with  constant  coefficients. 

9.     Let  the  Lagrangian  function  be 

L  =  LQ  +  Afi  +  A£  +  &c.  +  B/  +  B$  +  &c. 


where  the  coefficients  are  all  functions  of  t  and  the  co-ordinates 
0,  (f>,  &c.  have  been  so  chosen  as  to  vanish  along  the  steady  motion. 
We  have  therefore  for  the  steady  motion 

55'-4=o. 

&c.=0. 

The  oscillations  about  the  steady  motion  are  given  by  the 
terms  of  the  second  order.  Our  present  object  is  to  transform 
these  to  others  with  constant  coefficients  by  the  following  substi- 
tutions : 

6  =  p^  +  pjj  +  p3z  4-  &c., 


&c.  =  &c., 
where  the  p's,  qs,  &c.  are  functions  of  the  time  at  our  disposal. 

^  Substituting  and  equating  the  coefficients  of  a/2,  y'2,  &c.   to 
unity,  we  have  as  many  equations  of  the  form 

...  =  l  ...............  (1) 

4—2 
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as  there  are  co-ordinates.     Equating  the  coefficients  of  the  pro- 

m  _  1 
ducts  x'y  ',  x'z't  &c.  to  zero  we  get  n  —  ^—  equations  of  the  form 


...-Q  .........  (2), 

supposing  that  there  are  n  co-ordinates. 

Equating  to  the  constants  alt  or2,  ...  the  coefficients  of  xx'y  yy  , 
&c.  having  subtracted  the  differential  coefficients  of  (1)  we  have 
n  equations  of  the  form 


Adding  the  coefficients  of  xy    and  x'y  and   subtracting   the 

n—  I 
differential  coefficients  of  (2)   we  have    n  —  -  —  equations  of  the 

form 

'  ' 

Equations  (1),  (2)  and  (4)  give  n2  equations  to  find  the  nz 
quantities  p^p^  &c.,  q^2,  &c.  The  solution  of  these  equations  is  a 
purely  geometrical  problem.  If  we  construct  the  two  quadrics 


and  refer  them  to  their  common  conjugate  diameters,  by  writing 


&c.  =  &c., 

making  the  first  quadric  to  become  what  we  may  call  a  sphere  by 
projection;  the  values  ofp&p  &c.,  q^  &c.  thus  found  are  the  values 
required  to  make  some  of  the  coefficients  in  the  Lagrangian 
function  become  constant.  These  values  must  of  course  make  all 
the  other  coefficients  of  the  second  order  in  the  Lagrangian  function 
constants  also,  and  thus  we  have  n(n+1)  analytical  conditions 
that  the  motion  should  be  steady. 

It  might  be  supposed  that  greater  generality  would  be  obtained 
by  replacing  the  zero's  of  equations  (2)  and  (4)  or  the  unities  of 
(1)  by  arbitrary  constants.  This  may  be  convenient  in  practice, 
but  as  we  know  that  by  a  subsequent  real  change  with  constant 
values  of  Pip#  &c.,  q^qv  &c.,  we  can  render  them  zero  or  unity,  it 
simplifies  the  argument  to  perform  the  two  transformations  at 
once. 
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10.  The  geometrical  problem  just  alluded  to  admits  of  a  real 
solution  whenever  one  quadric  can  be  projected  by  a  real  projection 
into  a  sphere.  The  problem  then  becomes  that  of  finding  the 
principal  axes  of  the  other.  This  is  just  our  case,  since  the  ex- 
pression 


is  necessarily  positive  for  all  values  of  &<$ . 

It  is  unnecessary  to  describe  here  the  mode  of  solving  this 
problem.  It  is  sufficient  to  say  that  it  may  be  reduced  to  the 
solution  of  the  symmetrical  determinantal  equation 

=  0, 


where  Dn,  D12,  &c.  are  the  coefficients  of  $2,  0$,  &c.  in  the  second 
of  the  quadrics.  The  roots  of  this  equation  are  known  to  be  real 
when  the  suppositions  just  mentioned  are  satisfied. 

11.  In  order  to  examine  the  fundamental  determinant  in 
Art.  4  a  little  more  closely,  let  us  suppose  it  reduced  to  depend  on 
three  co-ordinates.  We  may  then  have  the  advantage  of  a  geo- 
metrical analogy.  Let  the  co-ordinates  be  £  77,  £>  and  let  the 
equations  of  motion  be  written 

d2 
kJ-A)?t' 


Let  a  geometrical  point  P  move  in  space  so  that  its  co-or- 
dinates referred  to  any  axes  are  f,  97,  f.  Then  the  position  and 
motion  of  the  point  exactly  give  us  the  position  and  motion  of 
the  system. 

12.  Looking  at  the  equations  of  motion  just  written  down  we 
see  that  they  are  similar  to  those  which  give  the  oscillations  of  a 
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system  about  a  position  of  equilibrium,  but  that  there  are  in 
addition  terms  E  -~,  F-^,  &c.  The  general  effect  of  these  terms, 

as  will  appear  from  what  follows  in  the  subsequent  chapter,  is  to 
increase  the  stability.  If  we  transpose  these  terms  to  the  other 
sides  of  the  equations,  we  may  regard  them  as  impressed  forces 
acting  on  the  system,  whose  resolved  parts  in  the  directions  of  the 
axes  f,  77,  f,  are 


v      F  r 

Y=Edt~GTt 


dt         dtl 


We  see  at  once  that 


dt 


dt 


dt 


so  that  these  forces  are  at  once  orthogonal  to  the  path  of  the 
representative  point  P  and  also  orthogonal  to  the  straight  line 
whose  direction  cosines  are  proportional  to  Et  F,  G.  These  forces 
are  therefore  of  the  nature  of  centrifugal  forces,  as  if  they  were 
produced  by  the  rotation  of  the  system  about  this  straight  line. 

13.  We  may  show  that  the  straight  line  (Ey  F,  G)  is  fixed  in 
space.  To  prove  this,  let  us  transform  our  co-ordinates  from  f  ,  77,  £ 
to  x,  y,  z,  where  x,  yy  z  are  connected  with  f  ,  77,  £  by  any  linear 
relations,  such  as 


Let  the  portion  of  the  Lagrangian  function  under  consideration 
(Art.  1)  be 

" 


Substituting  for  f,  77,  %  their  values  in  terms  of  x,  y,  z,  we  find  that 
the  difference  between  the  coefficients  of  xy  and  x'y  is 

G'=G 
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with  similar  equations  for  E'  and  F'.     But  if  /A  be  the  determi- 
nant of  transformation 


with  similar  equations  for  x  and  y.  The  ratios  of  J£,  jPJ  (7  are 
therefore  transformed  as  if  they  were  co-ordinates.  If  the  trans- 
formation be  a  real  transformation  of  Cartesian  co-ordinates,  let 
lengths  each  equal  to  unity  be  measured  from  the  origin  along  the 
axes  Og,  Orj,  0£,  thus  forming  a  tetrahedron  whose  volume  is  V. 
Let  a  similar  construction  be  made  for  the  new  axes,  forming  a 

V 
tetrahedron  of  volume  V.     Then  *  /*  =  -^ .     Hence  the  quantities 

E     F     G 

y,    y,    y  may  loe  transformed  as  if  they  were  lengths  measured 

W       Tf"      C1' 

along  the  axes  and  become  yf ,  -y, ,  -y, .  If  both  systems  of  co- 
ordinates are  rectangular  we  have  V=  ~  ,  V  =  ~  . 

o  o 

W      W      C* 

Let  ft>  be  the  resultant  of  -y,    y,    -~ ,  then  &>  may  "be  regarded 

as  a  fixed  length  measured  from  the  origin  along  a  straight  line 
fixed  in  space. 

Let  v  be  the  velocity  of  the  representative  particle,  6  the  angle 
between  the  direction  of  this  velocity  and  the  axis  whose  direction 
cosines  are  proportional  to  E,  F}  G.  Then  the  resultant  of  the 
forces  X,  F,  Z  is  easily  seen  to  be  2v  Fo>  sin  6  acting  perpendicular 
to  the  axis  and  to  the  direction  of  the  motion.  We  might  call  the 
straight  line  (EFG]  the  axis  of  the  centrifugal  forces. 

[*  Let  (£x  i7x  fj),  (£2  772  £,),  (£3 173  f3)  be  the  co-ordinates  of  three  points  A,  B,  C  re- 
ferred to  any  oblique  co-ordinates.  Let  us  find  the  volume  V  of  the  tetrahedron 
of  which  these  and  the  origin  are  the  angular  points.  Since  the  volume  vanishes 
when  any  angular  point  as  C  lies  in  the  plane  containing  the  origin  and  the  other 
two  A,  B,  the  expression  for  the  volume  must  contain  the  factor 


The  volume  is  evidently  an  integral  rational  function  of  the  co-ordinates  when 
the  axes  are  rectangular  and  the  plane  AOB  is  taken  as  the  plane  of  xy,  it 
easily  follows  that  this  is  true  for  all  axes.  Since  this  function  cannot  be  of  more 
than  the  third  order,  we  have  V'  =  M/j.,  where  M  is  independent  of  the  co-ordinates 
of  A,  B,  C.  When  the  points  A,  B,  G  are  on  the  axes  at  unit  distances  from  the 
origin,  let  Fbe  the  volume  of  the  tetrahedron.  In  this  case  /*=!,  and  .-.  M=  V. 
We  have  therefore  in  all  cases  V  =  T>. 

In  the  text,  let  the  extremities  of  the  unit  lengths  measured  along  the  axes  of 
x,  y,  z  be  called  At  B,  C.  Then  the  (£??£)  co-ordinates  of  A,  B,  G  are  (av  a2,  a8), 
(&i»  &2>  &s)»  (cx>  c2,  c3),  respectively.  Hence  by  what  has  just  been  said  F'=F/x.] 
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14.  The  expressions  for  the  co-ordinates  in  terms  of  the  time 
will  in  general  contain  as  many  periodic  functions  as  there  are 
co-ordinates.  If  the  initial  conditions  are  such  that  each  contains 
one  and  the  same  periodic  function,  the  motion  recurs  after  a  con- 
stant interval  and  the  system  is  said  to  be  performing  a  simple  or 
harmonic  oscillation. 

If  the  system  be  oscillating  about  a  position  of  equilibrium, 
with  a  Lagrangian  function 


we  know*  that  the  harmonic  oscillations  are  represented  by  rectili- 
near motions  of  the  representative  particle,  and  that  these   are 
along  the  common  conjugate  diameters  of  the  two  quadrics  : 
equations  are 


where  a  and  &  are  two  constants  chosen  to  make  the  quadrics  real. 
Let  us  consider  what  are  the  harmonic  paths  of  the  representative 
point  when  the  system  is  oscillating  about  a  state  of  steady  motion. 

In  any  harmonic  vibration   we   have   f  =  L  cos  (\t  +  a)    with 
similar  equations  for  rj  and  f.     Hence 


Substitute  these  in  the  equations  of  Art.  11.  Differentiate  and 
substitute  again.  Multiply  by  f,  77,  £,  add  and  integrate,  we  obtain 

f  \       2          /  I" 

£  /  Y        2 

where  c  is  some  constant.  The  harmonic  path  lies  on  this  quadric, 
which  has  a  common  set  of  conjugate  diameters  with  the  two 
quadrics  a  and  b. 

If  we  resume  the  result  of  the  substitution  of  -p,  &c.  in  the 

equations  of  Art.  11,  and  multiply  by  E,  F,  G  respectively  and  add, 
we  obtain 


[(BnE+  B^F  +  B^G]  %  +  &c.]  X2  +  [(AUE+  ...)£  +  &c.]  =  0, 


*  A  short  paragraph  in  Thomson  and  Tait's  Natural  Philosophy,  page  273,  is 
the  only  notice  of  this  which  the  author  has  discovered. 
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which  is  a  plane,  and  is  diametral  to  the  straight  line  (EFGf)  with 
regard  to  the  quadric  c. 

The  harmonic  paths  are  therefore  ellipses.  The  three  harmonic 
planes  are  diametral  to  the  same  straight  line  and  this  straight 
line  is  fixed  in  space,  being  the  axis  of  the  centrifugal  forces. 

If  we  eliminate  A,  between  the  equations  to  the  plane  and  the 
quadric  c,  we  get  a  cubic  surface  on  which  the  three  harmonic 
conies  lie. 

If  E,  F,  G  are  zero,  which  is  the  case  when  the  system  oscillates 
about  a  position  of  equilibrium,  the  quadric  c  becomes  a  cylinder. 
This  may  be  conveniently  shown  by  referring  the  system  to  such 
co-ordinates  that  the  coefficients  7?12,  J913,  B&,  A^,  A13,  A^  are  all 
zero.  In  this  case  the  diametral  plane  of  every  straight  line  passes 
through  the  axis  of  the  cylinder.  The  harmonic  oscillations  are 
therefore  rectilinear. 

If  R  be  the  length  of  that  semidiameter  of  the  quadric  (c) 
which  is  parallel  to  the  fixed  straight  line  (E,  F,  G\  it  may  be 
shown  that  the 

Product  of  the  axes]  _  R  2c 

of  the  quadric  c     }      JE*  +  F2  +  G2  '  ^  ' 


If  E,  F,  G  are  all  zero,,  and  their  ratio  is  indeterminate,  R  is 
any  diameter.  Hence  one  of  the  axes  of  the  quadric  (c)  must  be 
infinite  and  the  quadric  will  be  a  cylinder. 

[If  the  quadric  (c)  be  a  cylinder  and  E,  F,  G  are  not  all  zero, 
we  must  have  either  A,  zero  or  R  infinite.  In  the  latter  case  the 
axis  of  the  cylinder  will  coincide  with  the  axis  of  the  centrifugal 
forces.] 

The  quadric  (c)  has  also  the  following  geometrical  property. 
Let  the  lengths  of  semidiameters  of  the  quadrics  (a)  and  (b)  drawn 
parallel  to  the  axis  of  the  centrifugal  forces  be  p  and  p.  Through 
the  intersection  of  these  quadrics  describe  a  quadric  so  that  the 


Product  of 
its  axes 


fl 

— 2 2 

\P  P 


+$*+&     /(!_JLVJL__1 

This  quadric  is  similar  to  the  quadric  (c). 

15.  The  introduction  of  the  representative  point  to  exhibit 
the  motion  of  a  system  may  appear  somewhat  artificial.  If  how- 
ever we  properly  choose  the  co-ordinates  the  particle  moves  exactly 
as  a  free  particle,  and  we  might  reduce  the  problem  of  finding  the 
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oscillations  of  a  system  to  a  problem  in  Dynamics  of  a  particle. 
Refer  the  quadric 


to  its  principal  axes  and  let  the  equation  thus  changed  be 

3i'¥+-B«'1-  +  ...  =  &. 

Since  J3U',  B^,  &c.  are  positive  quantities,  we  may  put 

^y,    &c. 


The  quadric  has  thus  been  "projected"  into  a  sphere.  Let  x,  y}  z 
be  now  chosen  as  the  co-ordinates  of  the  system  and  let  the 
Lagrangian  function  be  expressed  in  the  form 


the  terms  of  the  first  degree  being  omitted  as  not  necessary  to  our 
present  purpose.  The  three  equations  of  motion  at  the  beginning 
of  Art.  1  1  take  the  form 


&c.  =  0, 

which  are  the  three  equations  of  motion  of  a  free  single  particle 
whose  co-ordinates  are  x,  y,  z  under  the  action  of  forces  whose  force 
function  is 


and  a  force  acting  perpendicular  to  the  path  and  also  perpendicu- 
lar to  a  fixed  straight  line,  the  force  being  proportional  to  the 
velocity. 

16.  As  an  illustration  of  this  theory,  let  us  here  make  a  short 
digression.  However  the  particles  of  light  may  oscillate,  whether 
in  a  rotatory  or  linear  manner,  we  know  the  motion  is  related  to  a 
certain  plane  called  the  plane  of  polarization.  It  may  be  shown 
that  any  harmonic  oscillation  about  a  position  of  equilibrium  may 
be  represented  by  a  rectilinear  oscillation  of  the  representative 
particle.  Let  us  represent  the  motion  at  any  point  of  the  ether 
by  a  rectilinear  oscillation  in  a  direction  perpendicular  to  the 
plane  of  polarization.  This  would  be  Fresnel's  Vibration.  The 
representative  particle,  as  just  shown,  would  not  necessarily  move 
as  if  it  were  a  free  single  particle.  But  let  us  assume  (and  a  proof 
is  not  necessary  to  our  present  purpose)  that  when  the  oscillation 
is  drawn  as  above  described  the  motion  in  the  plane  of  the  front 


IV.]  THE   REPRESENTATIVE  PARTICLE.  59 

is  the  same  as  that  of  a  free  particle,  while  that  perpendicular  is 
not  free.  On  this  assumption  we  see  that  Fresnel  in  his  theory 
of  double  refraction  is  justified  in  taking  actual  instead  of  relative 
displacements,  for  it  is  the  representative  particle  he  is  considering. 
He  also  neglects  the  force  normal  to  the  front,  for  the  particle 
moves  as  a  free  particle  only  in  the  plane  of  the  front.  These 
general  remarks  are  not  meant  to  explain  Fresnel's  theory,  but 
merely  to  show  how  the  representative  particle  may  be  used  to 
replace  a  complicated  motion. 

17.  [When  a  system  is  performing    a   harmonic  oscillation 
about  a  state  of  steady  motion  or  about  a  position  of  equilibrium, 
the  motion  repeats  itself  continually  at  a  constant  period,  that  is 
to  say,  the  values  of  the  co-ordinates  recur  at  this  interval.     This 
is  the  chief  peculiarity  of  a   harmonic   oscillation.      When   the 
oscillation  is  about  a  position  of  equilibrium,  the  representative 
particle  oscillates  in  a  straight  line  whose  middle  point  repre- 
sents the  position  of  equilibrium.    Thus  the  system  passes  through 
the  position  of  equilibrium  twice  in  each   complete   oscillation. 
When  the  oscillation  is  about  a  state  of  steady  motion  the  path 
of  the  representative  particle  is  an  ellipse  whose  centre  is  at  the 
point  occupied  by  the  system  in  steady  motion  at  the  same  in- 
stant.    Thus  the  system  does  not  in  general  ever  coincide  with 
the  simultaneous  position  of  the  system  in  the  undisturbed  or 
steady  motion.     When  a  system  is  disturbed  by  a  small  impulse 
from  a  state  of  steady  motion,  it  will  in  general  describe  a  com- 
pound oscillation  made  up  of  at  least  two  harmonic  oscillations, 
at  the  instant  of  disturbance  these  two  neutralize  each  other  so 
that  in  the  disturbed  and  undisturbed  motions  two  simultaneous 
positions  are  coincident.     But  it  is  clear  that  this  cannot  occur 
again  unless  either  the  periods  of  the  two  harmonics  are  com- 
mensurable or  the  period  of  one  of  them  is  infinite.] 

18.  [In  some  cases  the  ellipse  degenerates  into  a  straight  line. 
Thus  if  the  quadric  (c)  be  a  cylinder  the  diametral  plane  of  the 
axis  of  the  centrifugal  forces  will  pass  through  the  axis  of  the 
cylinder,  and  thus  the  harmonic  oscillation  corresponding  to  this 
particular  value  of  X  will  be  rectilinear.     In  this  case  the  system 
twice  in  each  oscillation  passes  through  the  position  it  would  have 
occupied  at  the  same  instant  in  the  undisturbed  motion. 

The  quadric  (c)  has  a  common  set  of  conjugate  diameters  with 
the  quadrics  (a)  and  (b).  Hence  if  (c)  be  a  cylinder,  its  axis  must 
be  parallel  to  one  of  the  three  common  conjugate  diameters  of  (a) 
and  (b).  If  we  refer  the  quadrics  (a)  and  (b)  to  their  common 
conjugate  diameters,  they  take  the  form 
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The  cylinder  which  passes  through  their  intersection  and  has 
its  axis  parallel  to  the  diameter  £is  found  by  eliminating  f*  between 
these  equations.  We  see  therefore  that  B^\z  +  A^  =  0.  If  then 
the  axis  of  the  cylinder  cut  the  quadrics  (a)  and  (6)  in  D  and  D' 
respectively,  we  find  that  for  this  oscillation 


2_      _ 

~OD2'    b    ' 

It  has  already  been  shown  that  when  this  value  is  finite,  the 
direction  of  ODD'  is  along  the  axis  of  the  centrifugal  forces.] 

19.  [In  some  cases  two  or  more  of  the  values  of  \  are  zero. 
In  these  cases  the  co-ordinates  will  have  terms  of  the  form  nt  -f  e, 
where  n  and  e  are  two  small  constants.     When,  as  explained  in 
Art.  3  of  this  Chapter,  there  are  several  parallel  states  of  steady 
motion,  these  terms  imply  that  the  motion  is  stable  about  a  state 
of  steady  motion  very  nearly  the  same  as  the  undisturbed  motion 
but  not  coincident  with  it.     The  actual  undisturbed  motion,  unless 
n  is  zero,  is  unstable  in  the  sense  that  if  a  proper  disturbance 
be  given  to  the  system,  the  system  will  depart  widely  from  the 
positions  it  would  have  simultaneously  occupied  in  the  undisturbed 
motion.] 

20.  In  many  cases  of  small  oscillations  it  will  be  found  that 
the  Lagrangian  function  T  —  V  is  not  a  function  of  some  of  the 
co-ordinates  as  0,  $,  &c.  though  it  is  a  function  of  their  differential 
coefficients  6',  $',  &c.     In  such  cases  the  steady  motion  will  be 
usually  given  by  constant  values  of  these  differential  coefficients, 
while  the  other  co-ordinates  as  f,  77,  &c.  are  also  constant.     It  is 
evident  that  the  determinantal  equation  of  Art.  4  is  needlessly 
complicated.     It  is  clear  that  there  will  be  as  many  pairs  of  roots 
equal  to  zero  as  there  are  co-ordinates  0,  (f>,  &c.     It  will  be  an 
advantage  to  eliminate  0',  (/>',  &c.  altogether  from  the  Lagrangian 
function,  and  to  find  the  remaining  roots  by  operating  only  with 
the  co-ordinates  (-,  77,  &c.     We  shall  thus  obtain  a  determinant  with 
just  as  many  rows  as  there  are  co-ordinates  of  the  kind  f  ,  77,  &c. 

Let  LI  be  the  Lagrangian  function  expressed  as  a  function  of 
#'<£',  &c.  £77,  £V,  &c.  Let  L2  be  its  value  when  &$'  are  elimi- 
nated, so  that  L2  is  a  function  of  £77,  %'r)',  &c.  only.  To  effect  this 
elimination  we  have  the  integrals 

dT_         dT_       , 
dO'^C*    d$  = 
where  ci;  C2,  &c.  are  constants.     Then 

dLt_dL       dL,  dff      dL,dJl_ 

d?~~d%  +  dff-di-  '+ 

dL         dff 
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dt  d£       d£      dt  d£'       d£        l  \dt  dj~      d£J 
But  _rf«|._^|  =  o. 

Hence  if  we  take 

dT  dT 

and  eliminate   6',  </>'  by  help  of  the  integrals  -j^-,  =  c  ,     .  ,,  =  ca, 

dc/  d<p 

we  may  treat  L'  Jws^  as  we  do  the  Lagrangian  function  L.  The 
equations  giving  the  small  oscillations  about  the  steady  motion 
will  be 

dt-df7"^"0'  ^c'  =  a 
The  function  L'  may  be  called  the  modified  Lagrangian  function. 

[It  should  be  noticed  that  this  is  equivalent  to  a  partial  use  of 
Hamilton's  transformation  of  Lagrange's  equations.  Sir  W.  R. 
Hamilton  eliminates  all  the  differential  coefficients  6',  <£',  &c.  by 

7/T7  j  rn 

the  help  of  equations  of  the  form  ^  =  u,  --,  =  v,  &c.  where  u,  v, 
&c.  are  made  to  be  new  variables*.  In  our  transformation  only 

*  The  Hamiltonian  transformation  of  Lagrange's  equations  bears  a  remarkable 
analogy  to  the  transformation  of  Eeciprocation  in  Geometry.  This  may  be  shown 
in  the  following  manner. 

When  the  system  has  three  co-ordinates  0,  0,  \f>,  we  may  regard  0',  0',  \f>'  as  the 
Cartesian  co-ordinates  of  a  representative  point  P.  The  position  and  path  of  P  will 
exhibit  to  the  eye  and  will  determine  the  motion  of  the  system.  Let  w,  v,  w  be  the 
Hamiltonian  variables,  so  that 

w=— *      v=^      w  =  —1 

where  Tl  is  the  semi  vis  viva  expressed  as  a  function  of  0,  <f>,  \f/,  0',  <£',  \f/.  Then 
u,  v,  w  may  be  regarded  as  the  co-ordinates  of  another  point  Q  whose  position  and 
path  will  also  determine  the  motion  of  the  system. 

If  the  semi  vis  viva  be  given  by  the  general  expression 


it  is  clear  that  the  point  P  always  lies  on  the  quadric  T±=U  where  U  is  the  force 
function  and  the  co-ordinates  0,  0,  \p  have  their  instantaneous  values.     The  point 
Q  must  therefore  lie  on  another  quadric  which  is  the  polar  reciprocal  of  the  first 
with  regard  to  a  sphere  whose  centre  is  at  the  origin  and  whose  radius  is  equal  to 
X/JJ77.     The  equation  to  the  reciprocal  quadric  is  therefore 
_       1      Ow      v       w      =U, 
2=~2A     u    Ju    A12    A13 
v    Alz    A22    A23 
w    A,*    A*,*    A™  \ 
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those  new  variables  are  introduced  which  would  be  constants  in 
Sir  W.  R.  Hamilton's  transformation. 

This  remark  suggests  an  extension  of  the  process.  If  L  be  a 
function  of  0,  <£,  &c.  as  well  as  of  &y  $ ',  &c.  the  quantities  clt  C2,  &c. 
will  not  be  constants.  We  express  this  by  writing  u,  v,  &c.  instead 
of  ct,  c2,  &c.  Suppose  we  wish  to  eliminate  some  of  the  differential 
coefficients,  viz.  &,  <f>,  &c.  and  to  retain  the  remaining  ones,  viz. 
f ,  if,  &c.  If  we  put 

L'  =  L-u6'-v<l>'-&c. 

we  may  easily  show  as  in  the  preceding  page  that 
rf  dL'     dL!__      & 
dt  d?      d£ 

where  A  is  tlie  determinant,  called  the  discriminant,  which  may  be  formed  from 
the  determinant  just  written  down  by  omitting  the  first  row  and  the  first  column. 

This  is  a  general  expression  for  the  Hamiltonian  function  and  agrees  with  that 
which  may  be  deduced  from  the  result  in  Art.  21,  when  all  the  variables  are  trans- 
formed by  the  Hamiltonian  process. 

Since  the  polar  reciprocal  of  the  polar  reciprocal  is  the  original  quadric,  it 
follows  that 

„,    dT.         ,     dT,        „     dT» 

e  =  d^>    *-•*''    *  =  <^' 
which  are  three  of  the  six  Hamiltonian  equations. 

We  may  also  show  geometrically  that  if  the  coefficients  of  T-.  be  functions  of 

jrr\  :j/p 

any  quantity  6,  then  —~  =  -  -.* .    To  prove  this  we  notice  that  if  x,  y,  z  be  the 

co-ordinates    of    a    point  P,   situated    on    a    radius    vector   OP'  of    a    quadric 

/  0P\a 
0  (x,  y,  z)  =  l  referred  to  its  centre  0,  then  0  (x,  y,  z}=(  -^  J  .    The  quadrics 

T!  =  1  and  T2  =  1  may  be  regarded  as  polar  reciprocals  of  each  other  with  regard 
to  a  sphere  whose  radius  is  *J2  and  whose  centre  is  the  common  centre  of  the 
two  quadrics.  Let  P  be  any  point  on  the  quadric  Tj  =  l,  and  let  the  radius  vector 
be  produced  to  Z  so  that  OP  .  OZ  =  2,  then  the  quadric  T2  =  1  touches  a  plane 
drawn  through  Z  perpendicular  to  OP  and  Q  is  the  point  of  contact.  Let  these 
quadrics  be  slightly  altered  in  consequence  of  a  variation  of  6,  so  that  their 
equations  are  now  T^dT^l  and  Tt  +  dTt  =  l.  Let  OP  and  OQ,  produced  cut 
these  quadrics  respectively  in  P'  and  Q'.  Then 


Now  if  Z'  be  a  point  on  OP  produced  so  that  OP' .  OZf=OP  .  OZt  the  quadric 
T2  +  dT2will  touch  the  plane  drawn  through  Z'  in  some  point  q  near  Q'.  The 
point  Q'  will  therefore  lie  very  nearly  in  the  tangent  plane,  so  that  by  similar 
triangles 

OQ_OZ  =  OP' 
OQ,'     OZ'      OP ' 

Since  each  of  these  ratios  is  indefinitely  nearly  equal  to  unity,  it  follows  that 
dT,=  -dTz. 

If  we  put  L  =  Tt  +  17  and  H=T2-U,  Lagrange's  equations  may  be  written  in 

the  forms  u'  =  -^  ,     v'  =  -=- ,    «>'=  —  .     Hence  we  have 
dd  d<j>  d$ 

,    dH  ,    dH  ,    dH 

-U=d6>     -V=d*>      -W=d*> 

which  are  the  remaining  three  of  the  Hamiltonian  equations. 
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We  have  thus  as  many  equations  of  the  Lagrangian  form  as 
there  are  variables  f  ,  77,  &c.  Also  since  u  =  -^  ,  &c.  we  have  by 
differentiation 


dL'       dL 


d& 


with  similar  equations  for  <£',  &c.     By  Lagrange's  equations  we 
obtain 

dL'       , 


Thus  we  have  as  many  sets  of  equations  of  the  Hamiltonian 
form  as  there  are  variables  6,  <f>,  &c.] 

21.     We  may  effect  this  elimination  once  for  all  and  find  a 
definite  expression  for  L'. 

Let  the  kinetic  energy  be 

m  rrr       "          .       /n         /v    i  / 


Then  the  integrals  used  will  be 

+  ...  =  cl- 


&c.  =  &c. 

For  the  sake  of  brevity  let  us  call  the  right-hand  members  of 
these^  equations  q  -  X,  c2  —  F,  &c.  Since  T  is  a  homogeneous 
function,  we  have 


F)  +  &c. 


If  we  substitute  in  the  second  line  the  values  of  ff,  f  ,&c.  found 
by  solving  the  integrals  just  written  down,  we  have 


-X      T6l 

F/Tl 
J-  0, 
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where  A  is  the  discriminant  of  the  terms  in  T,  which  contain  only 
6 ',  <f>,  &c.,  and  may  be  derived  from  the  determinant  just  written 
down  by  omitting  the  first  row  and  the  first  column. 

We  may  expand  this  determinant  and  write  it  in  the  form 


2A 


-  V 


0       c,        c 

1 

0 

X 

c,       T99      5 

?**    ... 

+  2z\ 

X 

Ts 

c*       Te,      ' 

r*    r 

Y 

Te 

1 

0       X 

Y      . 

A 

c,      Tee 

?:; 

> 

where  X,  Y,  &c.  stand  for 


&c. 

The  first  of  these  three  determinants  will  contain  only  the 
constants  cv  ca,  &c.,  and  the  co-ordinates  f,  ?;,  &c.  The  second  will 
not  contain  cz,  cv  &c.  but  will  be  a  quadratic  function  of  f  ',  77',  &c. 
The  last  determinant  will  contain  terms  of  the  form  f',  rj  with 
variable  coefficients  which  may  also  be  functions  of  cv  c2... 

Since  ff,  7777',  &c.  are  all  small  quantities,  it  is  clear  that  this 
expression  for  L'  when  expanded  will  take  a  form  precisely  similar 
to  that  given  in  Art.  2,  only  that  we  have  fewer  variables  to  deal 
with. 

22.     As  an  example,  let  us  consider  the  following  problem. 

A  body  has  a  point  O  which  is  in  one  of  the  principal  axes  at 
the  centre  of  gravity  G  fixed  in  space.  The  body  is  in  steady 
motion  rotating  with  angular  velocity  n  about  OG  which  is  vertical. 
Find  the  conditions  that  the  motion  may  be  stable. 

Let  OA,  OB,  00  be  the  principal  axes  at  0  and  let  OC  co- 
incide with  the  vertical  OZ  in  steady  motion.  Let  f,  rj  be  the 
direction  cosines  of  the  vertical  OZ  referred  to  OA,  OR  Let 
«!,  o>2,  o>3  be  the  angular  velocities  about  the  principal  axes  at  0. 
Then  to  the  first  order 
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Let  6  be  the  angle  ZO  C,  ty  the  angle  the  plane  ZO  C  makes 
with  a  plane  ZOX  fixed  in  space  and  </>  the  angle  it  makes  with 
the  plane  A  0  C  fixed  in  the  body.  Then 

f  =  —  sin  6  cos  </>) 
97  =     sin  6  sin  </>)  ' 

0)3  =  <£'  -f  ijr'  COS  9 

=  f+(x'-* 

putting  •%  =  <j>  +  -ty-.     We  easily  find 


If  then  A,  B,  C  be  the  principal  moments  of  inertia  at  0,  the 
Lagrangian  function  is 


W  +  x'B'  +    (-f 


where  M  is  the  mass  of  the  body,  and  h  =  00. 

Since  ^  is  absent  from  the  equation  we  have  the  integral 
dL 

_    ^—  /» 

7       /    -    ^1> 

dx 

which  gives 

%  =n  +  terms  of  second  order. 
Hence 


const  +        '2  + 


Using  this  as  the  Lagrangian  function  we  easily  find 

(A-C}n*  +  Mgh-Bm\  (A  +  B-C)nm 

-  (A  +  B  -  0)  ww,  (5  -  0)  n2  +  %A  -  Am 


=0. 


R.A.  5 


66  THE  MODIFIED  [CHAP. 

The  roots  of  this  equation  to  find  m  must  for  stability  be  of 
the  form  ±  jB*J  —  1.  Putting  mz  —  —  X2  we  have  a  quadratic  to 
find  X2.  The  roots  of  this  quadratic  must  be  real  and  positive. 

If  A  =  B,  as  in  the  case  of  a  top  spinning  with  its  axis  vertical, 
we  have 

C*n*  - 


The  motion  is  stable  or  unstable  according  as  (7V  is  greater  or 
less  than  kAMgh.  If  C2n2  =  ^AMgh,  the  equation  has  equal  roots 
and  as  the  first  minors  are  not  zero  the  motion  is  unstable. 

23.  [As  another  example  of  the  use  of  the  modified  Lagrangian 
function,  let  us  consider  a  case  discussed  by  Prof.  Ball  in  the 
Notices  of  the  Royal  Astronomical  Society  for  March,  1877.  In  a 
problem  in  Physical  Astronomy,  we  want  the  relative  co-ordinates 
of  the  system,  while  its  absolute  motion  in  space  does  not  concern 
us.  Lagrange's  equations  involve  both  the  relative  and  absolute 
co-ordinates,  and  are  therefore  not  particularly  well  adapted  for 
such  problems.  By  using  the  modified  Lagrangian  function,  we 
may  eliminate  the  absolute  co-ordinates. 

Let  the  system  have  n  co-ordinates,  let  us  choose  as  three  of 
them  the  co-ordinates  of  the  centre  of  gravity  of  the  whole  system, 
viz.  0,  $,  i|r.  There  will  then  remain  n  —  3  co-ordinates  which  are 
independent  of  these.  Let  T  be  the  kinetic  energy  of  the  system 
relative  to  its  centre  of  gravity,  V  the  potential  energy,  M  the 
whole  mass.  Then  the  Lagrangian  function  is 


In  problems  in  Physical  Astronomy  the  potential  energy  is  a 
function  only  of  the  relative  positions  of  the  bodies,  and  is  there- 
fore independent  of  6,  <£,  ty  and  their  differential  coefficients.  We 
have  therefore 

dL_         dL_         dL  _ 

def~''c*   d<t>/='~c*   d^'~Cy 

Hence  the  modified  Lagrangian  function  is 
L'=  T  -  V-  a  constant. 

It  is  independently  clear  that  we  might  take  this  as  the 
Lagrangian  function,  for  the  first  three  terms  of  L  do  not  enter  into 
any  one  of  the  Lagrangian  equations,  except  the  three  formed  by 
differentiating  with  regard  to  ff,  </>',  ^'. 

The  function  T'  is  made  up  of  two  parts,  (1)  the  kinetic  energies 
of  the  rotations  of  the  bodies  about  their  centres  of  gravity,  which 
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we  may  call  T7/,  and  (2)  the  relative  kinetic  energies  of  the  several 
bodies,  each  collected  at  its  centre  of  gravity,  which  we  may  call 
T2'.  Let  mv  ra2,  &c.  be  these  masses;  xlt  a?2,  &c.  the  abscissae  of 
their  centres  of  gravity  referred  to  the  centre  of  gravity  of  the 
whole  as  origin.  Then,  accents  denoting  differential  coefficients 
with  regard  to  the  time,  we  have 

'  4-  &c.  =  0. 


Let  us  square  this  and  write 

2<<  =  *;24-*;2-(<-<)2, 

If  we  examine  the  coefficient  of  any  power  as  xj*  we  see  that  it 
is 

m?  +  wi1(mg  4-  m3  +  &c.)  =  m£m. 
Hence  the  square  becomes 

—  #'2  —  0. 


Similar  expressions  hold  for  the  y  and  z  co-ordinates.  Hence 
on  the  whole  we  see  that  the  relative  kinetic  energies  of  the  several 
bodies  collected  at  their  respective  centres  of  gravity  is 


22m    ' 

where  v  is  the  relative  velocity  of  the  centres  of  gravity  of  the 
masses  rat,  mz.  If  we  express  this  in  any  kind  of  co-ordinates,  we 
may  use  the  Lagrangian  function  L'  to  find  the  relative  motion. 

The  expression  for  T{  agrees  with  that  given  by  Prof.  Ball,  but 
his  demonstration  is  quite  different.  Prof.  Cayley  has  given 
another  demonstration  in  the  same'number  of  the  Astronomical 
Notices. 

The  Lagrangian  function  thus  found  may  be  still  further 
"  modified."  To  avoid  symbols  of  summation,  let  us  consider  the 
case  of  three  particles  moving  in  one  plane  under  their  mutual 
attractions.  Let  the  separate  masses  be  m1}  mv  ra3,  and  let  //,  be 
their  sum.  Referring  the  system  to  ml  as  a  central  mass,  let  the 
distances  of  w2,  ra3  from  ml  be  respectively  r  and  p,  and  let  the 
opposite  side  of  the  triangle  be  It.  Let  the  interior  angle  be- 
tween r  and  p  be  <£  and  the  exterior  angle  between  r  and  R  be 
X-  Let  6  be  the  angle  r  makes  with  some  fixed  straight  line  in 
space.  We  easily  find 


where  A  =  m^n^r*  +  ??i1map2  -f 


C  =  •=- 


5—2 
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and  v  is  the  relative  velocity  of  the  masses  m,  m'  calculated  on 
the  supposition  that  rax  is  fixed,  and  that  the  straight  line  r  has 
no  rotation  round  m:.  Thus  A,  B,  C  are  all  functions  of  r,  p,  <p 
and  their  differential  coefficients  with  regard  to  the  time. 

If  we  only  want  the  changes  in  the  form  and  magnitude  of 
the  triangle  joining  the  three  particles,  we  may  eliminate  &  by 
means  of  the  equation 

Aff 


We  then  find  as  our  modified  Lagrangian  function 

L-Mc-l^-v 


p\"     2       A 
which  contains  only  the  three  co-ordinates  r,  p  and  <£].. 

24.  When  the  geometrical  equations  contain  differential 
coefficients  of  the  co-ordinates  (• ,  77,  f,  &c.  of  the  system  with  regard 
to  the  time,  we  cannot  express  the  co-ordinates  xy  y,  z  of  any 
element  of  a  body  in  terms  of  f ,  77,  £  &c.  by  means  of  equations  of 
the  form 


y 


It  follows,  as  is  pointed  out  in  our  books  on  Rigid  Dynamics, 
that  Lagrange's  equations  cannot  be  employed  in  the  form 

d   dT_dT__dV 
dt  d?     d%~     d£' 

In  many  of  the  most  interesting  problems  in  Rigid  Dynamics, 
it  so  happens   that  the   geometrical   equations   do   contain   ~t 

-~  ,  &c.     For  example,  let  a  sphere  be  set  rotating  about  a  vertical 

diameter  and  be  on  the  summit  of  a  perfectly  rough  surface  of 
any  form.  If  a  small  disturbance  be  now  given  to  it,  the  sphere 
may  roll  round  and  round  the  summit.  During  this  motion  the 
velocity  of  the  point  of  contact  is  zero,  and  our  mode  of  re- 
presenting this  analytically  in  terms  of  the  co-ordinates  will  give  us 
two  equations  of  the  form 


To  include  such  cases  the  equations  of  motion  must  be  modi- 
fied. If  L  be  the  difference  between  the  kinetic  and  potential 
energies,  all  the  Lagrangian  equations  may  be  written  in  the  form 
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where  8f,  Sij,  &c.  are  any  small  arbitrary  displacements  consistent 
with  the  geometrical  equations.  But  if  these  geometrical  equa- 
tions be  given  in  the  form 

--  =  °1 

h...=OK 
&c.  =  OJ 

these  arbitrary  displacements  must  satisfy 


&c.  =  0 

If  this  were  not  the  case,  the  geometrical  displacement  of  the 
body  given  in  applying  Virtual  Velocities  would  not  be  such  as  to 
cause  the  unknown  frictional  forces,  &c.  to  disappear.  Using  the 
principle  of  Indeterminate  Multipliers,  we  get 

d  dL     dL 
d 


&c.  = 
These  joined  to  the  geometrical  equations 

01F+0,i/+...  =  o' 

&c.  =0 

are  sufficient  to  determine  the  unknown  co-ordinates  f,  97,  &c.  and 
the  multipliers  X,  fj,,  &c. 

It  will  be  more  convenient  to  write  these  equations  in  the  form 
d_dL     dL       dQ        dH 
~~  *** 


ddL_dL        dG       dH 

dtdrj'  a^4A3?+^af  "* 

&c.  =  0, 
the  geometrical  equations  being 

0  =  0,       77=0,       &c. 

It  is  of  course  obvious  that  these  indeterminate  coefficients 
X,  p,  &c.  are  merely  the  frictions  or  other  resistances  introduced 
into  the  equations  in  a  convenient  form. 

25.  In  order  to  apply  these  equations  to  the  oscillations  of  a 
system  about  a  state  of  steady  motion,  it  will  be  convenient  to 
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change  the  co-ordinates  f  ,  77,  &c.  into  others  0,  $,  &c.  which  vanish  in 
the  steady  motion.  Let  L  be  thus  expanded  in  powers  of  6,  <£,  &c. 
as  explained  in  Art.  1,  and  let  P  and  Q  have  the  meaning  given 
to  them  in  Art.  5. 

Let  us  then  put 

^a  +  07,     7/  =  /3  +  <£',     &c. 
X  =  X0+X1,     fA  =  p0  +  fii,     &c. 


where  a,  /3,  \,  pot  &c.  are  the  values  of  f,  77',  X,  ^,  &c.  in  steady 
motion.     The  geometrical  equations  will  then  take  the  form 


G=G 

A      \  ^  i     \'  I       * 

&C.  = 

and  the  equations  connecting  S0,  80,  &c.  will  be 


+  &c.  =  0) 
&c.  =  OJ  ' 


In  these  equations  6^,  6r2,  &c.  are  functions  of  0,  <£,  &c.  Let 
a  square  bracket  indicate  that  the  value  of  the  inscribed  quantity 
in  steady  motion  is  to  be  taken.  Thus  [6rJ  means  the  value  of  6rx 
when  0,  </>,  &c.  have  all  been  put  zero. 

The  equations  of  steady  motion  may  then,  exactly  as  in  Art.  3, 
be  written 


&c.  =  0 

From  these  the  relations  which  exist  between  a,  ft,  &c.,  X0,  yu,0,  &c. 
may  be  found. 


The  equations  of  the  oscillatory  motion  may  be  written 
dP  +  d  VdG~  dH 


j 


with  similar  equations  for  </>,  ^,  &c. 


26.  The  final  determinant  written  for  two  variables  0,  <f>,  and 
two  geometrical  equations  G  and  H  in  the  notation  of  Art.  19, 
will  be 
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u 

«.[§ 


<•• 


u 

[§] 


tc 


-((?„-  eg 


r«*fln 
UfJ 


dG 


dd 


m 


[dH~\ .  \dH~\ 

\_  d9  J      \_d$  J 


m 


It  will  be  noticed  how  very  much  this  determinant  is  simplified 
if  the  values  of  X,  p  in  steady  motion  are  zero. 

27.  Let  us  apply  these  equations  to  the  solution  of  the  follow- 
ing problem. 

A  heavy  sphere  rotating  about  a  vertical  diameter  rests  in  equi- 
librium on  the  summit  of  a  perfectly  rough  surface  and  being 
slightly  disturbed  makes  small  oscillations,  find  the  periods. 

As  the  sphere  moves  about,  its  centre  always  lies  on  a  surface 
which  may  be  called  parallel  to  the  given  surface.  Let  the  high- 
est point  of  this  surface  be  taken  as  the  origin  and  let  the  axes  of 
x  and  y  be  the  tangents  to  its  lines  of  curvature  at  0,  so  that  the 
equation  to  the  surface  in  the  neighbourhood  of  0  is 


Let  P  be  the  centre  of  the  sphere,  PC  that  diameter  which  is 
vertical  when  the  sphere  is  in  equilibrium  on  the  summit.  Let 
PA,  PB  be  two  other  diameters  forming  with  PG  a  system  of 
rectangular  axes  fixed  in  the  sphere.  Let  the  inclination  of  PG 
to  the  axis  of  Z,  which  is  vertical,  be  0,  and  let  the  vertical  plane 
through  PC  make  with  the  plane  xz  an  angle  ^,  and  with  the 
plane  GPA  an  angle  0.  The  vis  viva  2T  of  the  sphere  will  then 
be  given  by 


T=  i  (x*  +  y'2)  +  i  F  {(</>'  + 


cos  ey  +  <9'2  + 
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Let    sin  6  cos  ty  =  £,     sin  9  sin  ^  =  77,     then    we   have  to  the 
necessary  degree  of  approximation 


Also  let  <f>  +  ^=X'  These  transformations  of  co-ordinates  are 
all  permissible,  because  they  do  not  involve  any  differential  coeffi- 
cients with  regard  to  the  time.  We  thus  find  if  L  be  the  differ- 
ence between  the  kinetic  and  potential  energies 


If  (ox)  cov,  a),  are  the  angular  velocities  of  the  sphere  about 
parallels  to  the  axes,  the  geometrical  conditions  are 


x  —a    ay  —  co,  —  }  =  0 


where  a  is  the  radius  of  the  sphere.     These  equations  by  well- 
known  rules  reduce  to 

-  -  +  <f>'  sin  ty  sin  6  +  &  cos-^  -  (^  +  </>'  cos  6)  ^  =  0 

-^-(^'cosiJrsintf  +  fl'sin^  +  ^'H-  </>'cos0)  -  =  0 
expressing  these  in  terms  of  our  new  co-ordinates  we  have 


The  position  of  the  system  has  now  been  expressed  in  terms 
of  such  co-ordinates,  that  the  coefficients  in  the  governing  func- 
tions L,  G,  H  are  all  constant.  See  Art.  2. 

The  steady  motion  is  given  by  x,  y,  f,  77  all  zero,  and  %  =  n. 
To  find  X0,  fju0  we  may  use  the  equations  of  steady  motion 


where  q  stands  for  any  one  of  the  co-ordinates.      Taking  q  =  x  and 
q  =  y,  we  see  that  \  =  0  and  fiQ  =  0. 

To  find  the  oscillation  we  may  substitute  in  the  determinant 
and  thus  form  the  equation  which  gives  the  periods.     As  most  of 
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the  constituents  of  the  determinants  are  zero,  it  will  be  more  con- 
venient to  form  each  equation  directly  from  the  standard  formula 


where  q  stands  for  any  one  of  the  co-ordinates.     Taking  q  in  turn 
to  be  x,  y,  %,  f,  77  we  find 


Pi 

x 

p, 


Putting  x  =n  these  with  the  two  geometrical  equations  are  all 
linear  and  ready  for  elimination. 

Eliminating  \,  ^  we  have 


Substituting  for  f ,  ij  from  the  two  geometrical  equations,  we 
have 

Jc*  4-  a2   „        a?  _      w&2    ,  n 


To  solve  these  put  se  =  JTcos  (_pi  +  ^)  y  =•  Ysin  (pt  +  q)  so  that 
p  is  the  quantity  required.     We  obviously  have 


which  is  a  quadratic  to  find  p*. 

If  plt  p2  have  opposite  signs  the  roots  cannot  be  real,  and  the 
steady  motion  must  be  unstable.  If  plt  pz  are  both  positive,  so 
that  the  sphere  is  on  the  summit,  the  motion  is  stable  only  if 


CHAPTEK    V. 


Certain  subsidiary  determinants  are  formed  from  the  dynamical  de- 
terminant, and  it  is  shown  that  there  must  be  at  least  as  many 
roots  inaicating  stability  as  there  are  variations  of  sign  lost  in 
t/iese  subsidiary  determinants,  and  must  exceed  the  number  lost  by 
an  even  number.  Arts.  1  —  5,  and  9. 

This  is  equivalent  to  a  maximum  and  minimum  criterion  of  stability 
with  similar  limitations.  Arts.  6  —  8. 

Effect  of  equal  roots  on  this  test  of  the  stability  of  the  system.  Arts. 
10,  11. 

Example.     Art.  12. 

1.  In  order  to  test  the  nature  of  the  roots  of  the  determi- 
nantal  equation,  let  us  apply  a  method  analogous  to  that  by  which 
Dr  Salmon  in  his  Higher  Algebra  proves  the  reality  of  the  roots 
of  the  equation  which  occurs  in  the  determination  of  the  secular 
inequalities  of  the  planets. 

2.  Let  A  be  the  determinant  which  forms  the  left-hand  side 
of  the  fundamental  equation,  let  AOT  be  the  determinant  formed 
by  omitting  the  j?th  row  and  qih  column.     Let  Ar  be  the  determi- 
nant formed  by  omitting  the  first  r  rows  and  r  columns.     Thus 
Ax  =  Au.     We  then  have  by  a  known  theorem  in  determinants 


It  has  already  been  noticed  that  if  we  change  m  into  —  m,  the 
determinant  A  is  changed  into  another  determinant  whose  columns 
and  rows  are  the  same  as  the  rows  and  columns  of  the  first  de- 
terminant. It  easily  follows  that  the  minor  —  A12  is  changed  into 
the  minor  —  A21  by  changing  m  into  —  m.  Hence  if 


then  A21  =  0  (m2)  -  mf  (m2) 


) 

}  ' 
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Hence  the  product  A12A21  is  necessarily  positive  for  all  negative 
values  of  m*. 

It  also  follows  that  if  A12  vanishes  for  any  negative  value  of  mz 
then  A21  also  vanishes  for  the  same  value  of  w2. 

3.     When  the  determinant  At  vanishes,  we  have 


so  that  A  and  A3  must  have  opposite  signs,  or  one  of  them  must 
be  zero.  Consider  then  the  series  of  determinants 

A,  A,,  A2,  A3,  ...... 

each  one  being  formed  from  the  preceding  by  erasing  the  first 
row  and  the  first  column.  We  thus  have  a  series  of  functions 
of  m2  whose  degrees  regularly  diminish  from  the  nih  to  the  first. 
As  we  may  suppose  the  determinant  A  to  have  a  row  and  a  column 
of  zeros  added  on  at  the  bottom  and  right-hand  side,  but  with  any 
positive  constant  in  the  right-hand  bottom  corner,  we  may  add 
to  this  series  of  determinants  any  positive  constant.  We  have  just 
proved  that  if  any  determinant  of  this  series  vanish  for  a  nega- 
tive value  of  m2,  the  two  determinants  on  each  side  have  opposite 
signs.  The  case  in  which  two  successive  determinants  vanish  for 
the  same  value  of  m2  will  be  considered  afterwards. 

We  may  then  use  these  determinants  in  a  manner  somewhat 
similar  to  that  in  which  we  use  Sturm's  functions,  provided  no  two 
successive  functions  vanish  for  the  same  negative  value  of  m2. 
No  variation  of  sign  can  be  lost  as  we  pass  from  m2  =  —  oo  to 
m2  =  0  except  by  the  vanishing  of  the  determinant  A  at  the  head 
of  the  series.  And  when  a  variation  of  sign  is  lost,  it  will  be 
regained  again  at  the  next  root,  unless  a  root  of  the  determinant  Ax 
separates  the  two  roots  of  the  determinant  A.  If  therefore  in  this 
passage  from  —  oo  to  zero,  as  many  variations  of  sign  are  lost 
as  is  indicated  by  the  highest  power  of  m2,  the  values  of  m2  found 
from  the  determinantal  equation  must  be  all  real  and  negative. 
It  will  also  follow  that  the  roots  of  each  of  the  series  of  determi- 
nants are  all  real  and  negative,  and  that  the  roots  of  each  separate 
the  roots  of  the  determinant  next  above  it.  If,  however,  the 
proper  number  of  variations  of  signs  be  not  lost  in  the  passage 
from  m2  =  —  oo  to  m2  =  0,  it  does  not  follow  that  the  values  of  m2 
are  not  real  and  negative. 

4.  If  the  proper  number  of  variations  of  sign  has  not  been 
lost  in  this  passage  from  m2  =  —  oo  to  m2  =  0,  this  proposition 
does  not  leave  us  without  information  as  to  the  nature  of  the  roots. 
We  infer  that  the  number  of  real  negative  values  of  m2  is  equal 
to  or  exceeds  the  number  of  variations  of  sign  lost  by  an  even 
number  and  unless  the  number  of  variations  be  even  the  system 
is  unstable. 
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If  we  do  not  object  to  the  labour  of  expanding  the  determi- 
nants, we  might  extend  this  theorem  to  determine  the  positions 
of  the  negative  values  of  m2  as  well  as  their  number.  The  number 
of  real  negative  values  of  m2  between  m2  =  —  a  and  m2  =  —  <3  is 
equal  to,  or  exceeds  by  an  even  number,  the  number  of  variations 
of  sign  lost  in  the  series  of  determinants.  In  this  form  the 
theorem  resembles  Fourier's  theorem  in  the  Theory  of  Equations. 

5.  The   converse   of  this  proposition   has  not  been  proved. 
If  in  the  passage  from  m2  =  —  GO    to  m2  =  0,  the  number  of  varia- 
tions of  signs  is  unaltered,  it  is  not  true  that  the  values  of  m2 
cannot  be  real  and  negative.     Thus  in  the  simple  case 

=  0, 

where  an,  a12,  a^  are  all  positive,  no  variations  are  lost,  yet  if 
a12  >  Van  +  Va^  the  values  of  m2  are  real  and  negative.  And  if 
at2  =  Van  +  VaM  the  roots  are  equal  and  negative.  It  will  be 
noticed  that  the  minors  in  this  last  case  are  not  zero. 

6.  In  order  to  discover  the  meaning  of  these  losses  or  gains 
of  changes  of  sign,  it  will  be  convenient  to  make  such  changes  of 
the  co-ordinates  as  will  simplify  the  dynamical  determinant  as 
much  as  possible.     Let  us  write 


If  we  now  change  our  co-ordinates  by  writing  for  6,  <£,  &c. 
linear  expressions  of  some  new  co-ordinates,  we  know  that  we  can 
clear  this  expression  of  all  the  terms  containing  the  products. 
We  know  also  that  this  can  be  done  in  an  infinite  number  of 
ways.  We  thus  have 


where  the  symbols  6lt  c^,  &c.  represent  the  new  co-ordinates. 
Again,  let  us  consider  the  expression 


the  coefficients  are  here  the  values  of  P,  Q,  &c.  in  the  general 
expression  for  T,  Art.  1,  Chap.  IV.  when  6,  </>,  &c.  are  all  put  zero. 
But  since  T  is  necessarily  positive  for  all  values  off,  ?/...f,  ??,... 
it  follows  that  WQ  is  positive  for  all  values  of  & ',  <f>  ...  Hence 
by  a  well-known  theorem,  we  may  by  a  real  linear  transform- 
ation of  the  variables  clear  the  expression  WQ  also  of  the  terms 
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containing  the  products  6',  </>',  £c.  and  can  make  the  coefficients 
of  the  squares  any  positive  constants  we  may  please.  We  thus 
have 


It  may  be  shown  that  we  cannot  in  general  clear  the  expres- 
sion for  jTof  all  the  terms  containing  0<j>',  0'$,  &c.  unless  Ca—C9l=0, 
&c.  by  substituting  for  (9,  </>,  &c.,  any  linear  functions  of  other 
variables.  As  this  is  only  a  negative  result  of  which  no  further 
use  will  here  be  made,  it  is  unnecessary  to  supply  the  demon- 
stration. 

7.  Using  these  simplifications  the  determinant  A  will  now 
take  the  simpler  form, 

A  = 


,=0, 

—  a12m,      m*  +  a22! 
-  am,       -  a23m, 


where  «12  =  <721  -  Cw  &c. 

If  we  form  from  A  the  series  of  subsidiary  determinants 
A,  A,,  A,..., 

terminating  with  any  positive  constant,  we  see  that  when  m2  =  —  oo, 
these  subsidiary  determinants  are  alternately  positive  and  negative, 
and  when  ra*  =  0,  they  become 


which  are  all  positive  if  a^,  aM...  are  all  positive.  Hence  if 
au,  «22...  are  all  positive,  the  proper  number  of  changes  of  sign  has 
been  lost,  and  therefore  the  roots  of  the  dynamical  equation  A  =  0 
satisfy  the  condition  of  stability. 

If  an,  a22,  &c.  are  all  positive,  we  see  that  F0  is  a  minimum  for 
all  variations  of  0lt  ^>x,  &c.  and  therefore  for  all  variations  of  the 
original  co-ordinates.  If  au,  «22,  &c.  are  not  all  positive,  there  will 
be  as  many  variations  of  sign  lost  as  there  are  positive  quantities 
in  the  series  an,  or22,  &c.  In  this  case  Ve  is  a  minimum  for  some 
variations  of  6,  <£,  &c.  and  not  for  others.  It  is  shown  in  the 
appendix  to  Williamson's  Differential  Calculus,  that  n  independent 
conditions  are  necessary  that  a  quadratic  expression  of  n  variables 
should  be  always  positive.  These  are  given  in  the  form  of  deter- 
minants and  may  be  briefly  summed  up,  in  the  statement  that 


78  MAXIMUM  CRITERION   OF   STABILITY.  [CHAP. 

is  always  positive  if  the  roots  of 

A19,       &c.    =0, 
An,       &c. 

&c.  &c.  &c.        &c. 

are  all  real  and  negative. 

8.  We  may  now  put  the  proposition  of  Art.  3  into  another 
form.     Let  L  be  the  general  expression  for  the  excess  of  the  kinetic 
energy  over  the  potential  energy  of  a  dynamical  system  in  terms  of 
its  n  co-ordinates  x,  y,  &c.     Let  this  system  be  moving  in  steady 

motion  with  constant  values  of  TT  ,  TT  ,  &c.     Then  if  L  be  a  maxi- 

d  u    dtj 

dx    (lv 

mum  for  all  variations  o/x,  y,  &c.  keeping  -^- ,  -£  ,  &c.  unchanged, 

d  tj     d  u 

then  that  steady  motion  is  stable. 

If  however  only  r  of  the  n  conditions  necessary  to  make  L  a 
maximum  be  satisfied,  then  the  number  of  roots  of  the  dynamical 
equation  which  satisfy  the  conditions  of  stability  is  equal  to  r  or 
exceeds  r  by  an  even  number.  There  cannot  be  stability  unless  n  —  r 
is  an  even  number. 

If  the  system  be  oscillating  about  a  position  of  equilibrium, 

fl  'T1         CM  ?/ 

-jT ,  Hr ,  &c.  are  all  zero,  and  this  leads  at  once  to  the  condition, 
at    at 

that  the  equilibrium  is  stable  if  the  potential  energy  is  a  minimum. 

9.  In  this  reasoning,  we  have  for  convenience  excepted  the 
case  in  which  two  successive  determinants  in  the  series  A,  Ax,  A2... 
vanish  for  the  same  value  of  m2.     But  this  exception  is  of  no  real 
importance,  for  we  may  change  these  determinants   into  others 
whose  constituents  are  very  slightly  different  from  those  of  the 
given  determinants  but  which  are  such  that  no  successive  two  of 
the  series  have  a  common  root.     In  the  limit,  therefore,  when 
these  arbitrary  changes  of  the  constituents  are  indefinitely  small, 
the  roots  of  the  series  of  determinants  will  still  be  real  under  the 
same  circumstances  as  before,  and  the  roots  of  each  will  separate, 
or  coincide  with,  the  roots  of  the  next  above  it  in  the  series. 

To  show  that  these  changes  are  possible,  let  us  consider  the 
row  of  determinants  beginning  at  the  last.  The  determinant  An  is 
a  positive  constant,  the  next  An_t  is  m2  +  ftnn.  Proceeding  thus, 
suppose  we  arrive  at  two  determinants  which  we  may  call  Az  and  A 
which  have  a  common  root.  If  we  now  change  the  constituents 
On,  a^,  a13,  &c.  into  au  +  Soj,,  <z124  Sa12,  &c.  we  do  not  alter  A1?  but, 
except  for  the  root  w2  =  0,  we  do  alter  A  in  an  arbitrary  manner. 
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When  for  example  a,3  is  altered,  we  alter  a  constituent  both  in  the 
first  row  and  in  the  first  column.     Since 


where  A'  is  the  determinant  formed  from  A  by  omitting  the  first 
and  third  rows  and  columns,  we  see  that  when  A  and  Aj  both 
vanish,  the  product  A13A31  and  therefore  by  Art.  2  both  A]3  and 
A31  must  vanish.  The  determinant  A  is  therefore  altered  by 
-  A2  (&a13)2  which  does  not  vanish,  since  A2  is  by  hypothesis  finite 
for  the  particular  value  of  m2  under  consideration. 

If  any  determinant  of  the  series  vanishes  when  ra2  =  0,  it  is 
clear  that  one  of  the  quantities  an,  a&...  must  be  zero.  If  we 
replace  this  by  any  small  positive  quantity,  the  argument  will 
apply  as  before. 

10.  It  is  important  to  consider  the  effect  of  equal  roots  on  the 
test  of  stability  given  in  Art.  8  of  this  chapter. 

In  this  case  we  know  that  the  roots  of  the  minor  An  separate 
the  roots  of  A.  If  therefore  A  have  two  negative  equal  roots,  it  is 
clear  that  An  must  have  one  of  them.  In  the  same  way  A22,  A^, 
&c.  must  also  all  vanish  for  this  value  of  m.  Since 


it  follows  as  in  Art.  2  if  A  and  App  both  vanish,  that  Ap?  and  A9P 
also  vanish.  Hence  all  the  first  minors  of  the  determinant  A 
vanish.  This  is  the  case  considered  in  Art.  5  of  Chap.  I.  The 
equal  roots  instead  of  introducing  into  the  expressions  for  the  co- 
ordinates terms  which  contain  t  as  a  factor  merely  render  two  of 
the  coefficients,  instead  of  one,  indeterminate. 

In  the  same  way  if  the  equation  A  =  0  is  satisfied  by  three 
values  of  mz  equal  to  the  same  negative  quantity,  the  equation 
A^  =  0  must  have  two  of  them,  and  its  principal  minor  must  have 
one  of  them.  Reasoning  as  before  we  see  that  all  the  second 
minors  of  A  must  be  zero.  This  is  again  the  test  that  there  should 
be  no  terms  which  contain  t  as  a  factor. 

The  presence  therefore  of  equal  roots  does  not  in  the  theorem 
of  Art.  8  affect  the  stability  of  the  motion. 

When  a  system  is  disturbed  from  a  position  of  equilibrium 
whether  stable  or  unstable,  the  roots  of  the  fundamental  determi- 
nant are  separated  by  its  minors  in  the  manner  described  in 
Art.  8  of  this  Chapter.  By  what  has  just  been  proved,  we  see 
that  if  the  fundamental  determinant  have  equal  roots,  whether 
positive  or  negative,  these  do  not  introduce  into  the  integrals 
terms  which  contain  t  as  a  factor. 
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11.  The  proposition  in  the  last  article  may  be  made  more 
general.  If  the  fundamental  determinant  be  reduced  to  the  form 
indicated  in  Art.  7,  we  shall  show  that  if  A  vanish  for  two  equal 
negative  values  of  m*  which  are  numerically  greater  than  the 
greatest  negative  quantity  in  the  series  an,  a?2,  «33,  &c.,  then 
these  equal  roots  will  not  introduce  any  terms  into  the  solution 
with  t  as  a  factor.  If  «n,  «22,  &c.  are  all  positive,  this  reduces  to 
the  proposition  proved  in  the  last  article. 

Following  the  same  notation  as  before,  we  have 


If  neither  Au  nor  A^  are  zero,  they  must  have  the  same  sign 
when  A  vanishes  for  a  negative  value  of  m2.  For  their  product 
is  equal  to  A12A21  which  has  been  proved  in  Art.  2  to  be  positive. 
Hence  all  the  leading  first  minors,  viz.  Au,  A22,  &c.  must  have  the 
same  sign  for  any  negative  value  of  m2  which  makes  A  vanish. 

By  differentiation  we  have 

^  =  20zAu  +  a12A12  +  .  .  .  -  c*12A21  f  2mA23  +  .  .  .  -  &c. 

But  we  have  also 

A  =  (m*  +  aj  An  +  auw  A18  +  ... 


A  =  -  al2m  A21  +  (m2  +  aj  A22  +  ... 
&c.  =  &c. 

Hence  if  n  be  the  highest  power  of  m  occurring  in  A,  we  have 

^A  -     A  2  A  2 

If  then  A  and  ^—  both  vanish  for  any  negative  value  of  m* 

greater  than  the  greatest  negative  quantity  in  the  series  an,  a22, 
&c.,  we  have  the  sum  of  a  number  of  quantities  all  of  the  same 
sign  equal  to  zero.  This  requires  that  each  should  be  zero.  We 
have  therefore  Au  =  0,  A22  =  0,  &c.  The  rest  of  the  proof  is  the 
same  as  before. 

By  differentiating  the  expression  for  ^—  and  substituting  for 

-Y-"  ,   —j-^ ,  &c.  their  values  in  terms  of  their  first  leading  minors, 

we  may  extend  this  proposition  to  the  case  in  which  the  funda- 
mental determinant  has  three  equal  roots,  and  so  on. 
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12.  A  sphere  is  suspended  *by  a  string  OA  from  a  fixed  point 
O,  and  is  set  rotating  about  a  vertical  diameter  which  is  in  the 
same  straight  line  as  the  string  with  an  angular  velocity  n.  A 
small  disturbance  is  given,  determine  if  the  steady  motion  is 
stable. 

Let  0  be  the  origin,  and  let  the  axis  of  z  be  vertically  down- 
wards, let  Ix,  ly,  I  be  the  co-ordinates  of  A,  the  point  at  which  the 
string  is  attached.  Let  C  be  the  centre,  and  let  af,  arj,  a  be  the 
co-ordinates  of  C  relative  to  A.  Then,  exactly  as  in  Chap.  IV.  Art. 
22,  the  Lagrangian  function  may  be  shown  to  be 


the  mass  being  taken  unity. 

Putting  %  =  n,  f ,  if,  x ',  y  all  zero,  we  see  that  L  is  a  maxi- 
mum when  x,  y,  £,  rj  are  zero,  the  steady  motion  is  therefore  stable 
for  all  values  of  n. 

If  we  put  &2  =  |  a2,  and  m2  =  —  X2,  so  that  x,  y,  &c.  are  all  repre- 
sented by  terms  of  the  form  ^A  cos  (A,£  +  a),  we  may,  by  the 
methods  of  the  last  chapter,  prove 


This  equation,  whatever  may  be  the  sign  of  n,  has  two  positive 
and  two  negative  roots.     All  four  give  stable  oscillations. 


R.A. 


CHAPTEE    VI. 


If  tJie  energy  of  the  system  be  a  maximum  or  minimum  under  certain 
conditions,  the  motion  whether  steady  or  not  is  stable.  Arts.  1  —  3. 

When  the  motion  is  steady,  it  will  be  also  stable  if  a  certain  function  of 
the  co-ordinates  called  V  +  W  is  a  minimum.  Art.  4. 

If  there  be  only  one  co-ordinate  which  enters  into  the  Lagrangian  function, 
except  as  a  differential  coefficient,  this  condition  is  necessary  and 
sufficient.  Arts.  5,  6. 

Additional  conditions  when  there  are  two  co-ordinates.     Art.  8. 

1.  When  a  system  is  oscillating  about  a  position  of  equilibrium, 
it  is  well  known  that  we  may  determine  the  stability  or  insta- 
bility of  the  equilibrium  by  what  we  may  call  the  "  energy  crite- 
rion." This  criterion  may  also  be  sometimes  used  when  the 
system  is  oscillating  about  a  state  of  steady  motion. 

Let  E  be  the  sum  of  the  kinetic  and  potential  energies  of  the 
system.  Then  throughout  any  motion  of  the  system  we  have 

E=li, 

where  h  is  a  constant  depending  on  the  initial  conditions.  If  6,  <£, 
&c.  are  the  co-ordinates  of  the  system,  E  is  a  known  function  of 
6,  &  ',  $,  <£',  &c.  Suppose  that  some  of  the  other  first  integrals  of  the 
equations  of  motion  are  known.  Let  these  be 


J£  (4  ?,*+',&&)  -0.  1 

&c.  ) 

the  time  t  being  absent. 

For  the  purposes  of  this  proposition  let  us  suppose  0,  ff,  $,  </>', 
&c.  to  be  separate  variables  unconnected  with  each  other  except 
by  the  equations  just  written  down. 

If  E  be  an  absolute  maximum  or  an  absolute  minimum  for 
all  variations  of  6,  0',  &c.,  those  corresponding  to  the  given  motion 
making  E  constant,  then  that  motion  is  stable  for  all  displacements 
which  do  not  alter  the  constants  Clf  (\,  &c. 
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If  this  proposition  be  not  evident,  it  may  be  proved  by  elimi- 
nating as  many  of  the  letters  as  possible.  If  0,  #',  &c.  be  the  re- 
maining co-ordinates  we  have 


Let  h  be  the  value  of  E  in  the  given  motion,  and  let  the 
system  be  started  in  some  slightly  different  manner  so  that 


If  E  be  a  maximum  along  the  given  motion,  then  any  change 
whatever  in  0,  6\  &c.  decreases  E.  Hence  0,  0',  &c.  cannot  deviate 
so  much  from  their  values  along  tlie  given  motion  that  the  change 
in  E  becomes  greater  than  §h. 

2.  Let  us  apply  this  principle  to  a  system  of  bodies  which 
moves  in  steady  motion  with  some  co-ordinates  0,  $,  &c.  such  that 
their  differential  coefficients  6',  <£',  &c.  are  constant,  and  the  re- 
maining co-ordinates  f,  77,  &c.  themselves  constant.  Let  us  further 
suppose  that  the  energy  is  a  function  of  6',  </>',  &c.,  but  not  of 
6,  (f).  By  Lagrange's  Equations  we  have  the  integrals 


-C     --C    &c 
dff~  °"    df~°2)i 

It  is  clear  that  the  system  can  describe  any  one  of  a  number 
of  steady  motions,  which  we  have  already  called  parallel  motions, 
and  which  are  determined  by 

0'  =  p,         <f>'  =  q,  &c. 
f  =  or,  77  =  /3,  &c. 

where  p,  q,  &c.  a,  /3,  &c.  are  constants  which  satisfy  all  Lagrange's 
equations. 

We  have  thus  as  many  relations  between  these  constants  as 
there  are  co-ordinates  f,  rj,  &c.  Let  the  system  be  started  with 
any  initial  conditions  we  please,  then  the  constants  Clf  (72,  ...  are 
given.  These  being  known  we  have  as  many  relations  between 
the  constants  of  steady  motion  as  there  are  co-ordinates  6,  </>,  &c. 
The  steady  motion  is  therefore  determined.  If  the  energy  of  this 
initial  motion  is  nearly  equal  to  that  of  this  steady  motion,  and 
if  it  be  a  maximum  or  minimum  as  explained  above,  then  the 
system  will  never  deviate  far  from  its  corresponding  position  in 
the  steady  motion,  and  this  steady  motion  may  be  called  stable. 

3.  Example.  A  top  is  set  spinning  on  its  point  on  a  perfectly 
rough  horizontal  ground,  with  its  axis  inclined  to  the  vertical,  find 
the  condition  of  stability. 

Let  6  be  the  inclination  of  the  axis  00  of  the  top  to  the 
vertical  OZ;  -^  the  inclination  of  the  plane  ZOC  to  a  vertical 

6—2 
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plane  fixed  in  space,  and  <£  the  inclination  to  a  plane  through  OG 
fixed  in  the  body.  Let  0  be  the  apex,  G  the  centre  of  gravity 
which  lies  in  OC,h=OGr.  Let  A,  A,  C  be  the  principal  mo- 
ments of  inertia  at  0,  and  M  the  mass  of  the  top.  We  have  then 


cos 


+  Mgh  cos  0. 


By  La-grange's  equations  we  have  the  integrals 

$  +  i|r'  cos  0  =  nt 
Cn  cos  6+  A  sin20tjr'  =  m, 

where  n  and  m  are  two  constants,  the  former  representing  the 
angular  velocity  of  the  top  about  its  axis,  and  the  latter  the 
angular  momentum  about  the  vertical.  If  we  now  eliminate 
<£'  and  i|r'  we  find  that  E  is  a  minimum  when  6  =  a,  if 

(7V  >  kMghA  cos  a, 
which  is  the  result  given  by  other  methods. 

4.     The  theorem  of  Art.  2  may  be  put  into  another  form. 
Let  the  kinetic  energy  be 


Then  since  0,  <£,  &c.  are  absent  from  the  coefficients,  we  have 
the  integrals 

-  T^f  -  &c. 

-  Ttff  -  &c. 
&c.  =  &c. 

For  the  sake  of  brevity  let  us  call  the  right-hand  sides  of  these 
equations  Q  —  X,  C2  —  Y,  &c.  Since  T  is  a  homogeneous  func- 
tion of  6',  0',  &c.,  we  have,  as  in  Chap.  iv.  Art.  21, 


If  we  substitute  in  the  second  line  the  values  of  6',  <j>,  &c. 
found  by  solving  the  linear  equations  just  written  down,  we  have 
the  determinant 

0,       Q  +  X,   CZ+Y,  &c. 
-X       T  T         &c. 


,-*; 

&c. 


&c. 


&c.       &c. 


2A 
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where  A  is  the  discriminant  of  the  terms  in  T  which  contain  0', 
&c.     This  determinant  is  unaltered  by  changing  the  signs  of  X, 
&c.  and  is  a  quadratic  function  of  Clt  <72,  &c.,  X,  Y,  &c.     Hence 
the  terms  C^Y,  CZY,  &c.  do  not  occur.     If  then  we  put 


0  C19  Cz)  &c. 

Clt  Te&,  T9h  &c. 

<72,  2%,  TH,  &c. 

&c.  &c.  &c.  &c. 

we  have 


where  B#,  £c.  are  independent  of  Clt  Cz,  &c.  Now  Tis  essentially 
positive  for  all  values  of  the  variables,  and  therefore  for  such  as 
make  Glt  (72,  &c.  all  zero.  Hence  the  terms  involving  f ',  77',  &c. 
are  together  a  minimum  when  f ,  77',  &c.  are  all  zero.  The  coeffi- 
cients Bgt  &c.  may  all  be  treated  as  constants  since  f ',  rj',  &c.  are 
all  small  quantities. 

If  V  be  the  potential  energy,  we  have  therefore  the  following 
rule.  If  W  +  Y  be  a  minimum  for  all  variations  of  %,  77,  &c.  then 
the  steady  motion  is  certainly  stable.  It  should  be  noticed  that 
W+  V  is  a  function  only  of  f,  97,  &c.  the  corordinates  which  are 
constant  in  the  steady  motion. 

5.  If  the- energy  be  a  function  of  one  only  of  the  co-ordinates, 
tJiough  the  differential  coefficients  of  all  the  others  enter  into  its 
value,  this  condition  is  sufficient  and  necessary. 

Let  f  be  this  co-ordinate.     Then  by  vis  viva  we  have 

D  r 


Differentiating  we  have 


This  equation  must  be  satisfied  by  the  steady  motion  repre- 
sented  by   f  =  cc.     The   second   term      ^    7j" — -'  must  therefore 

at; 

vanish  when  {?  =  «,  so  that    W+  Fis  a  maximum  or  minimum. 
To  fiud  the  oscillation  let  us  put  f  =  a  +  x,  we  find 
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where  the  square  bracket  implies  that  f  is  to  be  put  equal  to  a 
after  differentiation.  By  the  same  reasoning  as  before  B^  is 
necessarily  positive,  and  the  motion  will  be  stable  or  unstable  ac- 
cording as  (W  +.  V)  is  a  minimum  or  maximum. 

6.  If  we  refer  to  Art.  21  of  Chap.  IV.  we  see  that  this  function 
.W  +  Vis  the  value  of  the  modified  Lagrangian  function  U  when 
f,  77',  &c.  are  all  put  zero,  and  the  sign  of  the  whole  function 
changed.  It  therefore  follows  by  Chap.  V.  that  when  W+  V  is  a 
minimum  the  steady  motion  is  stable.  The  "energy  criterion  of 
stability,"  as  far  as  it  applies  to  steady  motion,  may  therefore  be 
deduced  from  that  given  in  Chap.  v.  Art.  8.  But  the  mode  of 
demonstration  adopted  in  that  chapter  gives  us  more  information 
as  the  nature  of  the  motion,  while  the  modes  of  application  to 
examples  of  the  two  criteria  are  quite  different.  The  energy 
criterion  may  also  be  sometimes  applied  to  determine  the  stability 
of  a  motion  which  is  not  steady. 

[The  relation  between  the  theorem  in  this  Chapter  in  which 
E—  T+  V  is  made  a  minimum  to  that  given  in  Chapters  IV.  and  V. 
in  which  L=  T—  V  is  made  a  maximum  may  be  more  distinctly 
perceived  by  the  following  statement. 

Let  x,  y,  &c.  be  the  co-ordinates  of  the  system,  and  let  L  be 
the  Lagrangian  function*  so  that  L  =•  T—  V>  then  by  Art.  3  of 
Chap.  IV.  the  co-ordinates  in  steady  motion  satisfy  the  equations 

dL     A        dL          p 

-y-=0,       -7-  =  0,  &c  ......................  (1). 

dx  dy 

Here  L  is  expressed  as  a  function  of  x,  y,  &c.  x,  y\  &c. 

Suppose  some  of  the  co-ordinates  as  0,  <£,  &c.  are  absent  from 
the  expression  for  L,  so  that  L  is  a  function  of  &  ',  (/>',  &c.,  the  re- 
maining co-ordinates,  viz.  f  ,  77,  &c.  and  their  differential  coefficients. 
Then  if  we  form  the  modified  Lagrangian  function  as  in  Art.  21 
of  Chap.  IV.  the  equations  (1)  of  steady  motion  become 


=  0,  =  0,&o. 


Here,  as  in  the  Hamiltonian  equations  (7,,  <72,  &c.  are  the 
6,  </>,  &c.  components  of  momentum,  and  L'  is  expressed  as  a 
function  of  (7t  ,  678,  &c.  f,  77,  &c.  f,  77',  &c. 

jrn 

*  Let  M,  v,  &c.  be  the  a;,  y,  &c.  components  of  momentum,  so  that  u  =  —lt  &c. 

and  let  H  be  the  Hamiltonian  function.     Then  H=T+  V  and  we  easily  deduce 
from  the  Hamiltonian  equations  that,  in  steady  motion, 

dH  dH 

-y~  =  0,  -   =0,  &C. 

dx  dy 

Here  H  is  expressed  as  a  function  of  x,  y,  &c.  t/,  r,  £c. 
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But  L  =  r-  C/  -  Crf  -  &C.  -  V 


by  Euler's  theorem  on  homogeneous  functions.     In  steady  motion 
f,  T/>  &c-  aU  vanish,  hence  the  equations  (2)  become 


_d(W+V)_  d(W+V)_ 

~~    0>      ~        • 


where   TF  is  the  value  of  T  when  f,  ?;',  &c.  are  all  put  equal  to 
zero.     Here  JFis  expressed  as  a  function  of  019  Oa,  &c.  f,  T;,  &c. 

It  is  shown  in  Chap.  V.  that  if  the  Lagrangian  function  ex- 
pressed as  required  in  equations  (1)  be  a  maximum  the  motion 
is  stable.  It  is  shown  in  this  Chapter  that  if  the  function  TF-f  V 
be  a  minimum  the  .motion  is  stable.] 

7.  To  find  the  condition  of  stability  when  the  Lagrangian 
function  is  a  function  of  two  only  of  the  co-ordinates,  though  the 
differential  coefficients  of  all  the  others  enter  into  its  value. 

Let  f,  rj  be  these  two  co-ordinates,  then  the  modified  Lagran- 
gian function  as  explained  in  Art.  20  of  Chap.  IV.  will  be  a  func- 
tion of  f,  ?;,  f,  r;'  only. 

Let  the  steady  motion  be  given  by  f  =  a,  vj  =  /:?,  with  the 
corresponding  values  of  the  other  co-ordinates  6,  <f>,  &c.  Then  a 
and  /3  are  constants.  Let  £  =  a  +  x,  ?7  =  ^-f-y,  and  let  us  expand 
the  modified  Lagrangian  function  in  powers  of  x,  y.  Neglecting 
the  terms  of  the  first  order,  as  they  only  give  the  steady  motion, 
let 


+  Cnxx+C^y'  + 

Also  let  E—  (721—  C7U.     Then  the  condition  of  stability  is  that 
the  roots  of  the  following  equation  should  be  of  the  form  ft  J—  1. 

=  0. 

I  B^nf  —  An  —  Em, 

If  A'  =  I 


and 
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be  the  two  discriminants  and  the  other  invariant,  this  leads  at 
once  to  the  conclusion  that  the  motion  is  stable  only  when 

(1)  A  is  positive, 

(2)  Ez-  0  is  positive  and  >  2 


If  AU^  +  A&y  +  A2  2  is  a  maximum  when  as  and  y  are  zero, 
the  two  conditions  are  obviously  satisfied. 

This  condition  may  be  otherwise  expressed  ;  if  L'  be  the  modi- 
fied Lagrangian  function,  a  steady  motion  is  given  by 

f  =°'  f  -o.  r-ft  v=o. 

This  motion  will  be  stable  if  for  the  values  of  £,  y  thus  found, 


d*L'd*L'  . 

(1)     -j-&  -j-r  —  JTT-  is  positive. 
d%     drf 


'  d*L' 
(d?  d^        drj2   d£*  ~ 


is  positive  and  greater  than 


?    dri' 

8.  The  nature  of  the  motion  when  thus  reduced  to  depend  on 
two  co-ordinates  may  be  illustrated  by  geometrical  reasoning.  Let 
the  position  be  defined  by  two  co-ordinates  x,  y  which  are  zero 
along  the  steady  motion.  Let  these  be  regarded  as  the  co- 
ordinates of  a  point  P  referred  to  any  axes.  Then  the  motion  of 
P  exactly  represents  that  of  the  system. 

Let  us  construct  the  conies 

~,2  ,,.2 


Then,  exactly  as  in  Arts.  11  —  14  of  Chap,  iv.,  it  will  be  found 
convenient  to  transform  the  co-ordinates  by  writing 


If  IJL  be  the  modulus  of  transformation  we  have  n  =  aj>2  —  aj\. 
It  is  easy  to  see  by  actual  substitution,  if  E'  be  the  difference  of 
the  coefficients  of  pq  and  pq,  that 
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If  the  transformation  be  from  one  set  of  oblique  co-ordinates 
to  another,  let  &>,  to  be  the  angles  between  the  axes.  We  then 
have 

E'         E 

sin  ft/     sin  &>  ' 

Transforming  the  axes  to  the  common  conjugate  diameters,  the 
conies  become 


the  signs  of  a  and  b  being  such  as  to  make  these  conies  real.     The 
equation  to  find  m  becomes 


=  0. 

It  is  therefore  necessary  for  stability  that  the  conic  a  should  be  an 
ellipse  as  well  as  the  conic  b.     It  is  also  necessary  that 

Ef 


both  roots  having  the  same  sign  and  the  inequality  being  nu- 
merical. 

Let  OP,  OP  ;  0  Q,  0  Q'  be  the  common  conjugates  of  the  two 
conies,  this  condition  then  becomes 

E      area  of  conic  b      OQ      OQ 
sin  a)'"       Tr^ab          >  OP+ OP" 

When  the  system  describes  an  oscillation  with  one  period,  i.e. 
an  harmonic  oscillation,  the  path  of  the  representative  particle  is 
easily  seen  to  be 

(5U'W  -  A")/  +  0&25>2  -  4B")  <f  =  constant. 

The  harmonic  paths  are  therefore  ellipses.  It  also  appears  that 
the  two  ellipses  which  represent  the  two  harmonic  vibrations  and 
the  two  ellipses  a  and  6  have,  all  four,  a  common  set  of  conjugate 
diameters. 
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Any  small  term  of  a  high  order,  if  its  period  is  nearly  the  same  as  that  of 
an  oscillation  of  the  system,  may  produce  important  effects  on  the 
magnitude  of  the  oscillation.  Art.  1. 

Origin  of  such  terms,  with  an  example.     Arts.  2 — 3. 

Supposing  the  roots  of  the  determinantal  equation  to  satisfy  the  conditions 
of  stability  to  a  first  approximation,  yet  if  a  commenstirable  relation 
hold  between  these  roots  it  is  necessary  to  examine  certain  terms  of  the 
higher  orders  to  determine  whether  they  will  ultimately  destroy  tlie 
stability  of  the  system.  Art.  4. 

If  a  certain  relation  hold  among  the  coefficients  of  these  terms,  they  will 
not  affect  the  stability  of  the  system,  but  only  slightly  alter  the  periods 
of  oscillation.  Arts.  5 — 7. 

Examples,  the  first  taken  from  Lagrange's  method  of  finding  the  oscilla- 
tions about  a  position  of  equilibrium.  Arts.  8 — 9. 

If  the  coefficients  of  the  equations  of  motion  should  not  be  strictly  con- 
stant, but  only  nearly  so,  the  stability  will  not  be  affected,  unless  the 
reciprocals  of  their  periods  have  commensurable  relations  with  the 
reciprocals  of  the  periods  of  oscillation  of  the  system.  Art.  10. 

1.  If  we  understand  that  a  motion  is  called  stable  when  any 
small  disturbance  does  not  cause  the  system  to  deviate  far  from  its 
undisturbed  motion,  it  is  clear  that  we  cannot  be  certain  of  the 
stability  without  examining  the  terms  of  the  second  order.  It  is 
possible  that  some  of  these  may  have  their  periods  so  timed  that 
their  effects  may  accumulate  until  the  motion  is  changed. 

Returning  to  the  equations  of  Art.  3,  Chap.  I.  we  shall  have  on 
the  right  hand,  instead  of  zero,  a  series  of  small  terms  of  orders 
higher  than  the  first.  To  find  a  second  approximation,  we  substi- 
tute the  values  of  x,  y,  &c.  given  at  the  end  of  Art.  3,  in  these 
terms. 

They  will  therefore  take  the  form  Neni  and  will  produce  in 

Nf 

y,  y,  &c.  terms  of  the  form  ^-r-r  ent,  where  N'  is  of  the  same  order 

at  least  as  the  term  considered,  and/(w)  has  the  same  meaning  as 
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in  Chap.  I.  These  will  have  to  be  expressed  in  trigonometrical 
real  forms,  but  it  is  unnecessary  to  exhibit  the  process,  for  we  see 
at  once  that  no  small  term  or  force  (whatever  it  may  be  called)  of 
a  high  order  can  affect  the  stability  of  the  motion  unless  it  makes 
f(n)  very  nearly  or  exactly  equal  to  zero.  In  this  case  its  period 
is  very  nearly  or  exactly  equal  to  one  of  the  periods  of  the  motion 
given  by  taking  terms  of  the  first  order  only. 

A  remarkable  use  of  this  principle  was  made  by  Captain  Kater 
in  his  experiments  on  the  magnitude  of  gravity.  It  was  important 
to  determine  if  the  support  of  his  pendulum  was  perfectly  firm. 
He  tells  us  that  he  had  recourse  to  a  delicate  and  simple  instru- 
ment the  sensibility  of  which  was  so  great  that  had  the  slightest 
motion  taken  place  in  the  support  it  must  have  been  instantly 
detected.  The  instrument  consists  of  a  steel  wire  the  lower  part 
of  which,  inserted  in  the  piece  of  brass  which  serves  as  its  support, 
is  flattened  so  as  to  form  a  delicate  spring.  On  the  wire  a  small 
weight  slides  by  means  of  which  it  may  be  made  to  vibrate  in  the 
same  time  as  the  pendulum  to  which  it  is  to  be  applied  as  a  test. 
When  thus  adjusted  it  is  placed  in  the  material  to  which  the 
pendulum  is  attached,  and  should  this  not  be  perfectly  firm  its 
motion  will  be  communicated  to  the  wire,  which  in  a  little  time 
will  accompany  the  pendulum  in  its  vibrations.  This  ingenious 
contrivance  appeared  fully  adequate  to  the  purpose  for  which  it 
was  employed,  and  afforded  a  satisfactory  proof  of  the  stability  of 
the  point  of  suspension.  See  Phil.  Trans.,  1818. 

2.  Since  the  term  Nent  is  obtained  by  compounding  the 
different  terms  in  the  values  of  x,  y}  &c.  it  is  clear  that 


where  p,  q,  £c.  are  positive  integers  whose  sum  is  the  order  of  the 
term.  It  is  therefore  only  when  the  roots  of  the  dynamical  equa- 
tion f  (ni)  =  0  are  such  that  a  linear  relation  of  the  form 

pm1  +  qm^  -f  ...  =  mi  very  nearly 

exists  between  them,  that  we  may  expect  to  find  important 
terms  among  the  higher  orders.  The  order  of  the  terms  to  be 
examined  will  be  p  +  q  +  .  .  .,  and  unless  this  be  also  small,  the  terms 
will  probably  remain  insignificant.  If  the  root  ml  should  occur 
twice  in  /(m)  =0  it  is  clear  that  the  divisor  /  (ri)  will  be  a  small 
quantity  of  the  second  order,  and  the  term  may  be  said  (as  in 
the  Lunar  Theory)  to  rise  two  orders. 

3.     To  take  an  example,  let  us  suppose  a  particle  to  be  describ- 
ing an  ellipse  about  a  fixed  centre  of  force  in  one  focus.     If  dis- 
turbed it  will  describe  a  slightly  different  ellipse.     Since 
r  =  a{l  —  ecos  (nt  4-  e  —  o>}+  ...}, 
0  =-nt+  e  +  2e  sin  (nt  +  e-  o>)  -f  ... 
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we  see  that  a  slight  change  in  the  elements  will  cause  variations 

in  r  and  6  of  the  period  —  ,  an  additional  variation  in  6  of  the 

n 

form  t&n  4-  Be  and  an  additional  variation  in  r  of  the  form  $a.  All 
these  variations  should  by  Art.  3  of  Chap.  i.  be  indicated  by 
expressions  of  the  form 


where  the  values  of  m  are  the  roots  of  the  equation  f(m]  =  0. 
The  roots  therefore  of  the  equation  /(w)=0  for  $r  are  m  =  0 
and  ±  n*J  —  1,  and  for  &d  are  m  =  0,  0  and  ±  n  V  —  1.  We  there- 
fore infer  that  any  small  disturbing  causes  of  the  second  order 
whose  periods  are  nearly  equal  to  that  of  the  particle,  will  cause 
important  inequalities  in  both  $r  and  BO,  and  (since  /  (m)  =  0  has 
two  roots  equal  to  zero)  any  term  of  long  period  will  rise  two 
orders  in  S0. 

4.     If  the  roots  of  the  subsidiary  equation  are  such  that  the 
relation 


holds  accurately,  the  solution  changes  its  character.  We  have 
now  in  the  value  of  a?  a  term  of  the  form  temit.  Unless  the  real 
part  of  mt  is  negative,  this  indicates  that  the  system  will  depart 
widely  from  the  motion  which  we  took  as  a  first  approximation. 
We  must  therefore  modify  our  first  approximation  (as  in  the  Lunar 
Theory)  by  including  in  it  the  terms  which  produced  these  im- 
portant effects.  We  may  then  enquire  how  far  this  modified  first 
approximation  indicates  that  the  motion  is  stable  or  unstable. 
When  these  terms  are  included  the  equations  to  be  solved  are  in 
general  no  longer  linear,  and  it  is  sometimes  impossible  to  find  a 
solution  sufficiently  accurate  to  serve  as  a  first  approximation 
throughout  the  whole  motion. 

5.     In    some    cases,   however,   the    oscillations   may   still   be 
represented  by  expressions  of  the  form 


&c., 

where  the  values  of  n^  nz...  differ  but  slightly  from  the  roots  of 
the  equation  f(m)  =  0.  Let  us  investigate  the  condition  that  this 
should  be  true,  and  also  determine  whether  the  changes  in  the 
values  of  mv  mz...  are  sufficient  to  affect  the  stability. 

Suppose  that  we  have  completed  our  first  approximation,  and 
find  on  proceeding  to  a  higher  approximation  that  the  terms 
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present  themselves  on  the  right-hand  side  of  the  first  set  of 
equations  in  Art.  3,  Chap.  i.  These  terms  are  supposed  to  have 
arisen  from  several  relations  of  the  form 


If  these  terms  can  be  included  in  the  first  approximation  by 
writing  nv  nz,  £c.  for  mlt  m^  &c.  we  have,  by  substitution  in  the 
differential  equations,  certain  equations  connecting  n,  M,  M',  &c., 
whose  left-hand  sides  are  the  same  as  those  used  in  Art.  3  with 
«!  written  for  m^  but  on  the  right-hand  sides  we  have  instead  of 
zero  the  quantities  Nv  NJ,  &c.  The  test  of  the  success  of  the 
process  is  that  these  modifications  in  the  values  of  ra  must  satisfy 
the  same  relations  as  before. 

Now  Nlt  Nl,  &c.  are  all  at  least  of  the  second  order  of  small 
quantities,  hence  up  to  terms  of  the  first  order  the  ratios  Mv  M^, 
&c.  will  be  the  same  as  before,  so  that  we  may  put 

M1  =  L1alt  M^Lfa  &c., 
following  the  same  notation  as  in  Art.  3,  Chap.  I.     We  also  have 


Let  n^  =  m^  +  8mx  we  find 

cm,  = 
Similarly 


&c.  =  &c. 

It  is  evident,  by  the  theorem  of  determinants  alluded  to  in 
Art.  3,  that  these  are  symmetrical  expressions. 

6,     We  may  conveniently  express  these  results  in  the  form  of  a 
rule. 

Suppose  we  have  to  a  first  approximation 


Eliminate  from  the  differential  equations  all  the  variables  except 
x  in  the  usual  manner.  This  may  be  done  by  performing  on  the 
several  equations  the  operations  represented  by  the  minors  a,  a, 

&c.,  -T-  being  written  for  m.     Let  the  equation  thus  found  be 
ctt 
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Then  all  these  terms  can  be  included  in  the  first  approximation, 
provided 

P  P 

* 


satisfy  the  relations 

p^m^  +  q$m2  +  .  ..  =  Smi, 

&c.  =  &c. 
which  exist  among  the  roots  of  the  dynamical  equation. 

7.  The  general  results  we  have  arrived  at  may  be  summed 
up  as  follows.     Though  some  of  the  terms  of  the  higher  orders 
may  affect  the  magnitude  of  the  oscillation,  yet  no  term  will  arise 
to  affect  the  stability  of  the  motion  unless  there  be  some  relations 
between  the  roots  of  the  dynamical  equation  of  the  form 

pm^  +  qm2  +  .  .  .  =  ml  , 

where  p,  q,  &c.  are  all  integers.  Even  if  such  relations  occur,  the 
lowest  order  of  the  term  is  p  +  q  +  ...»  and  if  this  be  considerable 
the  term  will  not  produce  any  important  effects  until  a  con- 
siderable time  has  elapsed.  If  a  certain  relation,  just  found,  hold 
among  these  terms,  their  only  effect  is  slightly  to  modify  the 
periods  of  oscillation,  without  altering  the  type  of  motion. 

8.  As  an  example,  let  us  consider  a  system  of  bodies  to  be 
oscillating  about  a  position  of  equilibrium.    We  know  by  Lagrange's 

feneral  solution,  that  the  equation  f(ni)  =  0  is  of  an  even  order. 
ts  roots  are  of  the  form 


—  aV—  1,     ?H3  =  /3     —  1,  &c. 

Whatever  the  numerical  values  of  these  may  be  we  have 
m1  +  m3  +  m4  =  mlt        m2  +  ms  +  m4  =  m2  ,  &c. 

so  that  the  small  terms  of  the  third,  fifth,  &c.  orders  might  affect 
the  stability  of  the  oscillation.  But  we  shall  now  show  that  they 
only  affect  the  periods  of  oscillation,  and  not  the  stability  of  the 
system. 

Since  both  sides  of  Lagrange's  equations  must  be  of  —  2  di- 
mensions in  time  and  the  impressed  forces  are  also  of  —  2  dimen- 
sions, it  is  clear  that  these  terms  must  consist  of  powers  of  x,  y,  &c. 

^>    Tjjf  >   &c-   and   products  of  an   even   number   of  factors  of 


-j-  ,    -~  ,   &c.     We  know  also,  by  Lagrange's  solution,  that  the 
co-ordinates  take  the  form 

ic  —  Mt  cos  (it  +  \)  +  &c» 

y  =  M[  cos  (it  +  \)  +  &c. 
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The  minors  a,  a,  &c.  are  also  all  even  functions  of  y  ,   hence 
the  equation  found  after  elimination  is  of  the  form 

X  =  P  cos  (at  +  X)  +  Qcos  (J3t+fi)  +  ... 


Replacing  aV—  1,   —  aV  —  1,  &c.  by  mlt  mz,  &c.  we  find  by 
Art.  6, 

,  P  ,  P 


Since  /"(»H)  is  of  odd  dimensions,  and  ml=—m9  we  clearly  have 
3??^=—  8w2,  and  therefore  the  test  is  satisfied. 

9.  [As  another  example  let  us  apply  the  rule  of  Art.  6  to 
some  very  simple  case  which  will  involve  no  algebraical  substi- 
tutions of  any  length. 

The  motion  of  a  simple  pendulum  under  the  action  of  gravity 
may  be  made  to  depend  on  the  equation 


(1), 


where  a  and  ft  are  two  constants  and  x  is  the  inclination  of  the 
pendulum  to  the  vertical  which  is  supposed  to  be  small.  The  first 
approximation  to  the  motion  is 

x  =  M,?™*  +  MjP*  .....................  ..(2), 


where  m^  and  ra2  are  the  roots  of  the  equation  m2  H-  a2  =  0.  Our 
object  is  to  ascertain  by  help  of  the  rule  given  in  Art.  6  whether 
the  small  force  represented  by  Boi?  renders  this  first  approximation 
unstable  or  merely  slightly  alters  the  numerical  values  of  ml  and 
mz. 

The  two  roots  are  connected  by  the  relation. 

?^  +  ^2  =  0  .............................  (3). 

Substituting  the  value  of  x  on  the  right-hand  side,  we  have 


Hence  by  the  rule  in  Art.  6 

g        SfIMM 
W'~   2m,Mt  '          m'- 

These  clearly  satisfy  the  relation 

Sm,  +  Bm,  -  0, 
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and  therefore  the  first  approximation  taken  above,  so  far  as  the 
disturbing  force  Bx*  is  concerned,  is  stable. 


If  the  small  force  had  been  ft   -        instead  of  fix3,  it  is  easy 
to  see  in  the  same  way  that 


.  gm  _ 

™~ 


so  that  the  relation  Sm1  +  Sw^  =  0  would  not  have  been  satisfied. 
The  first  approximation  taken  above  is  therefore  not  sufficiently 
accurate  to  serve  as  a  first  approximation  throughout  the  motion. 

In  this  example  we  have  considered  the  effect  of  a  small  force 
of  the  third  order  in  disturbing  the  stability  of  the  motion  given 
by  equation  (1).  The  same  equation  will  obviously  occur  in  many 
other  cases  of  motion.  For  example,  let  a  particle  describe  a  cir- 
cular orbit  about  a  centre  of  force  situated  in  the  centre.  If 
slightly  disturbed  the  equation  giving  the  disturbance  x  in  the 
radius  vector  takes  the  form 


where  a,  @  and  7  are  constants.     Similar  remarks  will  therefore 
apply  to  this  case  also.] 

10.  When  the  coefficients  of  the  equation  of  motion  are  not 
strictly  constant,  but  yet  do  not  vary  much,  then  we  may  transpose 
the  small  variable  parts  of  these  terms  to  the  right-hand  side  of 
the  eq-uation,  and  treat  their  products  by  the  differential  coeffi- 
cients of  the  co-ordinates  as  small  quantities  of  the  second  order. 
Suppose  the  variable  part  of  one  of  these  coefficients  to  be  p  sin  nt, 
where  p  is  small,  and  let/(ra)  =  0  be  the  equation  giving  the 
periods  of  oscillation  of  the  system  when  the  coefficients  are  taken 
constant.     Then  it  is  clear  that  unless  n  is  nearly  equal  to  the 
sum  or  difference  of  two  values  of  m,  this  term  cannot  rise  into 
importance.     On  proceeding  to  higher  orders  we  see  that  these 
terms  cannot  produce   important   effects   unless  some   commen- 
surable relation  between  m  and  the  roots  of  the  equation  /(w)  =  0 
should  be  very  nearly  satisfied. 

11.  It  should  be  remarked  that  when  the  coefficients  are  not 
constant  it  is  not  a  sufficient  test  of  stability  that  they  should 
always  satisfy  the  conditions  of  stability  obtained  by  giving  them 
their  instantaneous   constant  values.     Thus   if  the   equation  of 
motion  were 

J2*      1 

-T-O  +  Ti2#  =  0, 
2 
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the  coefficient  of  x  is  always  positive,  yet  as  the  equation  is  satis- 
fied by  x  =  a  *Jt,  x  may  become  as  great  as  we  please. 

Even  if  the  coefficients  are  nearly  constant,  we  must  yet  ex- 
amine, by  the  rules  just  given,  if  their  small  changes  are  so  timed 
as  gradually  to  increase  the  oscillation  until  the  divergence  from 
the  given  motion  is  no  longer  small. 

[Suppose  a  system  to  be  oscillating  so  that  its  motion  is  de- 
termined by  the  equation 


where  q  is  a  known  function  of  t,  which  during  the  time  under 
consideration  always  lies  between  /32  and  /3"2  the  latter  being  the 
greater.  Let  the  system  be  started  with  an  initial  co-ordinate  XQ 
and  an  initial  velocity  #0'  in  a  direction  tending  to  increase  x. 
It  may  be  shown  that  the  system  will  begin  to  return,  i.e.  x  will 

/~       #'* 

begin  to  decrease  before  x  becomes  as  great  as  *  /  #02  +  -j*  •     If 

+  m,  +m   be  two  successive  maximum  values  of  x}  we  may  also 

S' 
show  that  m  cannot  be  so  great  as  -^  m,  and  that  the  time  from 

one  maximum  to  the  next  lies  between  -5  and  5-,]. 

M  /j 


R,  A. 
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The  Hamiltonian  Characteristic  or  Principal  functions  when  found  deter- 
mine at  once  the  motion  of  the  system  from  one  given  position  to 
another,  and  whether  the  motion  is  stable  or  unstable.  Arts.  1  —  3. 

Examples  with  a  mode  of  effecting  the  integration  S  =  \Ldt  in  small  oscil- 
lations. Arts.  4  —  7. 

Tlie  Characteristic  function  supplies  the  condition  that  tJie  motion  is  stable 
as  to  space  only,  while  the  Principal  function  gives  the  conditions  that 
it  is  stable  both  as  to  space  and  time.  Art.  8, 

In  what  sense  the  motion  is  unstable  if  either  of  the  two  Hamiltonian 
functions  is  a  minimum.  Arts.  9  —  14. 

1.  If  we  had  any  convenient  methods  of  finding  the  Hamil- 
tonian Characteristic  or  Principal  function,  we  might  determine 
without  difficulty  the  conditions  of  stability  of  a  dynamical  system 
at  the  same  time  that  we  deduce  the  integrals  of  the  equations  of 
motion.     But  it  is  very  difficult  to  discover  either  of  these  func- 
tions by  an  ct,  priori  method.     We  have  indeed  differential  equa- 
tions which  they  must  satisfy,  and  Jacobi  has  taught  us  what  kind 
of  solution  will  serve  our  purpose.     But  the  difficulty  of  finding 
these  solutions  is  as  great  as  that  of  solving  the  equations  of 
motion.     For  these  reasons  it  does  not  seem  necessary  to  dwell  on 
the  uses  of  these  functions. 

2.  Suppose  the  Principal  function  8  of  a  dynamical  system  to 
have  been  found  in  terms  of  the  initial  co-ordinates  $0,  <£0  and  the 
co-ordinates  0}  <£  and  the  time  t.     Let  the  semi-vis  viva  be  given 


where  P,  Q}  R  are  known  functions  of  0,  (/>.     Let  P0,  QQ)  R0  be 
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the  values  of  these  when  #0,  <£0  .are  written  for  0,  (f>.     The  iinal 
integrals  of  the  equations  of  motion  are  then  given  by 


Let  the  system  receive  any  disturbance  at  the  time  t  =  0,  so 
that  while  starting  from  the  same  initial  position,  its  initial  ve- 
locities are  slightly  altered.  Let  x0',  y0f  be  these  initial  changes 
of  00'  and  c/>0'  and  let  6  +  x,  (f>+y  be  the  co-ordinates  of  the 
system  at  the  time  t.  Then  we  have 

d*8  dzS 


d*S 

WMy 

Here  XQ',  yQ'  are  small  arbitrary  quantities,  hence  x  and  y  will 

be  small  if  none  of  the  ratios  of  -£L  ,  J^,   J^_    or 

deQd6     dVjfy     dfadd 
to  the  determinant 


be  large. 

If  the  initial  position  as  well  as  the  initial  motion  be  altered 
we  may  find,  by  a  precisely  similar  process,  the  conditions  that 
x  and  y  should  be  small.  If  the  system  have  more  than  two  inde- 
pendent motions,  we  have  more  than  two  co-ordinates,  but  the 
conditions  of  stability  are  found  in  the  same  way. 

3.  If  x  and  y  be  small  throughout  the  whole  motion  from  the 
one  given  position  to  the  other,  not  only  does  the  system  not 
deviate  far  from  its  undisturbed  course,  but  the  system  at  any 
instant  is  also  very  nearly  coincident  with  its  undisturbed  place 
at  the  same  time.  It  is  important  to  notice  this,  for  the  word 
"  stability  "  is  sometimes  used  in  a  different  sense. 

This  condition  of  stability  may  be  put  under  a  form  in  which 
no  reference  is  made  to  S.  Let  u  and  v  be  the  components  of 
momentum  of  the  system  corresponding  to  the  co-ordinates  6,  </> 

.     ,      .      ,  dT          dT 

respectively,  i.e.  let  M  =  -™p  ,  V  =  "T77>  an°-  let  these  be  expressed 
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as  functions  of  00,  <£0,  6,  <j>  and  t.      Then  the  preceding  equation 
may  be  written 

du         dv 


du         dv 


IV  I    U 

where  a.  and  /3  are  small  arbitrary  quantities.     The  condition  of 
stability  is  that  the  values  of  x  and  y  thus  found  should  be  small. 

4.  As  an  example  let  us  consider  the  case  of  a  projectile. 
If  0  be  the  horizontal,  and  <£  the  vertical  co-ordinate  of  the  par- 
ticle, we  have 


The  equations  to  find  x  and  y  are  evidently 


Hence  the  system  continually  deviates  more  and  more  from  its 
undisturbed  place. 

5.  In  order  to  calculate  the  form  of  S  when  a  system  is  oscil- 
lating about  a  state  of  motion,  it  is  convenient  to  choose  as  co- 
ordinates some  small  quantities  x,  y  which  vanish  in  the  given 
state  of  motion,  Let  the  Lagrangian  function  be  written  in  the 
form 


where  Ln  is  a  homogeneous  function  of  x}  y,  x',  y\    Then  by  a 
theorem  of  Euler's 

dL9     .  dL, 


where  the  2's  imply  summation  for  all  co-ordinates. 
As  in  Art.  9  of  Chap.  IV.  we  have 

d  dL1  _  dL/i 
dt  dx       dx 

and  the  oscillations  are  given  by 

d^dL^dL^ 

dt  dx       dx 
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Hence  we  find 

d  dL9       ,  dL 


Integrating  we  have 


Thus  the  integration  has  been  effected,  but  in  order  to  express 
8  as  a  function  of  a?0,  y0,  x,  y  and  t,  it  will  be  necessary  to  find 
xd>  yd>  x'>  y  m  terms  of  these  quantities. 

6.  As  an  example,  let  the  position  of  the  system  depend  on 
one  co-ordinate  x  and  let 


where  the  coefficients  are  all  constants.     We  then  find  by  the 
process  just  indicated  that  A^  =  0  and 


n  (x* 


2  emt_e-mt 

^ 

where  m2  =  -      .     Applying  the  criterion  of  stability  we  find  that 


,     ,    will  finally  become  small  if  m  is  real.     The  motion  is  there- 

fore unstable  or   stable   according  as  An>  Sn  have  the  same  or 
opposite  signs. 

7.     If  the  position  of  the  system  depend  on  two  co-ordinates 
a?,  y,  let 


+  Cuxx  +  Cytctf  +  Cnyx' 
We  then  find 


where 

'        '  '        ' 
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If  we  now  express  of,  y,  x0',  yQf  in  terms  of  cc,  y,  #0,  yQ  and  t,  we 
find  for  cr  a  fraction  whose  numerator  is  a  homogeneous  quadratic 
function  of  x,  y,  #0,  y0,  the  coefficients  being  linear  functions  of 
exponentials  of  t,  and  whose  denominator  is  another  linear  function 
of  the  same  exponentials.  These  exponentials  become  sines  and 
cosines  when  the  motion  is  stable.  Thus  when  the  given  motion 
is  steady  the  simplest  inspection  of  the  form  of  cr  will  determine 
whether  the  motion  is  stable  or  not. 

Referring  the  motion  to  principal  co-ordinates  for  the  sake  of 
brevity,  and  writing  2  G  =  Clz—  (721,  we  find  that  cr  must  satisfy 
the  differential  equation 

da-         1     (da-           V        1      tdtr  V 

-JT  +  TTTT  I  -T-  +  CrV  1    +  7777-  1  -7 «w  }   = 

dt      ZBn  \dx        >J J       2j522  \dy          J 

This  equation  is  obviously  satisfied  by  such  a  function  as  that 
just  described.  The  solution  of  this  equation  may  be  reduced  to 
linear  equations  and  thus  a  may  be  found.  But  it  is  unnecessary 
to  dwell  on  this,  for  this  would  be  equivalent  to  returning  to  the 
Hamiltonian  equations. 

8.  If  we  wish  to  determine  the  condition  that  the  general 
course  of  a  dynamical  system  is  stable  without  requiring  it  should 
be.  near  its  undisturbed  place  at  any  the  same  time,  it  is  more 
convenient  to  use  the  Characteristic  function.  Suppose  that  the 
Characteristic  function  has  in  Jacobi's  manner  been  expressed  as 
a  function  of  the  co-ordinates  6,  c/>,  the  constant  h  of  vis  viva  and 
two  arbitrary  constants  alt  az.  Then 


The  relation  between  6  and  $,  which  may  be  called  the  equa- 
tion to  the  path  of  the  system,  is  given  by 

dV     . 

-^r  - 

where  6t  is  another  constant.  Let  the  system  be  disturbed  from 
the  same  initial  position  so  that  the  whole  energy  is  unaltered. 
The  change  in  $  corresponding  to  any  given  value  of  6  is  found 
from 


~T  —          i 

da*  da 

Let  A  be  the  initial  value  of  -—  ,  then 
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The  condition  that  the  path  should  be  stable  is  that  the  coeffi- 
cient of  Bat  should  not  be  large. 

We  might  also  use  the  function  called  Q  by  Sir  W.  R  Ha- 
milton, but  it  seems  unnecessary  to  dwell  more  on  this  subject, 

9.  The  instability  of  a  system  may  be  deduced  from  the 
Hamilton  ian  Principal  or  Characteristic  functions,  expressed  as  a 
minimum.  Suppose  a  dynamical  system  to  move  from  one  posi- 

tion A  to  another  B  in  a  time  t,  then  the  motion  may  be  found  by 

rt 
making  the  first  variation  of  8=  I  Ldt  equal  to  zero,  the  time  of 

J  o 

transit  being  constant.  The  constants  of  integration  are  deter- 
mined by  the  conditions  that  the  co-ordinates  have  given  values 
when  t  =  0  and  t  =  t.  To  determine  whether  S  is  a  maximum  or 
minimum  or  neither  we  must  examine  the  second  variation  and 
here  we  have  the  assistance  of  Jacobi's  rule.  The  determination 
of  the  constants  will  depend  on  the  solution  of  equations  and  may 
lead  to  several  different  kinds  of  motion  from  A  to  B.  One  of 
these  will  be  the  actual  motion.  Let  us  move  B  along  this  until 
one  of  the  other  motions  coincides  or  as  we  may  say  approaches 
indefinitely  near  to  this  actual  motion.  We  have  then  reached 
a  boundary  beyond  which  the  integration  must  not  extend  if  8  is 
to  be  a  maximum  or  minimum.  See  Todhunter's  History  of  the 

J1  T 

Calculus  of   Variations,  page-  251.     Further  -      ,  if  0  be  a  co- 


ordinate, is  positive  throughout  the  limits  of  integration,  so  that 
S  will  be  a  minimum  and  not  a  maximum. 

10.  When  there  are  several  co-ordinates  0,  $,  &c.  which  are 
to  be  found  as  functions  of  the  time,  we  may  easily  show  that 
Jacobi's  condition  is  a  necessary  one,  and  this  is  all  that  we  require 
for  the  next  proposition.  If  the  system  can  move  in  two  ways 
from  A  to  B,  then  B8  =  0  along  each,  and  therefore  when  these 
two  are  adjacent  we  have  both  ££  =  0  and  S(S+BS)  =  0.  This 
shows  that  the  second  variation  can  be  made  to  vanish  by  taking 
one  variation  through  the  other.  This  second  variation  will  then 
be  the  same  as  the  quadratic  term  of  the  series  obtained  by 
changing  the  co-ordinates  0,  $  into  6  +  BO,  $  +  S$,  because  we  can 
take  8*0  =  0  and  S2<£  =  0.  Hence  as  the  sum  of  the  terms  of  the 
third  order  does  not,  in  general,  vanish  for  this  displacement,  it  is 
clear  that  S  cannot  be  either  a  maximum  or  a  minimum. 

Let  the  actual  motion  be  from  A  to  B,  and  let  a  neighbouring 
motion  starting  from  A  lead  the  system  to  a  position  C  reached  in 
the  same  time  along  the  actual  motion  before  reaching  B.  Then 
we  can  show  that  a  variation  of  the  actual  motion  from  A  to  B 
can  be  found  which  makes  S*$  of  any  sign.  Let  P  be  any  position 
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on  the  neighbouring  motion  before  reaching  (7,  and  Q,  one  on  the 
actual  motion  after  passing  C.  Then  considering  P  and  Q  as  fixed 
and  also  the  time  of  transit,  the  motion  along  PCQ  cannot 
make  S$=  0;  for  this  condition  is  known  to  lead  to  the  ordinary 
dynamical  equations,  and  it  is  clear  that  (impulsive  forces  being 
set  aside)  no  actual  motion  can  be  discontinuous.  But  there  is 
discontinuity  at  (7,  for  otherwise  when  the  system  is  started  from 
B  towards  A,  two  courses  would  be  open  to  the  system  on  arriving 
at  G.  Hence  the  first  variations  of  $for  an  imaginary  motion  along 
PCQ  are  not  zero,  and  therefore  may  have  any  sign.  But  since 
the  discontinuity  at  C  is  of  the  first  order  of  small  quantities,  this 
first  variation  is  of  the  second  order.  Now  the  value  of  8  for  the 
actual  motion  is  equal  to  that  along  the  neighbouring  motion  to 
C  and  then  along  the  actual  motion  to  B.  Hence,  P  and  Q  being 
still  fixed,  variations  of  the  actual  motion  from  A  to  B  can  be 
found  which  make  &S  of  any  sign. 

11.  Let  us  apply  this  theory  to  determine  the  stability  of  a 
given  state  of  motion.  First  let  us  suppose  the  given  motion  to 
be  steady  and  to  depend  on  only  two  variables.  If  we  use  the 
function  S  there  will  be  one  co-ordinate  and  the  time,  if  V  two 
co-ordinates.  Let  the  system  be  disturbed  at  any  moment  by  an 
alteration  of  the  velocities  of  its  several  parts,  so  that  the  initial 
position  of  the  disturbed  motion  is  an  undisturbed  position.  If 
the  motion  be  stable  the  system  will  oscillate  about  the  un- 
disturbed motion,  the  oscillation  repeating  itself  at  a  constant 
interval.  It  follows  therefore  by  Jacobi's  rule  that  $  or  V  cannot 
be  a  minimum  for  a  period  longer  than  the  time  of  a  half-oscilla- 
tion. If  therefore  8  or  V  be  a  minimum  for  all  variations,  starting 
from  A  and  ending  at  B,  where  B  is  a  position  on  the  steady 
motion  reached  by  the  system  at  an  interval  as  long  as  we  please, 
then  the  motion  is  unstable. 

If  we  give  a  meaning  to  the  word  "  stable"  somewhat  different* 
from  its  usual  signification,  we  may  extend  this  proposition  to 
determine  a  test  of  the  stability  of  any  motion,  whether  steady  or 
not.  All  we  have  assumed  is,  that,  if  the  motion  be  not  altogether 
unstable,  there  are  some  disturbances  which  will  cause  the  system 
periodically  to  assume  the  same  positions  as  it  would  have  done  if 
it  had  been  undisturbed,  but  the  interval  of  these  periods  may  be 
any  whatever  provided  the  first  be  finite.  If  we  use  the  Character- 
istic function,  these  disturbances  must  be  such  as  not  to  ajter  the 
constant  of  vis  viva,  and  if  the  Principal  function,  they  must  be 
such  as  to  bring  the  system  to  an  undisturbed  position  in  the  same 
time. 

*  This  meaning  does  not  always  agree  with  the  results  of  Art.  14,   Chap.  iv. 
[See  also  Arts.  17  and  18,  Chap,  iv.] 
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12.  Next  let  us  suppose  that  as  the  system  proceeds  from  the 
initial  position  A  along  the  actual  motion,  S  ceases  to  be  a 
miirimum  at  some  position  B.  The  conditions  for  a  minimum 
are  of  two  kinds.  Suppose  the  system  to  depend  on  two  co- 
ordinates 6,  </>,  and  let  L  be  the  Lagrangian  function,  then  (1)  we 

d?L       ,  d*LdzL        d*L     ,    A        _  ,    ,ON  .4 

must  have  -j^-*  and  -7^  -rrn  — TZTTTT  both  positive,  and   (2)  it 
civ  cW   a<f>        at)  a<p 

must  be  possible  to  choose  three  arbitrary  constants  which  enter 
into  a  very  complicated  expression,  so  that  this  expression  may 
never  become  infinite  between  the  limits  of  integration.  The  first 
condition  is  clearly  always  satisfied  since  the  vis  viva  of  any 
system  is  necessarily  positive  for  all  values  of  6'  and  <£'.  The 
second  condition  will  fail  if  there  are  two  neighbouring  motions  by 
which  the  system  can  proceed  from  A  to  any  position  between 
A  and  B.  If  this  be  the  mode  of  failure,  it  is  clear  from  the 
reasoning  of  Art.  2  that  the  conditions  of  stability  are  satisfied  for 
one  kind  of  disturbance,  and  that  therefore  some  at  least  of  the 
harmonic  motions  are  stable  or  oscillatory,  though  the  motion  may 
be  unstable  for  a  different  kind  of  disturbance. 


13.  [These  conditions  become  much  simpler  when  the  position 
of  the  system  is  determined  by  one  co-ordinate,  or  when  the 
Lagrangian  function  can  be  reduced  to  depend  on  one  co-ordinate. 
Let  this  co-ordinate  be  so  chosen  that  it  vanishes  along  the  given 
motion,  and  let  us  also  suppose  that  both  it  and  its  differential 
coefficient  with  regard  to  t,  are  small  for  all  neighbouring  con- 
strained motions.  Let  this  co-ordinate  be  called  0  and  let  the 
Lagrangian  function  be 


Then  since  the  Lagrangian  equation  of  motion  is  satisfied  by 
hypothesis  when  0  —  0,  we  have  Al  =  B{  where  the  accent,  as 
usual,  denotes  differentiation  with  regard  to  t. 

If  the  system  be  now  conducted  from  the  initial  position  A  to 
any  other  position  B,  both  on  the  given  motion,  by  any  neighbour- 
ing mode  of  motion,  we  have 


s-  L0dt  +     Baer  +  cne&  +  Aur  dt. 

If  6  =  u  be  any  solution  of  the  Lagrangian  equation 


R.  A. 
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we  may  write  the  function  S  in  the  form* 


The  second  term  is  essentially  positive,  since  Bn  must  be 
positive.  Hence  S  is  a  minimum  along  the  given  motion  unless 

we  can  so  choose  the  arbitrary  displacement  6  as  to  make  —  -  =  0. 

This  gives  Q—CM  where  c  is  some  constant.  But  0  must  vanish 
at  the  two  limits  A  and  B,  hence  this  choice  of  6  is  excluded 
unless  there  is  some  neighbouring  mode  of  motion  by  which  the 
system  could  move  freely  from  the  given  initial  position  A  to  the 
position  B.  The  result  is  that  S  cannot  cease  to  be  a  minimum 
before  the  first  instant  at  which  some  neighbouring  motion  will 
bring  the  system  (starting  from  A)  into  coincidence  with  some 
contemporaneous  position  on  the  given  motion.  If  the  given 
motion  be  steady,  it  follows  that  S  cannot  cease  to  be  a  minimum 
before  a  time  which  is  half  that  of  a  complete  oscillation. 

We  thus  have  a  test  of  stability.  If  the  system  depend  on  one 
co-ordinate  and  if  S  be  a  minimum  when  the  limits  of  integration 
are  from  the  initial  position  A  to  all  positions  on  the  actual  motion, 
that  motion  is  unstable.  But  if  S  cease  to  be  a  minimum  at  some 
point  C,  then  the  actual  motion  is  stable  from  A  to  C.] 

*  [Following  Lagrange's  rule  we  may  write  the  second  term  of  S  in  the  form 
Su6»+(Cn+2\)06'+      ^u+X'    6*    dt. 


The  quantity  outside  the  integral  sign  is  to  be  taken  between  the  given  limits 
and  is  zero,  since  6  vanishes  at  each  limit.    Let  us  now  put 

u 


It  is  clear  that  this  value  of  X  cannot  be  infinite  between  the  limits  of  integra- 
tion unless  u  vanishes.  For  by  hypothesis  u  and  u'  are  both  finite  and  the  co- 
efficients -fiu  and  Cu  in  the  Lagrangian  function  are  also  finite.  It  then  easily 
follows  from  the  equation 


that  (<7n  +  2X)2= 

Hence  the  second  term  of  S  becomes 


which  is  the  result  in  the  text.  This  might  also  have  been  deduced  from  Jacobi's 
general  transformation  with  one  independent  variable  given  in  Prof.  Jellett's 
Calculus  of  Variations  or  Prof.  Price's  Differential  Calculus. 

If  w=0  between  the  limits  of  integration  this  transformation  fails,  but,  as  is 
evident  from  the  argument  in  the  text,  we  choose  0  =  u  to  represent  a  neighbouring 
free  motion  such  that  u  =  0  just  before  the  system  reaches  A  ;  also  B  is  so  placed 
that  the  next  instant  at  which  u  —  0  is  after  the  system  has  passed  B.} 
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14.  [We  shall  conclude  the  chapter  with  the  application  of 
this  criterion  of  stability  to  some  simple  case. 

A  particle  describes  a  circular  orbit  about  a  centre  of  force 
situated  in  the  centre.  It  is  required  to  deduce  the  conditions  of 
stability  as  to  space  from  the  Characteristic  function  V. 

Let  a  be  the  radius  of  the  circle,  n  the  angular  velocity  of  the 
particle  about  the  centre  0.  Let  (/>  (a)  be  the  law  of  force.  Let 
A  and  B  be  two  points  taken  on  the  circular  orbit,  and  let  the 
particle  be  conducted  from  A  to  B  by  some  neighbouring  path 
with  the  same  energy  as  in  the  circular  orbit.  Let  r  =  a  +  p  be 
the  radius  vector  of  this  path,  corresponding  to  any  angle  6. 

If  v  be  the  velocity  at  any  point  of  this  path,  we  easily  find 


f        p_ 

I  •*•  l/\  •       ***    I   IT* 

[        a 


If  s  be  the  arc  of  the  path,  we  have 
ds  1  fdc. 


If  the  angle  A  OB  =  /3,  we  therefore  have 

dO, 


If  the  neighbouring  path  be  a  free  path  described  with  the 
same  energy,  its  equation  is 

p  =  L  sin  pO, 

where  0  is  measured  from  the  radius  vector  OA,  and  L  is  an 
arbitrary  constant.  This  free  path  will  cut  the  circle  again  in 
some  point  C.  If  the  angle  A  OC=  7,  we  have  py  =  IT. 

If  B  coincide  with  C,  we  find  by  substituting  this  value  of  p  in 
the  expression  for  V,  that  the  second  term  of  V  is  zero.  If  B  be 
beyond  C  but  such  that  the  angle  COB  is  less  than  7,  draw  two 
free  paths  one  from  A  as  before  and  the  other  backwards  from  B 
to  meet  the  former  in  some  point  P.  Then  the  angle  A  OP  —  ^/3. 
If  the  particle  be  conducted  from  A  to  P  along  one  path  and 
from  P  to  B  along  the  other,  we  find  that  the  excess  of  the 
action  over  the  action  in  the  circular  arc  AB  is  equal  to 


Since  p/3  is  greater  than  TT  and  less  than  2-7T,  this  excess  is 
negative.     The  action  along  the  circular  arc  is  therefore  not  a 
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minimum  if  B  be  beyond  the  first  intersection  of  a  neighbouring 
free  path. 

Lastly  we  may  show  that  the  action  is  a  minimum  if  B  lie 
between  A  and  C.  To  prove  this  we  write  the  integral  in  the 
expression  for  V  in  the  form 


The  first  of  these  two  terms  is  zero  since  p  vanishes  at  each 
limit.    Following  Lagrange's  rule  we  make 


The  integral  then  becomes 


dd. 

This  is  always  positive  and  the  action  along  the  circular  arc  is 
a  minimum.  The  argument  however  requires  that  X  should  not 
become  infinite  between  the  limits  of  integration.  It  is  easy  to 
see  that 


where  E  is  a  constant  to  be  chosen  at  our  pleasure.    But  if  /3 
exceed  -f     it  is  impossible  to  choose  E  so  as  to  keep  \  finite  be- 

tween the  limits  of  integration.     Hence  the  action  is  a  minimum 
only  if  the  angle  A  OB  subtended  by  the  limiting  positions  at  the 

77* 

centre  is  less  than  -  .     The  circular  orbit  is  therefore  stable  if  p2 

P 
be  positive. 

If,   however,  p2  be  negative,  the  expression  for  \  changes  its 
character.    Writing  —  q2  for  p2,  we  find 


thus  the  value  of  X  can  be  chosen  so  as  not  to  be  infinite  for  all 
positive  values  of  0.  In  this  case  the  function  F  is  a  minimum 
for  all  arcs  AB  however  distant  B  may  be  from  A.  The  circular 
motion  is  therefore  unstable.] 
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